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PREFACE 


In this book I have tried to give a fairly complete 
account of the practical methods of analysing any 
given curve or set of observations into its harmonic 
constituents, if it be a periodic function, and of detect- 
ing any periodic component in it if it is not a periodic 
function. 

The book will, I hope, be readily intelligible, and 
readable, to anyone familiar with the elements of 
trigonometry and the calculus: for this purpose, an 
introductory chapter has been added dealing primarily 
with imaginary quantities, but including also, for 
convenience, a few theorems on the integrals and sums 
of trigonometric expressions that are required later in 
the book. 

I have ventured to use the term “‘ artificial function ”’ 
to signify a function which is defined by one mathe- 
matical expression over one part of its range and by 
another expression over another part of its range. 
This may displease some modern writers on pure 
mathematics, who will contend that there is no 
necessity for such a term. But such functions are of 
such wide occurrence in nearly all branches of mathe- 
matical physics that it seems desirable to introduce 
such a term to ‘cover them; and the difference 
between such functions and analytic ones is very 
obvious to all elementary students and_ practical 
workers. Fourier’s theorem, in fact, may be regarded 
as the method par excellence of representing an 
artificial function in the form (or semblance) of an 
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analytic function. And it is partly because of this, 
as well as its property of analysing functions into their 
harmonic constituents, that it is of such great value 
in mathematics. 

The view taken of Fourier’s theorem in this work 
is rather different from that sometimes adopted. 
There are many known ways of expanding an arbitrary 
function in a definite series of analytic functions which 
expansion is valid for all points between two definite 
limits, and Fourier’s theorem may be looked upon as 
merely one of such cases. In the other cases, how- 
ever, we neither know nor care what the series repre- 
sents outside these limits; while the Fourier series, 
outside these limits, represents continual repetitions of 
its behaviour between the limits. We have therefore 
taken Fourier’s theorem as applying only to periodic 
functions, either artificial or analytic. However any 
function behaves over any finite range, we can, by 
repeating this behaviour an indefinite number of 
times to both minus infinity and plus infinity, obtain 
an artificial periodic function from it. 

One advantage in looking upon Fourier’s theorem 
in this way is that the Fourier’s series (apart from 
the constant term) of any artificial periodic function 
can be directly written down from an examination 
of the discontinuities of the function and its succes- 
sive differential coefficients. Moreover, the Fourier’s 
series representing an analytic periodic function can 
be deduced in the same manner as the limit of the 
case when the discontinuities are infinitesimal; 
although this is rarely the easiest way of obtaining 
the expansion. For instance, the harmonic analysis 
of the function 2cos?z can be deduced from the 
analysis of the artificial periodic function of period 27, 
and represented from 0 to 2m by 2 cos? gt, by putting 
q=1. This artificial function possesses, both for 
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itself and for all its differential coefficients, discon- 
tinuities when t= 0; and from the formule on 
p- 59 the Fourier constants are readily found to be 
given by 
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Putting g = 1 gives zero for all values of » save when 
nm = 2 for which the expressions become indeterminate. 
If in this case we put g = 1 + 6 where 6 is very small, 
me castly sce that a, = 1 and 5, =°0, and thus that 
the Fourier’s series representing 2 cos? ¢ is cos 2¢ plus 
a constant. 

Chapter III deals with the harmonic analysis of 
various simple artificial functions, for the most part 
built up of straight lines, parabolas or sine curves. 
This chapter will afford the student excellent practice 
in the evaluation of trigonometric integrals. I have 
been surprised at the very small proportion of students 
who can evaluate such integrals correctly after having 
been through an ordinary course of the Integral 
Calculus. This chapter accordingly will form a useful 
item in the student’s education. In fact, the whole 
book has been written with the intention of making 
it constitute a valuable chapter in the education of 
any student of applied mathematics, as well as of 
making it valuable to the actual worker engaged in 
any kind of harmonic analysis. 

Chapter V deals with the harmonic analysis of any 
periodic function defined by the co-ordinates of a 
number of points and will be, for many readers, the 
chapter of most practical importance in the book. 
I have tried to include therein all methods that I 
think of practical value, while giving full attention to 
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the most satisfactory ones. Attention may be drawn 
to what is, I think, a new method of harmonic analysis, 
which consists in joining the given points by arcs of 
parabolas of the second or third degree and then 
obtaining the exact harmonic analysis of the artificial 
function defined by this series of arcs from the dis- 
continuities. 

Chapter VI, on theoretical considerations, I have 
added for the sake of completeness, and also because 
I find that most students of physics and engineering 
do show a considerable interest in the proofs of the 
various mathematical methods and operations taught 
them. Very few, in the course of their education, 
simply want to be given mathematical results that they 
can use and betray no interest in the establishment of 
the results. 

In Chapter VII I deal with Fourier’s integral 
theorem, the chief application of which is in the 
science of Optics; but since it is the limiting case of 
Fourier’s theorem when the period is made infinite, 
it could not logically be excluded, and it is, more- 
over, wanted for the next chapter. 

Chapter VIII deals with the detection and deter- 
mination of any periodic components in any non- 
periodic function of which observations have been 
made over some finite range of time. ‘The method 
adopted has been to apply the theorem of Chapter VII 
to the artificial function defined to be equal to the 
observations over the region for which these extend, 
and to be zero for all values outside this region. 

The detection of periodic components in pheno- 
mena not themselves strictly periodic is required 
in analysing many meteorological and astronomical 
phenomena. 

Finally, a very brief historical survey of the subject 
has been added in Chapter IX. 
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I need hardly say that I shall be grateful to any 


reader who will be kind enough to point out any 
mistakes or misprints. I shall also be grateful for 
any suggestions, especially from experienced practical 
workers, for incorporation in a further edition of the 
book should one be called for. 
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The University, 
Manchester. 
May, 1925. 
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PetwACTICAL. CL REATISHZON 
FOURIER’S THEOREM AND 
HARMONIC ANALYSIS 


INTRODUCTION 
ON IMAGINARY QUANTITIES 


Berore proceeding with the proper subject-matter 
of this book, it is desirable, we think, to give a very 
brief account of the so-called, and very badly called, 
“imaginary quantities,” as we shall frequently have 
to make use of them. It is almost impossible, in 
mathematical work, to deal with sines and cosines 
without imaginary quantities inevitably presenting 
themselves; so it is very important at the outset to 
get clear ideas about them and about their relations 
to the trigonometrical functions. 

Let O and A represent any two points in a plane and 
letay fepresent the ~° vector”? OA, 2.¢, Virepresents a 
step of magnitude OA in the direction OA and in the 
sense from O to A. This being so, it is evident that 
2V will represent a step of magnitude equal to twice 
the distance from O to A in the direction from O to 
A; and that in general if a is a positive quantity aV 
will represent a step of magnitude a times OA in the 
same direction and sense. Further, —V will clearly, 
with the ordinary meaning of the minus sign, represent 
something which will cancel or neutralise V, and 
therefore represents a step of magnitude equal to OA 
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in the sense from A to O; and similarly — aV repre- 
sents a step of a times OA in the direction and sense 
of AO. Hence the effect of the factor (— 1) is to 
reverse the direction of a vector without altering its 
magnitude. But reversing the direction of a vector 
is equivalent to rotating it through two right angles 
or a radians. Similarly, (— 1) times — V represents 
— V rotated through z radians or V rotated through 
2x radians. Proceeding in this manner, we see that 
if m is a positive integer (— 1)” V represents V turned 
through an angle of mm radians. Hence we may say 
briefly that the factor (— 1)” represents a rotation of 
nm radians. 

Now, exactly as in elementary algebra when dealing 
with indices, we define them first for positive integral 
exponents and then find we must interpret them for 
fractional or negative values because our mathematical 
operations give rise to them, so we must do the same 
here. If(— 1)”+” is always to represent a rotation 
of (m+ n) m radians whenever m+ n is an integer 
(— 1)” must represent a rotation through ma radians, 
whether is positive or negative, integral or fractional. 

Hence we take (— 1)” as representing a rotation of 
nx, to be universally true. Writing mx = @ or n= 6/z, 
we see that the factor (— 1)9/* applied to any vector 
rotates it through an angle 8. In particular, the factor 


(— 1)! or ,/—1, which is generally denoted by 3, 
: M4 . . 
represents a rotation of —, 7.¢. rotation through a 
2 
right angle. If we denote the vector V turned 


through an angle 6 by Vs, we have the result that, in 
general, 


Vo = (— 1% V, 
and in particular, 
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We will now obtain another expression for V» in 
terms of V. Let OA in Fig. 1 denote the vector V 
and OB denote Vs. Draw BC perpendicular to OA, 
then OC denotes the vector V cos @ or (cos 6).V. 
Similarly, CB is a vector whose magnitude is equal to 
that of (sin @).V, but since the vector CB has been 
rotated through a right angle away from the direction 
of V it is represented as a vector byisin@.V. But the 
vector OB is equal to the vector OC plus the vector 


7B 


CB, which is merely asserting that the step OB is equal 
to the step OC plus the step CB. Hence we have 


Ve=cos@.V+2 sin @0.V 
or Ve = (cos 6 + 7 sin 0).V. 
But we had previously found that 
Vo={(— 1)% V, 
hence we see that the two expressions (— 1)%/" and 
cos @ + i sin # are equal and both represent rotation 


through an angle 6. If nis any integer and we replace 
6 by n@, we have 


noé 
(— 1)* = cos 79 + 7 sin 10, 


no 
but (— 1)* = {(—1)%"}" = {cos 0 + 2 sin 0}". 
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We thus arrive at the result which is known in 
trigonometry as De Moivre’s Theorem, viz. that 


(cos 8 + 7 sin 0)” = cos nO + 7 sin 6. 


De Moivre’s Theorem is generally regarded as belong- 
ing to the advanced part of trigonometry, yet, as we 
have seen, it is nothing more than the expression of 
the fact that a rotation of 78 is equal to n rotations of 8. 

We can, by De Moivre’s theorem, easily obtain the 
power series for the sine and cosine. If we write 
6 = x/n we get 


o 6 x coo. BONE 
cos ¥ + 21s1n x =[ cos — + 2 sin ) 3 
n n 
Now let 2 tend to infinity, while » is finite; so that 


cos —~+ 7 sin ~ represents a rotation through an 
n n 


infinitesimal angle ~, and thus reduces to I + 7— in 
n n 


the limit. We then have 
cos # + 2 sin # =( + ix)’ pen (6) 
nN 


when is infinite. Expanding this by the Binomial 
‘Theorem, we obtain 


Sto : Te\. 
COTY SI = T + ix + (x — Dine + 
n 


[2 
(Dre MDE =D 
3 4 


when 1 is infinite. 

If we remember that 7? = — 13; 13 = — i, etc., and 
that, when this is done, the terms devoid of the factor 7 
represent displacements in the line of V; and that 
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those which contain the factor 7 represent displace- 
ments at right angles to V, we see that 


cos ¥ = I -(: = Jiee = Nae 


and 


sing =x —(1-1)(-2 AG a 5 
n n} |\3 


when 7 is infinite. 

Now both these series, when m is very large, are 
such that for any given value of x a finite number of 
terms will give their sums to any desired degree of 
accuracy. Hence we may safely put n= © orI/n=o0 
in these series. For in any term of finite order there 


are only a finite number of the factors (: = *) Pir 
n 


( — *), etc., and these give exactly unity when 1/n 
n 


is exactly zero. Hence we have the important results 
that 


6 oe. en 22) 
and 
iB (3) 


It is shown in books on Algebra that the limit, when 
nm is infinite, of the expression 


e) 


ee, where ¢ = 2°71821 ... and is the base of the 
Naperian logarithms, hence we see by (1) that 
CGce tr Paresiie View Oa on md) 


This is consistent with the previous value, viz. 
(— 1)*/", found for the left-hand side since it only 
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implies that (— 1) = ¢"=cos a+ sin a. If we 
change 8 into — @ in the equation cos 8 + Sin == nee, 
we get COSut== 4 SING =p Oe 


adding and subtracting these equations we have 


(QOx ei = 
2 
and 
: 29 — p—i0 
Since se = 
21 


These results are known as the exponential values 
of the sine and cosine. Their meaning can easily be 


D 


Fic. 2. 


erasped from the diagram, Pig. 2. “Dhevfactonga™ 
since it is the same as cos # + 7 sin @, represents a rota- 
tion through an angle 8. Hence in Fig. 2 OA, OB and 
OC represent respectively V, eV and 2— Vs. Clearly, 
the sum of the latter two vectors is a vector equal 
to 2 cos @ times V. Similarly OD) represents) the 
vector — ¢—* V and the resultant-of OB and OD is 
clearly a vector of magnitude equal to 2sin@ V, but at 
right angles to it, and hence represented by 2 7 sin @ V; 
this is equivalent to the equation e — ¢— = 27 sin 0. 


Any expression of the form a + d7 is called a Complex 
Quantity; a is said to be its “real” part and bz, or 
briefly b, is said to be its “imaginary” part. Much 
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a 
better and more appropriate names would be “ axial ” 
part and “non-axial” part. From what has been 
said above it is clear that (a + b7) Vrepresents a vector 
with a component aV along V and a component bV 
at right angles to V._ It thus represents a vector equal 
in magnitude to »/a* + b? V, and making an angle 
tan—+*b/a with V. But this vector can be represented 
by 2°,/a? + b? V, where 6 = tan—! b/a. 

Hence we see that 

a+ bi = re’?; where r = /a?+ 8 
and § —-tan—" b/a; a result that can be obtained 
immediately by writing a + bi in the form 
i a Die 
I Tae Jah oy 

The quantity r is called the “ modulus” of the 
complex quantity a+ i, and the quantity 9@ its 
“ argument.” "The modulus of any quantity Q is often 
denoted by | Q|; thus if Q is a real negative quantity 
| Q| denotes its numerical magnitude taken positively. 


Some TRIGONOMETRICAL THEOREMS 


It is convenient to place here some results which 
we shall require later so as to prevent the argument 
then from being interrupted by lengthy digressions. 

The formule for the resolution of the product of 
two sine or cosine terms into the sum of two terms are 
so frequently required that we will first state them 
here for reference. They are 

cos A cos B = $cos (A + B) + ¥cos (A — B), 
sin A sin B = $cos(A — B) — ¥cos (A + B) 
and sin AcosB = $sin(A + B) + $sin (A — B). 


If m is any integer we have the well-known results that 


2m 2m 
\ cos wxdx = 07%. (A) and | Gil VAN ==) Oboe CA) 
ce) 


° 
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sin 1x 
For the indefinite integrals are respectively ———~— 


cos nx 


ehate le = , and both these expressions take the same 


values when x = 22 as when x = o. From these results 
we can easily deduce that, if m is also an integer, 


Phi. 
\ cos mx cos nxdx = 0 if m + n and 


i =nifm=n .. . . (B) 
27 
\ sin mx sin nxdx =o if m = n and 
a ee ee LO 
2 
and \ Sin MecOS 7vdx == © invallicasesi.m ee) 
A 


To prove the first, we write cos mx cos nx = 
% cos (m+ n)x +4 cos (m— n)x; and, since 
m + n and m — n are both integers, the result follows 
at once from (A) unless m7. ‘Im thisecase, thes 
last term is $ and the integral of this from o to 
2n isa. ‘The result (C) follows in an exactly similar 
manner. For (D) we write sin mx cosnx = $sin(m 
+ n)x + % sin (m — n)x, from which it is seen by” 
(A) that this is zero in all cases, even when m= 7; 
for in this case the second term vanishes. 


Some summation results, of which the integral 
results just given are the limiting case, will now be 
given. If p is any integer we have 


- 


cos —~ + cos 42 4208 = 4, casero 
2 P p 


or as it is professionally written, 


p 27g 


COS 
q=t 


SO Gh oa a ere ey 
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The ¢g under the & indicates which letter the 
summation is with respect to, and the 1 below and 
the p above indicate the limits of g. This result is 


obvious geometrically, for ~ is the external angle of a 


regular polygon of p sides; so that the sum in question 

represents the sum of the projections of the sides, each 

of unit length, of a closed regular polygon on its base, 

which we have taken as the zeroth or pth side. 
Similarly, 


EN eas 
Seine owe ee (BE) 
q=1 
for this is the projection of the same sides on a line 
at right angles to the base. 
Further, if ~ is any integer, 


q=lI q=1 


unless 7 is a multiple of p, when the sum is pg. For if 
n is prime to 9, the first sum still represents the sum 
of the projections of the sides of the same polygon on 
its base line, but taken in the order mth, 2uth, 3nth, 
etc., and in this sum each side is included once. If 
n and p have a common factor 7, it can easily be seen 
that this sum will mot include all the sides of the 
polygon, but that every 7th side will be taken r times. 
But the p/r sides that are taken would, by themselves, 
form a closed regular polygon of p/r sides, so that the 
sum of their projections is still zero. The only 
exceptional case is when m is a multiple of p, then 
each term in the summation represents the projec- 
tion of the base and therefore the sum is # in the 
first case and zero in the second. 

We can now prove that if m and are integers not 


greater than p/2 both 
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p 
Dia pe a esta ME ah (H) 
paeb es p 
P 

and y sme! sn =o, (I) 
q=1 


if m + n, while both sums are equal to p/2 if m = 1, 
save in the exceptional case when m= n= p/2, 
when the sum in (H) is p and the sum in (I) is zero. 
These results can readily be obtained by replacing 
the products by the sum or difference of two cosines 


of the form cos an(m + )4 m)q 


The restriction on m 


and m is necessary to prevent m+ n becoming equal 
to p or a multiple of it when m + n. 
Similarly, we can prove that with the same limita- 
tions on m and n 
Po. 2nm 270 
> sin phe! cos palsies 


a , uence a) 


in all cases. 


Lastly for this Introduction the definite integral 


sin bx : : ‘ : 
(2? dx is very important in the theory of Fourier’s 


4 x 


series. It may be evaluated as follows : 


. . . cow 
The indefinite integral \ee= of course, 
¢ 


whether c be real or complex. Let us write c= —a 
+ 1b, where a and b are positive, then we have, 
oF : ee (— a+ ib)x a 
- le— (cos bx + 71 s8in bx)dx = a = pi 
\e ( ) [5]: —a+tib 
(e) 


a+ ib a+ ib 


7 @4+ib(a—ab) att BF 
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Equating real and imaginary parts we have, 


\e-# Poser de == oe 

fe) a* +. 5? , 

and \e- % sin bx dx = ee. 
fe) a’? + 5? 


In the first result, the integrand is a “ function ” of 
b and the given result holds for all values of b. Let us 
suppose, then, that we write down this result for 
an infinite number of values of b between o and B 
each differing by an equal infinitesimal amount db 
from the preceding value and add the whole set of 
equations after multiplying by db. ‘The expression 

B 


under the integral sign then becomes jena cos bx db, 


[e) 


which = aes and the right-hand side becomes 


B 
\ : adb_ Which = ‘tan ~* Bia. 


Hence we have 


(oe) 
e—% sin Bx 
\ dx = tan = 2:3 /a: 
- x 


If now we make a = o, the integral is still conver- 
gent and the result becomes 


co 
ye By 


oO 


IU 
dx = tan-!l0 = 


Example.—By integrating both sides of this result 
with respect to B from o to 2b obtain the result 


2 ah 
$1n. ww IU 
\ ; dx = —b. 
ie 2 


CHAPTER I 
PERIODIC FUNCTIONS IN GENERAL 


Classes of Functions. Mathematical functions 
may, for the purposes of the practical mathematician, 
be divided into three classes, viz. (a2) What are known 
as “ Analytic Functions ” which we may perhaps call 
‘Natural Functions 2 such as) 47,)Sin. v,0cm loc 
etc., and combinations of them. ‘These functions exist, 
not only for all real values of x, but also for all complex 
values of x of the form x + 1y.* 

(b) What we may call “Artificial Functions.” These 
are functions which are arbitrarily defined as being 
equal to one analytic function in one part of their 
range and equal to another one in some other part. 
For instance, the function defined by 


f(*) =o whenx< —1 
f(x) =hbwhen—-1<«x*<+1 
and f(*) = o when « > + 1 


is an artificial function which is hereby defined for all 
real values of x.+ Artificial functions must only be 
taken as existing in the regions for which they have 
been defined ; and thus must never be taken as existing 
for complex values of « unless they have been specific- 


ally defined for such complex values. (c) The third 


* Some complicated analytic functions only exist over definite 
regions of the xy plane, but this does not affect the argument. 

+ Strictly speaking, we have not defined the function for the values 
x= -+1. We could have got over this by saying f(x) = 4 when 
—1I<«<-+1, where < means is “less than or equal to.” For 
practical purposes, however, this is an unnecessary bit of hair splitting. 

iZ 
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class are known as “‘ Empirical Functions’; these are 
functions which we arrive at as the result of some 
experiments or observations where we obtain the 
result in the form of a graph or a set of points. In 
either case, the ideal function which we had in mind, 
and which it was the object of the experiment to deter- 
mine, is only known to a limited accuracy owing to 
necessary experimental errors. We are supposed not 
to know the “natural” or “ artificial ”’ function which 
represents this ideal function, as otherwise we should, 


Fic? 3. 


of course, use this representation, instead of the experi- 
mentally obtained values, in any mathematical opera- 
tions we had to perform on it. 

It is almost entirely to artificial and empirical func- 
tions to which we desire to apply Fourier’s Theorem 
and the method of Harmonic Analysis. 

Discontinuities of Artificial Functions. The 
most important distinction between an artificial func- 
tion and an analytic one, is the possession by the former 
of what are called ““discontinuities,”’ the nature of which 
we must now briefly study. The heavy broken line in 
Fig. 3 represents an artificial function which possesses 
discontinuities at the points where « = a, b, ¢ or d. 
At « = a there is an abrupt change in the magnitude 
of the function represented by A, A, ; at x = d, there 
is an abrupt change in slope but no change in magni- 
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tude, while at C there is an abrupt change in the rate 
of change of slope, or curvature, but no change either 
in magnitude or slope. Thus at any point of discon- 
tinuity of a function there is an abrupt change, either 
in the magnitude of the function, or of some of its 
differential coefficients. In general, if there is an 
abrupt change in the magnitude of any of these quan- 
tities, there is also an abrupt change in the magnitude 
of all the higher differential coefficients at the same 
points unless these are zero. It is accidental, so to 
speak, if a function changes in magnitude, but not in 
slapesras at va nig a3 

The discontinuities are measured by the magnitude 
of the abrupt changes at them. ‘The ordinate of the 
point A, is denoted by f(a — o) and that of A, by 
f(a + 0) and the magnitude of the change, f(a + 0) — 
(4 — 0), we shall denote by I, so that Ig stands for 
the magnitude of the abrupt change in f(#) at « = a. 
Similarly the slope of the arc A,B at B is denoted by 
f'(6 — 0) and the slope of BC at B is denoted by 
f(6+ 0), thus f(b + 0) — f(b — 0) is the magnitude 
of the discontinuity of f(«) at « = b, and we shall 
denote it by I’. Similarly, the magnitude of the 
discontinuity of f’(%) at « = will be denoted by 
I’, and so on. 

Periodic Functions. In the great majority of 
practical cases the independent variable in these 
functions is either the time or an angular magnitude. 
It will thus be more convenient to denote it by ¢ than 
by x. A function f(2) is said to be a periodic function 
when a constant J exists, such that f(T + #) = f(2) 
for all values of t. Changing z¢ into J + t we have 
f(2T + ) = f(T + 4) = fd); and in general f(nT + 2) 
= f(z) where is an integer. The least value of T 
for which the equation f(T + #) = f(#) is always 
satisfied is called the period of the function. An 
multiple of the period is said to be a period of the 
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r 


function. The behaviour of a periodic function over 
the range from t= — o to t=-+o is thus a 
continual repetition of its behaviour over the range 
fom-o-to J, 

It is generally convenient, in discussing periodic 
functions, to change the scale of the ¢ axis so that the 
period is 27 as this will be found to make a great simpli- 
fication in all the expressions. Actually, of course, in 
a diagram, we should mark the period divided into 
360° and not into 2” radians. Sometimes it is 
advisable to take an aliquot part of 2m as the period. 
For instance, the length of the radius vector from the 
centre of an ellipse has an angular period of z, and it is 
better to keep it such than to change the variable. 
But we shall suppose in future that the period of a 
periodic function is always 27 unless otherwise expressed. 

The student must grasp clearly the difference 
between a natural (or analytic) periodic function and 
an artificial one. If a periodic function of period 
2m is represented by one and the same analytic expres- 
sion over a whole range of 27, the function will never- 
theless be an artificial periodic function, unless the 
same analytic expression also represents it for all 
values (to + «) outside this range. For instance the 
periodic function, of period 2m, represented from 
t= 0 to t = 2m by e™ is an artificial periodic function 
with discontinuities at t= 0, + 27, + 4m, etc.: so 
too is a function of the same period and represented 
from — x to + 2 by sin pt unless p 15 an integer ; for 
it is only in this case that the periodic function we 
have defined is also represented outside the range from 
—zxn to +a by sin pz. As another illustration, the 
function of period 2m defined by 


I 


Ki) = 


2 + cost 


is an, analytic periodic function ; and so, it will easily 
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be seen, is the function defined by any rational 
expression which only involves sines or cosines of 
integral multiples of ¢. 

Four Special Classes of Periodic Functions. 
Besides the general periodic functions, which have no 
symmetry, there are four important special classes 


Fic. 4. 


characterised by different degrees of symmetry which 
enables them to be easily distinguished from one 
another by the eye. We have: 

(a) Class I called “ sine-harmonic functions.” A 
periodic function belongs to this class when it is 


oT -T ie) vis 27 3a 
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possible to choose the origin so that f(z) = — f( —2) for 
all values of t. ‘The function is then an odd function. 

Writing a + ¢ for tin the equation f(t) = — f( —d), 
we have fa + 2) = —f(—a— 2), = — fa — 2d, 
since 2% is a period, and hence the function is also an 
odd function when the origin is shifted half a period. 
An example of this class of periodic function is illus- 
trated in Fig. 4. 
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’ 
(>) Class II, called ‘ cosine-harmonic functions.” 
A periodic function belongs to this class when it is 
possible to choose the origin, so that f(t) = + f(— 2) 


for all values of t. "The function is then an even func- 


Fic. 6. 


tion. Writing x + ¢ for ¢ in this equation, we have 
f@+ t) =f(—2x—- 4, = fa — 29), since 2z is a period ; 
and hence the function is also an even function when 
the origin is shifted half a period. An example is 
illustrated in Fig. 5. 


0 as 0 ZAS 
(=) (4) 
am as 
fe) 7 27 O T 
(¢) (2) 
Fic. 7 


(c) Class ITI, called “‘ odd-harmonic functions.” A 
periodic function belongs to this class when f(t + 2) 
= — f(t) for all values of t, when = is half the period. 
This relation is independent of the origin. An 
example is illustrated in Fig. 6. 

These three classes of functions can be obtained by 

c 
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repeating any given half wave in different manners, 
see Fig. 7, where (a) is supposed to represent a half 
erlod wave. 

If the second half wave be obtained by rotating the 
first half through 180° about the point 2, we obtain a 
Class I function as in Fig. 7 (6) ; while if it is repeated 
by reversing it fore and aft as in Fig. 7 (c) we have a 
Class II function. If the repetition is obtained by 
rotating the half wave about the axis of #, and then 
displacing it half a period along this axis, we have a 
Class III function, as in Fig. 7 (d). A mere repetition 
of the half wave, without any turning about at all, 


Fic. 8. 


would, of course, result in the period being no longer 
an, but x. 

(d) Class IV, called“ bi-symmetric functions.” These 
are the periodic functions with the highest possible 
degree of symmetry. Each quarter period is geo- 
metrically similar to every other quarter period. An 
example is illustrated in Fig. 8. It will be seen that 
the functions of this class belong to all the three 
previous classes. With the origin as shown the func- 
tion is an odd function, if the origin were shifted to 


ff 8 . 
the function would clearly be an even function; 


moreover, the function is also obviously a Class III 
function. ‘The student should prove to himself that 
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any function that belongs to any two of the first three 
classes necessarily belongs to the other one also, and 
is a Class IV function. 

It is important to be able to recognise to which of 
these four classes any periodic curve that occurs in any 
problem belongs. Let us consider, for instance, the 
open circuit voltage wave form of an alternator. 
Since in practice a “north” field pole is made 
geometrically similar to the adjacent “south” one, 
and since both are symmetrical about their centre 
lines, it is clear that the voltage wave form must be a 
Class IV function. This is not true, however, when 
the machine is on load, for then, as is well known, 
the field due to the armature current weakens that 
due to the field poles on the sides at which the armature 
is approaching the poles and strengthens it on the 
leaving sides. Nevertheless, it is clear that after half 
a period, the induced voltage will still be numerically 
equal, but of opposite sign, and thus that the voltage 
wave form is a Class III function. 


Example 1.—Prove that any periodic function can be represented 
as the sum of an odd and an even periodic function. 

Rose (=O Se ) , fl tae identically. 

The first term merely changes sign on changing ¢ into — ¢ and the 
second one is unaltered thereby. Hence the first term represents a 
Class I function and the second a Class II one. 


Example 2.—Prove that any periodic function can be represented 
as the sum of a Class III function and another of half the period. 


Pe eee) ake 


The first term merely changes sign on changing ¢ into t+ z, since 
2m is the period. The second term is unaltered by this change and 
hence x must be a period of this term. 


Example 3.—Prove that any Class III function can be represented 
as the sum of two Class IV functions. 
For if the function be represented as the sum of an odd and an 
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even function by Ex. 1 it will be seen that both these functions are 
still Class III functions and must hence be Class IV functions. 


Example 4.—(a) Prove that ag-++ a, cos t-+ a, cos 2t-++ dg cos 3¢ 
._. Tepresents a Class II function of period 27. 
@) Prove that 4,sin t+ dysin 2¢-+ dg sin 3t-++ . . . represents 
a Class I function of period a 


(c) Prove that a, cos t+ a, cos 3t-+ a, cos 5¢-+ . . . +d, sin¢ 
+ b, sin 3t¢-+ b, sin St-+ . . . represents a Class III function of 
period 27. 


Each of these series is supposed convergent. 


Example 5.—Prove that either of the series 
(2) 4, cos t+ ag cos 3¢-+ as cos St-+ . . . 
and (6) b,sin t+ bg sin aap b,sin 5t-+ .. . represents a Class IV 
function eI period ae 
For the first series represents a function belonging to Classes II and 
III and the second a function belonging to Classes I and III. 


If the origin in the second series be shifted to the point t = s by 


‘writing t= — = t’ the series becomes 


b, sin (20) bia 3(° ANE ps sin (= aa ') + ae 
which 

. = b, cos t’ — by cos 3t/ + bs cos 5¢7— . . . 
showing that both these series represent the same kind of functions 
referred to different origins. 


Example 6.—Sketch the curves given by 


(a) f(t) = cos t + *2 cos 22. 

(6) f() = sin t+ °2 sin 2¢. 

(c) f@ = sin t+ ‘1 sin 32. 

(2) f(t) = sin t — ‘I sin 32. 

(e) f@ = sin t-+ ‘I cos 3¢. 
and verify in each case that your sketch belongs to the correct class 
of function. 


CHAPTER II 


FOURIER’S THEOREM 


Statement of Theorem. Fourier’s Theorem 
states that any periodic function, whether “ natural ” 
or “ artificial,” of period 2 a can be expressed by the 
series— 

* 
Ke = 4 + a, cost -+ a, cos 2t + a, cos 3¢ + 
+ b,sint + b, sin 2t + 6, sin 3¢+ 


(1) 
where 
1 (2% 
an = — \ Ff (2) cos aH 
% Jo 
pe (2) 
and bn =~ Ff (2) sin nt dt} 
w fe) 


There are a few unimportant restrictions on the 
behaviour of f(¢) which we will notice later in order 
that the theorem should be true. 

We cannot at this stage prove Fourier’s theorem, 
for this requires the proof that the representation Of 
the function by the series given in (1) is possible; 
but we can easily prove that 7f f(t) can be represented 
by the given series the coefficients must have the values 


given by (2). ‘The two sides of (1), being supposed 


as a. . 
* The reason for writing —° instead of a, is so that the value of a, 
2 


will be correctly given by putting =o in the expression for ay. 
21 
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identical for all values of ¢, must remain identical 
when the same operation is performed on both sides ; 
unless this renders the series divergent. Let us 
multiply both sides of (1) by dz and integrate from 
P= O° tO) t= 20. eeby. equations, (A),»pa7.cal stae 
terms on the right-hand side vanish except the first, 
so we are left with 


2 


Wer a 2a dt = na, 


2 
fe) 


which is the value obtained by putting » = 0 in (2). 

Again multiplying both sides of (1) by cos nt dt and 
integrating from 0 to 2x we see that all the terms on 
the right-hand side vanish by (B) and (D), p. 8, 


except the a, cos mt term; thus we get 
27 2m 

\ WE)COS HE di an cos? nt dt 
O 10) 


2m 

a a\ I + cos 2nt)dt = man, 
2 ° 

which is the value of ay given by (2). Similarly, we 

have 


27 2m @ 
\ F(Z) sin nt dt = dy \ sin? nt dt 


° oO 


by, 2 

= — \ (1 — cos 2nt)dt = aby, 
2 te) 

also in agreement with (2). Hence the theorem is 

proved on the assumption that the expansion is 

possible. 

In the expansion (1) the various terms are called the 
‘““harmonics” of f(z), the cosine terms being called 
the “‘ cosine harmonics ” and the sine terms the “ sine 
harmonics.” ‘Thus a, cos ¢ is called the “ first cosine 


ye 


| 
| 
| 
| 


| 
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harmonic” and 2d, sin 2¢ is called the “second sine 
harmonic.” Both the zth sine harmonic and the zth 
cosine harmonic, of course, go through n complete 
cycles in the period of a functions. "The two ath 
harmonic terms a, cos nt + b, sin nt can be combined 
together in the form 


r/ an? + bn? cos (nt — ¢), 


where tan ? = by/an. 

This expression is called the resultant or complete 
mth harmonic; »/an? + 5,2 is its amplitude and — ¢ 
its initial phase. The “‘ first harmonic ” is often called 
instead, the “‘ fundamental harmonic.” The a’s and 
b’s in (1) are often called the “ Fourier constants” of 


I). 


Two or three proofs of Fourier’s Theorem will be given later. The 
student familiar with the theory of linear differential equations with 
constant coefficients will be interested in the following. If f(#) has a 
period 2m we have f(t+ 2n)= f(t). By the symbolical form of 'Tay- 


lor’s theorem f(t-+ 2m) = e?"P f(t) where D =“, ; hence (e?™D_1) . f(2) 
= 0, which is a differential equation for f(t). We know that if 
m is a root of ¢2™ — 1 =o, the solution consists of the sum of a 
number of terms of the or For two imaginary roots, - im, 
the terms take the form A cos™mt- B sin mt, A, B and C being arbi- 
trary constants. Now e?™P — 1=0 has obviously no real root save 
D = 0; to search for imaginary roots we put D = 7A and get e#™A = 1 
or cos 2mA=1 and sinz7mA=o simultaneously; from which it is 
obvious that 4 can be any integer. Hence we see that f(#) can be 
expressed in the form 
f@) = C+ A, cos t+ A, cos 2¢-+ Ag cos 3¢ + 
+ B, sint + B, sin 2t + Bg sin3¢ + . 

The constants must be chosen to satisfy the initial conditions, which 

in this case consist of the behaviour of the function from o to 27, 


>) 


Simplification of Fourier’s Theorem for the 
Four Special Classes of Functions. If f(t) is 
an odd function, so that f(t) = — f(— 4) for all 


values of ¢, we have 


C2 
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fo 4 4, cos t + a, cos 2¢ + ee SEEDING Ee els 
2 
= — {22 + a, cos t+ ay cos 24+. ae 
2 
— >}, sint— b,sin 2i—... 
since cos(— #)=cosz and sin( — #)= — sind. 
This reduces to 


A, 2ay COS =- 20, COS 2b eee 
for all values of t; hence all the a’s are zero, so that any 
odd periodic function can be expanded in the series 
Ad) = b, sint + bz sin 2+ bgsin3¢+ . ~. (3) 
Similarly, if f(t) = + f(— 2) for all values of #, so 
that f(#) is an even function, we find that all the 3’s 
are zero and that 


f(t) = $2 + a, cos t + ay cos 2t + a) 


Again, if f(#) satisfies the condition that f(t + 2) = 
— f(2) for all values of ¢, we have 


a 
“oye ES aa BOS 


b,sint+ bsin2t+ ... 


a 


= — “9 — ay cost + m) — ay cos 2(¢ + 2) — 


Gg COS 3(¢ 3) 
— db, sin(t +2) — db sin 20¢ +a) — bs sin3(¢+2a).. . 
= — “2+ ay cos t — ay cos 2¢ + ag COSe3 t= 
+ b, sin t= Db, sin 2t 4d, sin 37 — . . 
and hence 
a, + 2d, cos 2t + 2a, cos 4t + 2a, cos 6t + 
++ 2b, sin 2t + 2b, sin 4¢ + 2b, sin 6¢ + 


= 0 for all values of t, whence 
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Cea Gd. = O'and Vy —=10, == b,c. == 0; 
and the function can therefore be represented by 
T(t )= a, cos t + as cos 3¢ + 
Cast Veesii feet. 4 8 NCS) 

It will be seen that these three cases correspond to 
Class I, Class II and Class III functions respectively, 
and the reason for the names “ sine-harmonic,” “ cosine 
harmonic ” and “ odd-harmonic ” will now be evident. 
We will now give simplified expressions for the 
Fourier constants in these cases. ‘The expressions 
previously given for the Fourier constants, viz. : 


I 2 I 2m 
an = \ (4) cos mt dt and by, = al f(t) sin nt dt 
fe} 


TJ) o 


can be written 
Ree Lue 
Cae a f() cos nt dt and b, = ‘ f(#) sin nt dt ; 
7% maa) mt ih} 


for, since the expressions under the integral sign have 
2m as a period, the result will be the same so long as the 
range of integration is 2m, independently of the posi- 
tion from which the range commences. 

For a Class I function with the origin suitably 
chosen, f(z) is an odd function and since cos u¢# is even, 
the product f(#) cos mt is an odd function; hence 
the integral of it from — a to + a is zero. ‘This also 
shows that in this case all the cosine terms vanish. 
Again, since sin mt is an odd function and the pro- 
duct f(z) sin nt is an even function, its integral from 
— 2 to + 2 is twice its integral from 0 to z. 

Hence the result that a Class I function of period 
2n may be represented by 


fd =bsint+hsinatet. . . (6 
where hee 2 fi)sinntdt . . . (9) 
1@) 
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Similarly for a Class II function with the origin 
chosen so that it is also an even function; J» vanishes 
since f(#) sin mt is an odd function, while a, is given 
by double the integral over half the range. 

Hence the result that a Class II function of period 
2n may be represented by 


fd) = 52 + ay cos t + ay cos 2t + PPO ES) 


where an = ‘i TCD ECOSHE aioe en) ee) 
4 6) 


Consider next a Class III function and break the 
integrals broken into two taken over the ranges — a 
to o and o to a respectively. If we consider corre- 
sponding elements in these two sections at a distance 
of x apart we see that f(z) has equal numerical values 
but opposite signs in the two cases, since f(t + 2) = 
— f(z). Further, since ~ is equal to an integral 
number of periods of the even harmonic terms, the 
product f(#) cos nt or f(z) sin nz will, when 1 is even, 
have equal numerical values but opposite signs at all 
corresponding elements and hence the integral from 
— za to o will cancel that from o tox. If, however, 
n is odd cos nt and sin mt will have an odd number 
of half wave-lengths and simply change sign in the 
range of , so that f(z) cos mt and /(2) sin nt will have 
the same values at all corresponding elements and 
we may therefore take the integrals from 0 to a and 
double them. 

Hence the result that a Class III function of period 
2n may be represented by 
f(t) = a, cost + ag cos 3t+ a,cosst+. . .) 

+ b, sin t+ d, sin 3¢ 4+ 5, sin st +. Jee 


where 
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an = A f(D cos nt dt 
oe ie) 
(11) 
and bn = A\ f@ sin nt dt 
ay [e) 


Lastly, let f(t) be a Class IV function with the 
origin chosen so that it is an odd function. Then 
since it is also a Class III function we have 


FE) a0 Sins si Yt be 12) 


where 
2 Tw 
by = “| F(Z) sin nt dt. 
wy Oo 


. IQ ve « 
If f(z) is referred to a new origin at t = —, it becomes 
2 


an even function with respect to the new origin; so 
too does sin nt when a is odd, for if n= 2m+ 1 it 


becomes sin (2m + ne + r) which = cos {mm + 


(2m + 1)t’} = (— 1)” cos (2m + 1)t’, hence the given 
integral from o to @ is equal to twice the integral 


7 
fromsdito =. 
2 


Hence the result that an odd Class IV function of 
period 2 can be represented by 


Wh) = Ope Oe ei ate, 
where 
4 m/2 i 
Ups A Wieaeiny Whar sees) 
e oO 


In an exactly similar manner it can be proved that 
an even Class [V function of period 2m can be repre- 
sented by 

f@ = .& cos t+ a, cos 34-4 a, cos 5¢ ++. « (14) 


where 
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m/2 
an = 4 TCE COS Ld am ere GES) 
Mo 


The student should guard against the mistake of 
working out these latter integrals when n 1s even and- 
then because the result is not zero putting in even 
harmonic terms in the expansion of f(z). The only 
reason why the range of integration is halved is because 
the values of a, and by are zero when n 1s even and can 
be obtained by doubling the integral over half the 
range when 2 is odd. 

In functions of Classes I, II or IV, only one kind 
of integral has to be evaluated, viz. [f cos nt dt or 
[fa sin nt dt; but in the Class III functions, and in 
the case of the general function, both types are required. 
In these cases it is most convenient to multiply the 
expression for by by 7 and add it to the expression 
Lor ay, so that 


Tv 
an + thn = it f(i[cos nt + 7 sin nt]dt 


Die in S 
= ty fe vdt; 


and then after evaluating this integral to separate it 
into its real and imaginary parts which are equated 
to dn and by respectively. 

If the period of f(t) is not 2x but an aliquot part of 


2m ; 
it, say — the formule we have given for the coeffici- 


ents still hold, since 2m is still a period of the function. 
In this case all the harmonics vanish except those 
which are multiples of the mth,* so that 


* Actually, of course, what appears here as the mth harmonic is 
physically the first or fundamental harmonic. 
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g 
ij = 3 + a, cos mt + a, cos 2mt + 
+b, sin mt+ b,sinzmti+ . . (16) 
where 
2 


I 
an =-— \ f(2) cos mnt dt; 
oe oO 


but since f(t) and cos mnt both have a period of 


O76 ; 
“— this may be written 
m 

27 
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O° 


and similarly for dy. 


Example.—Prove that if the period of f(t) be T, Fourier’s expansion 
becomes 


fl) = 2 + a, cos 7 + a, cos T+ ee ed 


ab 
a 2th . Amt . Ort 
+ b, sin 7 + 5b, sin bg sin = fy ee ee (LS) 
2 a 2mnt 
where an= Z| f(b) cos 2 ai Jes Me eee S) 
° 
2 - -. 2mnt 
and rn | 0) sin —— dt 1 je he ee) 
e} 


We may conveniently close this chapter with 


Some APPLICATIONS TO ELEecTRICAL ENGINEERING. 


Isolation of Particular Classes of Har- 


monics in any Periodic Function. If f(t) = “3 


“ig, coset... + bj sin ¢-+5..%.2., and if we 
displace it by a distance x, the displaced function, 


f(t + 2), is given by 
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feta) = 2+ a, cost +2) + Sle Dpto ere ee 


ae 
2 


— a, cos ¢ + a, cos 24 — d,cos3t + ...- 


— b, sint + & sin 2t — bg sin 3¢ + 
If we add this to the original function we get 


f(0) + fb + 2) = 2] 22 + ay 008 28 + a, C08 4 + ese 


+b sin 2t + Bysin get... 


that is, twice the sum of all the even harmonics. More 
generally, if we displace the function by the distances 


A ie da ae Abies pee in turn and add these 
{AY Tam P 

p —1 displaced functions to the original, we get, 
since the sum of all the sine and cosine terms which 


are not multiples of the pth vanish by (G), p. 9, 
the Sum = p | % + ap cos pit + dgpCOS2 pit ... 


+ bpsin pt-+ bypsin ape +. . «| 
rete 


that is, we get p times the sum of all the harmonics 
which are multiples of the pth. Again, if we subtract 


f(: + =) from f(z) it will be seen that all harmonics 


which are multiples of the pth will cancel out. As an 
example, let us find 

(1) The Terminal Voltage of a Three- 
Phase Generator. Let f(t) be the voltage in- 
duced in one phase of a three-phase star connected 
alternator. ‘Then the voltages induced in the other 
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phases are given by f( + 2) and f(: + =) the 
latter of which may be equally well represented by 
2m 
1— 2), 
f(r— 2 
The voltage between the second and third phases is 
f(: +. = — f(3 — 2) Since f(#) is necessarily a 


Class IIT function (p. 19) we have 


f(4) = a, cost + ag cos3t-+ a,cosst+ ... 
+ 5b, sint+ b,sin3t+ ..., 


and hence 


f(e+ = aie = 


which easily reduces to 
Mah 4, sint + a,sin 5¢ — a,sin 7t + a, sin IIt. 3 


+ /3{d, cos t— b, cos §t + b, cos 7t — b, cos11t+..}. 


If in this we write sin t = — cos (: a *) and 


cos ¢ = sin (: + *), Ctce awe (eet 
0/3 {4 cos(1+2)— a, COS 5(++2)-+2, cos 7(++2)— : 


+d, sin (147) —2, sin 5(74+2) +0, sin 7(1+7)- , 4; 
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If we compare this with f(z) itself we see that all 
the harmonics which are multiples of 3 have dis- 
appeared while the amplitudes of all the others are 
multiplied by ./3. We note also that the coefficients 
of the 5th, 11th, 17th, etc., harmonics have their 
signs changed, which may be interpreted as a change 
in phase of the harmonic in question amounting to 
half its period. 

(2) Voltage Wave Form of an Alternator of 
given Field Form. By the “field form” of an 
alternator we mean the voltage wave form induced 
in a single conductor. ‘This is not in general the same 
as the voltage wave form at the terminals of the 
machine, since the armature generally contains a large 
number of conductors connected in series which 
occupy a displaced position relatively to one another 
on the armature. We will suppose that the windings 
are uniformly spread over an angle of 2« radians, or 
om2e “electrical radians’ if the machine has) more 


than two field poles. If 
f G) == a, COs Fe a, COsia tae 
+ bd, sint+ 6, sin 3#+ ... 
represents the voltage induced in the central con- 
ductor, then 
f(t + 9) = a, cos (¢ + 6) + as cos 3(¢ + 6) 4+ . 
+ bd, sin (¢ + 6) + ds sin 3(¢ + 6) + . 
will represent the voltage induced in a conductor at 
an angular distance @ in advance of the central one, 


and the total induced voltage in all the conductors 
will be proportional to 


| + 6)d6, 


Chateis6to 
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104 a 
a, | cos (t + 8)dB + ay | cos 3(¢+ 6)d0+ ... 
—& 


+b sin(t+ 6)d@+ ... 


This is equal to 
a,{sin (¢ + «) — sin(t — «)} + alsin 3(t + «) — 
sin 3(¢ — «)} + b,{cos(¢— a) —cos(t+a)}4+ ...,/ 
or to 


: 2a: 
Za, sin &% cos ¢ + —#sin 3a cos3#+ .. . 
3 


+ 2b, sin asin t+ 29 sin 3a sin 3¢ + Anse 


If we divide this expression by 2 sin a, so as to make 
the amplitude of the first harmonic equal to its 
amplitude in the field form, we see that the amplitude 
of the mth harmonic in the voltage wave form is 


reduced (relatively) to ae of its value in the field 


a 


form. 
If, for instance, a = 2/3 or 60° all harmonics that 
are multiples of 3 will vanish and the field form 


TOY 4, 008 tod, COS 3r4-. 
+ b, sin t+ 5, sin 3¢+ ... 


will produce the voltage form 


Bicos t= 8008 SEY 0s 76 — oo ae 


+d, sin ¢— “sin st 2 sin 7t — as 


D 
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Example —lf «= ia or 72° so that the spread of the windings is 


144°, prove that the voltage wave form of a star connected three-phase 
alternator will be proportional to 
a, cos t+ *o8gaq cos 7t + ‘09144, cos I1t — *0484,3 cos 13t +. 
+ b, sin t+ ‘o89b, sin 7t-+ *ogtd,, sin 11¢ — ‘048d,, sin 13t +. 
where 
a, cos t+ ag Cos 3t- GetCOS |S | ets 
-+ 6, sin t-+ 6, sin 3¢-+ 3d, sin 5¢-- . . . 
represents the field form. 
[There is a difference of phase between these two expressions as 
written, but this does not affect the amplitudes of the various 
harmonics. } 


(3) The Root Mean Square (R.M.S.) Value 
of a Periodic Curve. We frequently, especially 
in electrical engineering, require to know the square 
root of the mean square of the value of a periodic 
function; f(#) being a periodic function of period 2z, 
this quantity is clearly by definition, 


els iP 

Es eke cde) 

Let us multiply both sides of the equation 
ft) = 2+ a, cos t + da, cos 2t + 


+ b, sin ¢+ b sin 2¢ + 
by f(#) dt and integrate from 0 to 2z. 


2 
Since \ f(d) cos nt dt = nay 
ie) 
27 
and \ fd sin nt dt = aby 
[e) 
we get 


ics he —(na,) + 4,(%4,) + de(a,) + i 
+ D,Gtb,) + DeGd) + 1S . 
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and hence 


Tai, Uma nO) tae gh 05? 
Pe en Nowra) | 2 So 1 1 2 2 
x 1 dt ri pee re 


If we denote the constant term %2 by c, and the 
2 


amplitudes of the Ist, 2nd, 3rd . . . harmonics by 
Eom nGst a ue) We esce- that: the. R.M.S; valuevot-a 
periodic function is given by 
A pug Came Tana aes 
2 

(4) The Form Factor. Another characteristic 
of a periodic function which is required in electrical 
engineering is the “‘ form factor” of a voltage wave. 
This is defined as the ratio of the R.M.S. value to 
the mean value of the rectified wave. 

Since the voltage wave of an ordinary alternator 
-can only involve odd harmonics, we limit ourselves 
to that case. Let the voltage be given by 

V=4,sint+ b,sin 3#+ . 
47 G4 COS f 4- dz COS 37 +> 
so that its R.M. S. value is 


aes +) 7) ae ae Tee 22 


We will suppose that . = 0 when ¢# = 0 and that 
b, is positive so that the positive half of the wave 
extends from t=otot=a. ‘This of course neces- 
sitates the relation a, + 4, +4,+ ...==0 among 
the cosine coefficients. 

The mean value of the positive half of the wave 


is ‘\ V dt, which is readily seen to be equal to 
7 Oo 


2 bs, Os ) 
ator et...) 
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So that the form factor is 


sade CEE Ue py ope C8 a Oe = 
2/2 b+ 6/3 -+4/5 +... 
The numerical factor is equal to I-111I approx. 
Whenever the equation of a periodic function is 
known, both its R.M.S. value and its Form Factor 
can be deduced directly from the integrals defining 
them: it is quite unnecessary then to use the expres- 
sions for these quantities given above in terms of the 
Fourier constants. 


CHAPTER “IIL 


ON THE REPRESENTATION OF ARTIFICIAL FUNCTIONS 
BY FOURIER’S SERIES 


WE shall in this chapter obtain the Fourier Series 
representing a variety of fairly simple artificial func- 
tions. ‘The Fourier series representing a function is 
called briefly the “Harmonic Analysis” of the 


function. 


Problem [.—Obtain the harmonic analysis of the 
eyenwiunction represented by /(7) = ) from +='0 
Po rand by f(t) 0 om tea tO f= 7.” 

Since the function is even, we have 


f(t) = 2 + ay cos t + ay cos 2¢ + is Seal) 


2 wT 
where an = \ T(t) cos nt dt. 
fe) 


If we substitute the given expressions for f(t) over 
their appropriate ranges we get 


2 o4 
ag N 4 cos nt dt, 
ch [@) 


: 2b . 
which gives an = — sin na. 
7 


* Tt will be assumed in all cases that the period is 27 unless other- 


wise stated. 
37 
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This readily gives all the a’s save a, which takes an 


5 . Oo ° . 
indeterminate form of —. If, however, we imagine 
fo) 


n to be very small but not quite zero, we see that the 


2hna _ zha ; 
— = ——, when # is zero. 


TUN Tt 


expression tends to 


2 © We Ye bd: 


Fic. 9. 


Hence substituting the values of the a’s in (1) we 
get 


{OQ = ze +10 > cos #4 — COS at 
sin 3a ] 
ean 3 COS'37 =P" atta: fo 2 Soak (2) 


In particular, ifa =” we get 
2 
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fa = ne -+ cos ¢ — 


cos 3f , cos St | 


F ale (3) 


A oe of this function is shown in Fig. 9(a) 
where the dotted line indicates the value of the first 
term. The way in which the successive terms 
approximate to the portion ABCD is shown enlarged 
in (6) where the curves marked 1, 2, 3 and 4 indicate 
the sum “ol the first 2; 3, 4 and 5 terms of +(3) 
respectively. 


+- 


Fic. 10. 


a. . 
We may note here that the constant term -2 is 
2 


always the mean height of the function ; for the mean 


2m 
height is = \ 40 dt, which is the value of me 


Problem II.—Analyse an even Class IV function 
which is equal to h from t =o tot=a and to zero 
7 


from t= a tot =>. 
Here 
f(t) = a 08 t + 43 c0s 34 +. 


where 


7/2 
dn = 4\’ f(t) cos nt dt = 4\" hb cos nt dt = ie sin na. 
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Hence 
_ 4h{sina sin 3a 
Hey rieen Cora aaleg cos 3¢ 


UN Cosas e ee ae or Pier ame 7) 
The graph of the function is shown in Fig. Io. 


The part from = to x is constructed from the part 


HirGemtite 


from =e to o by means of the relation f(t + z) 
= — f(#), which holds for all Class IV functions. 


Problem ha —Analyse an odd function represented 
by fi) = 7 from t= 0 to t=a and by = 
from t=a to t=a. This function is represented 
by Fig. 11. 

Unlike the two previous functions, this has no dis- 


continuities in magnitude but only in slope. 
Since the function is odd, we have 


f(t) = b, sint + bd, sin 2¢ + bg sin 3+ ... 


where bn = 2{" F(Z) sin nt dt. 
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If we integrate this expression by parts we get 


7 =e ole oe fs al I cos nt dt. 


cue dy 


This step would obviously be justifiable if f(z) was 
a single function. As it is not, we must break the range 
of integration up into two, one from o to a and the 
other from a@ tom, If we do this we see that the value 
of the first term at the upper limit of the first range 
is the same as its value at the lower limit of the second 
range and hence these cancel, leaving only the values 
at o and a as written. We might have said simply, 
that the first term only takes cognisance of the mag- 
nitude of f(t); and since f(#) has no abrupt changes in 
magnitude it behaves, as far as the first term is con- 
cerned, as an ordinary analytic function. Now /f(?) is 
Zero both when,+ =O and when t==«. Hence the 
first term entirely vanishes. ‘The second term, the 


: j b 
integral, breaks up into two; for f(t) = — from o tow 
a 


and = — - b from « to zx. 
vb oe 


Hence we have 


(o4 Tw \ 
b= aiAl cos nt dt — = {cos nt at| 
O 1 a 


mn : 


2) 20 (e no. sin no _ 2hsin no 
Y eabnrarrages oF 
nl ne na — a) a(t — a)n 


Hence 


= 


ing . Z : 
yi cast sin ¢ + sin 20 sin 2¢t we ine sin 3¢ 
vis Za 3 


eas 
) 
Sites Gad nell. Fase een Sin 
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Some special cases of this result are worth noticing. 
IU 
Vin og 5 We get 
£) = Sefsin s — Sn 3! 4 BSE _ TE +} 
gp 32 52 7 cS 4 
- (6) 


so that the function is a Class IV function, which is 
geometrically obvious from Fig. 11 whena = 2/2 

Again, let us put a = o so that the function takes 
the form shown in Fig. 12.5 


Fic. 12. 


The expression (§) requires evaluating as an indeter- 
minate form, since both numerators and denominators 


- sin 2% 
are zero.’ However, on replacing ——— by mn when 
a 


& = O, we get at once 


fi) = 2hf in, 4 Ze sin 3# , sin gt , i 
x | 2 3 4 J 
ee. ) 
Here the question arises, What is the value of this 
series when t= o? If we put t= 0 we see that the 
sum of the series is formally zero; but it really takes 
an indeterminate form, as is seen by imagining 1 
infinitesimally small when we get 


# FOURIER’S THEOREM 43 


Os 


which is of the form/o X 0} 

The value of the function at t = o is also indeter- 
minate: it is obvious from Fig. 12 that when ¢ = 0 
the function possesses all values from — h/ to + b. 
We shall prove subsequently (Chapter VI) that the 
series 1s also indeterminate when ¢ = o and that the 
range of indeterminateness of the series is 1-179 times 
that of the function. 


Problem IV.—Wet us analyse an even function 
which is represented by the same equations from 
o to z as the function in Problem III. 


ete (2) on a,cost+ a4,cos2t+ ... 


Tw 
where yee \ f(t) cos nt dt. 
de ie) 


If this be integrated by parts as in the last case, 
we get 


Cd a 2 | a = = ole f@ sin nt dt. 


Nn 


As before, the first term vanishes at both limits, 
and the second term gives 


io.4 Tw 
ly = — 2h (i sin nt dt — : sin nt dt 
eo | qa bi 


HANS Vo p i) 


2 cos ee ee 
an | CA % — % 

This takes different forms according as to whether 

m is odd or even owing to the cos mm term. If m is 


4b sin® na /2 
a(m — a) n® 


even it reduces to an = — 
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and if m is odd, to 
i= — I + cos a, 


Peers Lea ee Ue y 
he ae a) | n® | 


The value of a, which is twice the mean height of 
the function is obviously 5. ‘The shape of the function 


is represented by Fig. 7(c), p 


Problem V.—From the two previous problems we 
can at once write down the analysis of the Class Ill 


‘ a ae ht 
function which is equal to from ¢ = 0 toa and to 
a 


h(x — - 

Age. — 4) from t=atoz. For it must be represented 
Te CL 

by 


f(t) = a, cost + agcos3#+ .. 
+ db, sin t+ b, sin3t#+ ... 


and the values of these coefficients have all been 
found in the last two Problems, giving the result : 


ee ia sin 3a A 
7@) = olor Nigel. Sriy sin 3+ ... 


2a 2a 
— 1+ cosa 1 COGN. 


12 ae 
XT COST be an |. 
The shape of the function is represented by 
Fig. 7(d), p. 17. If in this series we put a = 7 


all the cosine harmonics vanish and the function 
becomes, as is geometrically obvious, a Class IV 
function. 
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Problem VI.—Analyse an odd Class IV function 


which is equal to HE from t = o toa, and to / from 
a 


wf 


~=a to gel 
2 
Here 
f(t) = db, sin t + b, sin 3t + . 


where 


1/2 
bi | fC) sin nt dt. 


Fic. 13. 


Integrating by parts we get 


7/2 7/2 
t) cos nt , 
Dy, as a| ie) ss n i -|- 4 ie (2) cos nt dt. 


Tale: 


This is justified, since f(t) has no discontinuities in 
magnitude. ‘The first term vanishes at both limits; 


y h : 
also f(t) = — from o to a and is zero afterwards : 
a py] 


hence 
b A sin no. 
b, = vl .sin na = me: 
mn an nan 
so that we have 


(Oe eee n¢4 30 3a 


32 sin 3¢ 


sin Pe 


+. 


Si Ghee iz. (8) 


—S 
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5 
7 : ip 
If, for example, we take a = —, we get sina = v3 : 
hic 
clink 2, = OR Gn SS ae etc., and so 


6,/2 : Sift Sf. 4 Sineyha Stuer a? 
fa) = w3 p{sint — a i os "+ ee 


or f(t) = bf{1-056 sin ¢ — -042 sin 5¢ + -022 sin 7¢# 
—-o0gsiniiz+ ... }. 


This function together with its first harmonic is 
shown in Fig. 13. 
Problem VII.—Analyse an even Class IV function 
2 
given by f(z) = (x — ) [POM =e OetOwr toe eS 


Fleres/(7)'—=«a, cost 174,.C0S 31.-- ae 


m/2 
where an = 4\ I(t) cos nt dt. 
Le oO 


Since f(z) has no discontinuities in magnitude we 
have, on integrating by parts, 


ee fe sin ud ‘ a ic Sint at un nO) 


TU nN tn O 


The first term vanishes at both limits ; f(z) vanishing 
Tt : : 5 ‘ 
at 5 and sin mt at zero. Since f(t) has no discon- 


tinuities in magnitude, we may integrate (9) again 


by parts, giving 
Pe een. ee cos | ne ha S45 


Tn —n man? 


7/2 
I(2) cos nt dt. 
fe) 


Again the first term vanishes at both limits. The 


integral can easily be evaluated, since f’(#) = — ge 
JU 
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rd 


LA 
from o to > Hence we get 


2h (7/2 2h 
a= orl cos nt dt = E 2 singe 
men J, 2 
The sine term is +1 and —1 alternately for 


i153, etc. and so 


iO 32h ee {cos # — we 4 £08 5¢ _ cos 7# a i (10) 


53 23 
This function is shown by the heavy line in Fig. 14. 


Differentiation of a Fourier Series. If we 
differentiate (10) we get 
sin 3t _ sin St 


“() — 32) Pt 1 5¢ | 
jab) 3 | oer my 
But f(z) is represented by the thin line in Fig. 14. 
The maximum numerical value of this function is 
pee: ’ 8t 7 
easily seen to be £ ; sinees{ (1) = — 3 from o to > 


If we call this maximum value H, the expression 
for f(t) becomes 
) 8H{.. sin 3f, sin St } 
pa) == (sin ¢ — 2 aa siege i hi 
which is in agreement with (6), p. 42, for the same 
function. 
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If we differentiate this again we get 
Debs eels ola © C05. 3h Con 5? es 
a= PH cos # : + : 
but f’(#) is the function represented by the dotted 
lines in Fig. 14, consisting of lines parallel to the 


axis and at distances H es above and below it. If 


this distance be called H’, we have 
” =| eee cOS Ge Tee 
14) = ees ; + | 


which is in agreement with (3), p. 39, if the axis be 
shifted parallel to itself so as to make the constant 
term zero. 

These examples illustrate the theorem that we may 
differentiate a Fourier Series term by term and 
provided that the differentiated series 1s convergent the 
result will represent the expansion of the differential 
coefficient of the given function. This result hardly, 
we think, requires further proof than is afforded by 
this example. Conversely, we may in a similar manner, 
integrate a Fourier series and the result will represent 
the integral of the given function. No qualification 
about convergency is necessary here, since the inte- 
grated series is necessarily more convergent than the 
original one. 

All the functions we have analysed so far have 
belonged to one of the four special classes. We now 
give one example of the general type. 


Problem VIII.—Analyse the function given by 
f(t) = bt/x from t= 0 to m and by f(#) =o from 
ft =a to 22. 

Since there is no symmetry about this function, we 
have 
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P. 
f(t) = <2 + a, cost + a, cos 2¢ + ort 
+ b,sint+sin2zt+ ... 
2m 
where pipe hn \ Fa at 
Tt (e) 


TT 
AON 
bi a 


Integrating by parts, we get 


int wT 
Oe aoe oa PAS iL ee on dt 


2 2 
TU in o «CU in 


in 

—— b : we e h enn __ a) 
pot SH} Sa ide 

emempertne that er == 1) cand el!" == cos tm 

(— 1)”, we get, on equating real and imaginary parts, 


bf 19 == i} and by = 2 (— Ieee 


hinge nn? | 


pcs 4 b 
The constant term “° is obviously -, so we have 
2 4 


2. ho (eS cos 3t , cos 5¢ i 
Se eee 


| ies sin 2t . sin Ae a) 

a\ I 2 ls 
Hitherto our artificial functions have been com- 
posed of straight lines or parabolas; we will now 
take some cases where they are composed of segments 
of sine curves. In the next two problems it will 
be convenient to take the period a sub-multiple 


of 22. 
E 
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Problem IX.—Analyse the rectified sine curve. 
The function consists of a repetition of the positive 
halves of the sine curve y = Asin#t. See Fig. 15. 

Here the period is x and the function is an even 
function, so that 


f(t) = 2 + ay cos 2t + ay cos 4t + as 


Tw 
2 ; 
where an = — \ A sin t cos nt dt. 
Ls 
(e) 
-7 (9) 7 277 
Fic. 15. 


But the curve / sinz is symmetrical about an axis 


of y through ¢ = = and so also is cos mt when 7 is 


5 Tt 
even; hence we may take the integral from o to ~, 


giving 


m/2 
an = ae \ sin t cos nt dt 
ate 


oO 
1 m/2 
= = ) {sin (nm + 1)t-— sin(n — 1)t\dt 
m/2 
_ 2b i (1 — 1)t  cos(m+ ve 
Tt n— I n+iI a 


Since ” is even, m — I and m + 1 are both odd, and 
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Pa 
therefore both terms vanish at the upper limit so that 
we have 


Ga 2 ume + z } ss Loe 4h 
mn A ARSE a(n? — 1) 
This assumes no indeterminate form when 1 — oO 
hence 
Mie ee en enol) 
: (11) 
The student of electrical engineering will recognise 
the rectified sine curve as the voltage wave form of a 


simple dynamo with a single coil in the armature and 
a two-piece commutator. The next problem will 
give the voltage wave form of an armature with three 
coils at 120° apart and a three sector commutator.” 


Problem X.—The positive halves of the sine wave 

= hsint are drawn at intervals of 27/3 and the 
ordinates added where they overlap. Analyse the 
resulting curve. 

This curve is shown in Fig. 16. 

In order to add the overlapping curves, we will 
draw the axis of y through a point of intersection as 
shown. The equation of the arc AB is then y = 


b sin (: +7) and of CD is y = Asin G- t). The 


* The “field form” is supposed to be sinusoidal in both these 
problems. 
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5 oo 
sum of these two is y = 2h sin = cos t= / cos #, and 


6 


and is the equation of 


2 4 aU 
this holds from — 6 to ; 


the arc EF. But the equation of the arc FG would 
be y = bcos t if the axis of y passed through its mid 
point; so that both arcs are similar and hence the 
required periodic function has a period of 2/3 or 27/6. 
The function is also an even function, so that, using 
the formule for the period an aliquot part of 2 
on p. 29, we have 


== as a, cos Of + a, COS 122 =. 
12 7/6 
where an = — \ f(2) cos nt dt 
he oO 


12h (7/6 
ee pees \ cos ¢ cos 6nt dt 
TU @) 


\ {cos (6n — 1)t + cos (6n + 1)thdt 


1 


_ 6b [2 (6m —1)¢ , sin (6n + va gk 
oa 6n — 1 6n+1 J, 


2 {an (nx — 2/6) . sin (mm + ) 
n 6n — I — 6n +1 ; 


Both sine terms are numerically equal to sin 2/6 
Or /2: By both are of opposite sign, so they give 
h 
(36? — 1) 
nately + and — starting with a, + ; a, is also obtained 
by putting = 0 in this expression. 
Hence we have 


an numerically, but with terms alter- 
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= al cos 6f cos 12 | cos 18¢ _ 
a ta ey, ate ta metiypmaite set 
= h@-955 + 0:055 cos 6f — 0-013 cos 122 
+ 0:006 cos 18 — .. . } 


which shows that the amplitude of the lowest harmonic 
is less than 6°% of the amplitude of the constant term. 


We will now take two mechanical problems leading 
to artificial functions composed of arcs of sine curves. 


Fic. 17. 


Problem XI.—A body performs simple harmonic 
motion in a fluid in which the resistance is proportional 
to the square of the velocity. Analyse the periodic 
function representing the resistance to the motion. 

If the velocity is represented by & sin z, the numerical 
value of the resistance will be given by Asin? z; but 
since the resistance always opposes the motion, it will 
change sign whenever the velocity reverses, so the 
resistance will be a Class IV function represented 
graphically by Fig. 17 where also the thin line shows 
the first harmonic. 

Here f(t) is represented by / sin? ¢ from t = 0 toa; 
hence 


fay = bjssin t+ 8, sin 38's 
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m/2 
where bn = 4b \ sin? ¢ sin nt dt 
7 0) 


1/2 
2 : 
= zy (1 — cos 22) sin nt dt. 
ve O. 
The first term of this integral gives 
7/2 

2h [s a f 

te ? 

aL—nd, 
which, since m is necessarily odd, is equal to 2h/xn. 
The second term gives 


h n/2- , 
zak | {sin (n + 2)t + sin (w — 2)t\dt 
ie) 


m/2 


2 lee (m+ 2)é , cos (m — 23) 


I n+2 n— 2 S 


Remembering that x is odd, this gives — Z( : 


z\n+ 2 
ie ) Sede Re noni. 20. ake 
u— 2 an* — 4) 
first term we get by = CEE and so 
oe ote PAaeSIN ha mein She erate ye oes \. 
m3 133-5 0 932577 957-9 


The loss of symmetry in this expression due to all 
the terms being negative save the first is in appearance 
only: really the factors in the first denominator are 
— I, I and 3. 


Problem XII.—A reservoir is divided into two by a 
vertical partition with a small hole in it near the 
bottom. The height of liquid in one part is kept 
constant and may be taken as zero; the height in the 
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other part varies according to kcos¢: analyse the 
function representing the velocity of flow through 
the hole. 

The velocity of flow is proportional to the square 
root of the numerical difference of heights and changes 
sign whenever this changes. So we have to analyse 
an even Class IV function of period 2% which is equal 


ia Iv . ° e 
to 4./cost from t= 0 tot = = This function is 


shown by the heavy line in Fig. 18,* the thin line 
representing the first harmonic. 


Fic. 18. 


Elere ifs Gy Coste Gs COS 3h} =~ .. 


4 t/2 
where an = Z| cos 4¢ cos nt dt. 
a 
fe) 


A Reduction Formula may be found for this integral 
which will make a, depend upon day-2, and thus finally 


upon 4,. Consider £ (cos 3/2¢ sin nt), where 2 is even. 


It is equal to 


— 3 cos 44 sin ¢ sin mt + 1 cos 3/2¢ cos nt 
2 


* By an oversight it is represented as an odd function in the 
diagram, 
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==7COS inf cos ¢ cos nt — : sin ¢ sin nt 
= so (on + 3) cos(m + 1)t + (2%— 3) cos(# — re}. 


Now integrate both sides of this from ¢ = 0 to t = 


Nia 


7/2 
The left-hand side becomes [cos 3/2¢ sin nt] which 
@) 


is zero since it vanishes at both limits. Hence we 
have 


7/2 2n—3 7/2 
\ cos ## cos (m+ 1)t td= -2=3{ cos #¢ cos (n—1)t dt. 
fe} 


fe} 
I 
Putting m= 2 we have @; = —-= 4, >“ putting 
7 
; I 
i — Ae we Obtained. — — 2 OE a AS Se and 
II 72 it 
so on. 


Elences)(¢)7=— ay cos t — : cos 3¢ + =f 2 cos 5¢ 


a) ee 
ate Te bees serent: } 
m/2 
where 2, = #{ cos 3/2¢ dt. 
we oO 


This integral can easily be evaluated by Simpson’s 
tule, taking about 5 intermediate ordinates. It will 
be found that Ge bli sp. 


* The student familiar with the Gamma Function will know that 
1/2 Vx T 
[. cos3/? 940 — =. ree (See Williamson’s Integral Calculus, 
Chapter VI and the Table of log I(p) at the end of the chapter.) We 


may point out that when 7 is large aq approximates to - es 
TH 
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If we integrate the expression for f(#) we shall 
obtain a result proportional to the total flow through 
the hole up to the time #, viz. 


sin 3f , sin §¢ 3 sin 7¢ | 
Fae) 530) 

A General Formula for the Fourier Constants 
in terms of the Discontinuities. We shall now 
obtain a very general formula for the Fourier con- 
stants from which the expansions in nearly all the 
previous examples can be written down straight 
away. Let f(t) be a periodic function of period 2x 
represented in different portions of its range by 
straight lines or by arcs of parabolas of any finite 
depree.* 

Let the function, or any of its differential coefficients, 
have discontinuities of magnitudes I., Ig, I, . 
V’., I’g Hast be pee eeat ethic  POINts 1051/0. ames 
using Bie notation of Chapter I, p. 14. We will 
suppose, to save extra words in the course of the 
proof, that there are no discontinuities at t= 0; 
although, since the origin is arbitrary, it is obvious 
that such merely require treating like any of the 
other discontinuities. 


We have f() = “2 + a cost+ ... +h sine + .. 


1-113h{sin t— 


where 
I 27 
ay Ag = fe dt. 
Le Oo 


Integrating by parts we get 


1p foent a es \ rem aren ae G3) 


in 


* A parabola of the pth degree is the curve given by the equation 
yoo yt cot? +... + opi? 
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Let now f(#) be represented by one analytic function 
from o to a and by another from a to 2x. The first 
term when taken between limits gives 


a fe — o)eina — f(o)eino + flan)einen — fla + o)eine| 
By supposition there is no discontinuity at o and so 
flan) = flo); -also ¢?*"==) 2°, so, we are left with 


pica erha 
| fla + 0) — fla — 0) 
1 
that is, with Bes 
wn 


If there are also other discontinuities at 8, y, etc., 
it is clear, by breaking up the range of integration 
into ranges from o to a, a to f, etc., that the first 


term of (12) will become + STyeina, where the 
un 


summation covers all the discontinuities of f(z) in 
magnitude. 

If we similarly integrate the second term of (12) 
by parts we get 


2 i] Fagen etl fet dt». (13) 


NIN in 242 


In an exactly similar manner the value of the first 
term of this expression is 
=e >) 
un" 
where I’, is the discontinuity of f(z) at t= a and 
the summation extends to all such discontinuities. 


Similarly, if we integrate the second term of (13) by 
parts we get 


ate [/ OS) bf | fer dt. 


mn? in min? 
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The value of the first term of this is similarly 


1 ; 
an? 


Proceeding in this manner, we get the final result 
that 


: 1 I 1 
n(an + ibn) = —2T,e"* — sre a te = eon 
n n n 


+ 21 acine gn CA) 


This series terminates; for if p is the highest 
degree of any of the parabolic arcs, f(+"(t) is every- 
where zero, so there are no discontinuities, I?+* or 
any higher ones. 


Separating real and imaginary parts we get 


I : I I ‘ 
Lan — —2 1, sin na — —2S 1g cos na + 21". sin na 
nN n* n 


I ttt 
+ il GCOS Mo. Ue ie ae CLA 
and 


I I : I 
04 = -21, cos na — =21, sin na — 1". cos na 
n ii n 


+ ae Se ey ee LO) 
n 


As an example, we will apply this last equation to 
the function in Problem I. Here we have a discon- 
tinuity + 4 at +a and one of —/ at +a; and 
there are no discontinuities in any of the differential 
coefficients, these all being zero throughout. Hence 


Eisai I : 2b. 
es ae, = =(— f) sin wa = 2 sin =n 
Man oP sin (— na) mat ) Ee a, 


the same as previously found. 
Again, let us take the function in Problem III. 
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Here only f(#) has any discontinuities; and has a 
b ha 
or 


discontinuity of magnitude — + ——— eee 
a a—a als — a) 


at +a. 


at — a and one of magnitude ests 
a(7 — o) 

Hence by (16), 

eee 


I ha 3 
SS ee = — $1n na 
n? a(x — a) 


n? a(n — a) 
ee) sin na 
a(% — a) "42 


sin (— na) + 


aby = = 


or by 


again in agreement with the value previously found. 
Formule (15) and (16) can be simplified in the case 
of functions belonging to the four special classes by 
consideration of the relation of the discontinuity at 
—a to that at + a: we leave this to the student. 

The formule (15) and (16) may be employed even 
when they do not terminate, provided that they are 
convergent. As an example, we will consider the 
rectified sine curve dealt with in Problem IX, save 
that we will now make the period 22 instead of x. 
The function is then represented by / sin ¢/2 from 
t=0O to ¢= 2. 

The function itself has no discontinuities in mag- 
nitude, but its differential coefficient has one of mag- 


: 2h ’ 
nitude = at? = 0; sothatI,= 4. Its second order 


differential coefficient has no discontinuities in mag- 


nitude, its third has one of magnitude — 2 at 
~=0; N16) that ie a =F similarly, ie — Es iRee ans 


= a etc., whence we get by (15) 
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tee (poe Nee a? 
tay, = elaine ) (4)... 


ee ie eae +...) 


n An? ~16n* 647" 
a Ae ee, eS 
n( ees Ane at 

4n” 


° o. 
Since the constant term is 0. we get 
A 


ape ae COstt mr COS 27.8 COS 37 } 
t = —_— —— oa rie 
) I.3 “Doe 5-7 


also in agreement with the previous result. 


EXAMPLES. 
1. Prove that the odd function of period 2m which is equal to b 


LLOMey— ORTORP—v70) 15 
sin 7a sin 5% sin 77 
a era et 


4° sin 
Ie tA: : ; 


2. Prove ae - even pert of period 2m which is equal to 7° 
from ¢= 0 to 7 is 


4 {fr  ©08 2 cos 3t | cos 4t Ne 
JO 4) 75 OI \s a 32 -- me sey 


3. (a) By putting t= 7 in Example 2 obtain the relation c= 


I i = 1 
ere a ar ae 


2 
(b) By putting t=— in (6), p. 42, obtain the result >= 


I I I 
er i ae eneae 


3 I 


(c) By putting to in (10), p. 47, show that a Te 38+ 
I I 
oe zr Poe 

(d) By putting t= ; in the result of Example 4 (below) show that 
a I lpeeet 
Soo aaa ar aa 


62 A PRACTICAL TREATISE ON % 


4. Prove that the odd Class IV function represented by f(t) = 

Breede ae 

o(B — +, ) from t= 0 to — is 

ener ot 2 

Sin 32 * ©sillS?) (asin .77 
ean ELEN MeN ier 
5. Four circles are drawn on the sides of a square as diameter. 

Prove that the polar equation of the quatrefoil defined by the four 

external semicircles is 


fo) = PEA sin + fk ie 


ne 4aV2{t , cos 40 cos 80 __ cos 120 ae \ 
fare eS a oat Ges) Meigs. 
where a is the side o&the square. 


6. Prove that the periodic functions of period 2m represented by 
sin mt and cos mt from — x to m are 


P bs sin t 2 sin 2¢ 3 sin 3¢ a 
sin mt == — sin me { - . - Rat 2 
a hem) ae Ra) ane ee 
and 
Be I m cos t Mm COs 2t ‘i 
cos mt = — sin mn{ = : ee fs 
1 Amar aE ZB pee J 
respectively. 


7. A long uniform chain hangs over a number of small smooth pegs 
at a distance 2m apart in a horizontal line. Each segment hangs in 
a catenary whose equation is y = A cosh #/A. 


Prove that the complete form of the chain is given by 
Les a cos x C0324 cos 3x \s 
ee ae rp! 1 4e to. | See 
Obtain this result, (2) by evaluating the integrals in the usual way 
with the help of the integral ee cos ax dx; and (b) by (15), p. 59, 


from the discontinuities of the function at x. 


8. Prove, by considering its discontinuities, that the odd Class IV 
function of period 2x which is represented from o to « by the line 
joining the points (0, 0) and («, 4); is represented from « to B 


by the line joining the points («, 4) and (8, &), and from 8 to : by a 
line parallel to the axis, is given by 


S(t) = 4, sin t+ by sin 3¢-+ bg sin 5t+ ... 


t= (C-feD swme (sin) 


where 
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g. Obtain the values of the coefficients in Ex. 8 when « = e B= 


I V/ Ser : 
b= and k= “2 so that the points lie on a sine wave. 


10. Prove that the even Class IV function of period 2m which is 
equal to cos ?¢ from ¢ = 0 to t= 7/2 is given by 


= eve eenk ews 
= C{cos 1 eg sae sare Rec OF 
pee Pes as ple ore P 
Baie Sat pe 7 ep 


COSTE clam oa 
where cat |e cos 2 + 1 dt,* 


sq ; ; Bier 
[Consider e (cos? + 'y sin n#) in a similar manner to that adopted 


in Problem XII.] 

Note that when ¢ is odd this series terminates; also that the given 
function is then represented by cos ?¢ throughout its whole range so 
that we have expressed cos ?# when is odd in a Fourier Series. 


11. A periodic function of period 2m has no discontinuities of any 
kind except at t= 7; that is, it is represented by one and the same 
analytic function from — 7 to-+ 7. Prove that its Fourier constants 
are given by 


I yi iM! 
man = (— PFE (yt t# St (3 
and 
I Ie iY 
nbn = (— 1)" 4 (—1" 1 (Ft 
where 1”) sy rn) =i (rn). 


Obtain the expansion when f(t) = cet!. 


st 
* The value of this integral in terms of Gamma Functions is 


es) 
a) 


(Williamson’s Integral Calculus, Chap. V1.) 
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12. Prove that the odd Class IV function of period 2m represented 


; To. 
by 4 sin ct from = 0 to t= > Is 


f= * at sin ¢ Bos sin $t = 


Toe g—¢2 * 25— 22 


By making ¢ tend to zero while he remains finite, obtain the expan- 


Tt 
sion when f(t) = ¢ from ¢= 0 to = 


13. Prove that the harmonic analysis of the even function which 
is equal to /cost from t= 0 to i= and to zero from t=" to 
fe ToS 

2b (1 cos 2t cos 4t _, cos Of i 
i) = — = ‘COs ¢ a6 648 
JO = =A eft amy 3.5 6 5.7 J 


all odd harmonics sate Ree itis the first. 


14. The positive halves of the curve / cos ¢ are repeated at intervals 
27 ' : 

of —. Prove, by the result of the previous example, that the harmonic 
n 


analysis of the sum of these curves is 


ee cos 2nt cos 4nt — 
T (2n—1)(n+1) (4n—1)4n4+ 1) 
cos 6nt i 
CFE Sy 
if m is odd and is 
2nh niz—1f cos nt cos 2nt 
antl a) Uin—1)@+1) (2mn— 1)(2n+ 1) 
cosgut 
2 (3u—1)3e+1) 7° } 


if m is even. 
Hence note that the voltage from a commutator with 27+ 1 


sectors is exactly the same as that from one with twice the number 
of sectors. 


“4: . ° Tw 

1$. I'wo sets of equidistant ordinates at distances of — apart are 
m 

drawn, the first set commencing at ¢= 0 and the second set at t= 


T : 
B50 that the members of the latter are midway between those of 
the former. ‘Through all the points of intersection of the first set 
with the curve f(t) sin¢ short lines are drawn, parallel to the 
axis of t, terminating at the adjacent ordinates of the ‘second eet, 
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¢ 
Prove that the harmonic analysis of the step-like sine function 
defined in this manner is 


ae sin (2m—1)t sin (2m-+ 1)t 
fl) = sin Ff sin s— RIN 8 ae 
sin (4m — I)t , sin (4m- 1)t _ | 
eee hie ee 
16. A bedy performs oscillations in a medium in which the resist- 


ance is proportional to the square of the velocity so that the equation 
of motion is 


ZS ey dia 
at2 + k di oe Cy == O, 
the sign changing every half oscillation. 

Prove, from the result of Problem XI, that, if the damping is so 
small that its effect on the period may be neglected, the motion is 
given by 

eee sin (nt + y) , 
ape 4A kunt 
3r 
» where A, and y are arbitrary constants. How do you explain the 
fact that when ¢ is large this result is approximately independent of 
the initial amplitude A,? 


17. Show that the result of Ex. 16 may be written 


ie ee eee 5 
ag CEE a ae 


where — 1, is the time at which the amplitude would have been infinite 
if the result was valid for large amplitudes. 


CHAPTER IV 


ON THE REPRESENTATION OF PERIODIC ANALYTIC 
FUNCTIONS BY FOURIER’S SERIES 


Ir is clear that any algebraic function of sines or 
cosines of @ and/or integral multiples of 6 must, since 
each of the sines and cosines is periodic in the interval 
2a, likewise be periodic in that interval, and therefore 
must be expressible in a Fourier’s series. ‘The Fourier 
constants will, of course, be given by the integrals in 
the usual manner; but they may be obtained in a 
much simpler way in these cases, for the Fourier 
expansion is nothing more than another representation 
of the same function which is already defined for all 
values of @ by the single given expression. It must 
be possible, therefore, to transform the given expres- 
sion into the Fourier series by mere algebraical manipu- 
lation. ‘Thus instead of obtaining the expansion by 
evaluating the integrals we may use the expansion to 
obtain the values of the integrals. 

To effect this algebraical transformation, we replace 
the sines and cosines by their exponential values ; 
then, for brevity, denote ¢® by a single letter, say w, 
and then arrange the expression in positive and 
negative powers of u, grouping uw” with u-". We 
finally replace u” + u-” by 2 cos n@ and u” — u-” by 
21sin v8. A few examples will make the method clear. 


Problem I.—Express, in the form of a Fourier’s 
series, the inverse square of the distance between the 


two points given by (7,0) and (a,0) in Polar Co-ordinates. 
66 
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The distance squared is a? — 2arcos@ + r%, so 
that we require to represent 
ees: = I 
a® — 2arcos@+ 7% a*¥—arutu-Nt+r 
u 
(42 — ur)(au — r) 
in ascending and descending powers of u. By the 
method of Partial Fractions, the expression can be 
written in the form 
P Q 
-|- 


a) 
A SUI. au = FT 


where P and Q are independent of uw. 
Adding these fractions and equating to the given 
expression we get 
Qa — Pr=o0 
and Pa — Qr = 1. 


Solving these for P and Q we get for the required 


expression 
I a im r ) 
a* — r*\a — ur au—r 


If 7 is less than a this may be written 


I I r I 
Spe POY ee : Pen bared fae 
a@—r\i—urla au 1—rfau 


the expansion of which is 


ali et les Os 
(ers Olen] 
pe grea) 8) aoe! 
Oi | 
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which is the Fourier series required. It is convergent, 
since r/a is less than unity. If7r/a > 1, we can obtain 
in a similar manner the expansion 


a : ile ae 2(“) cos 6 + (2) COS 20s la ers } . 
foes | r r 


Problem II.—TYo obtain the Fourier series repre- 
senting cos” @ and sin” 0. 

These expansions take a different form according 
to whether x is even or odd. 


We have 


= ae + 2nu2"-2 + 2n(2n = ey as 
22 | 2. 
[28 AR 2m 4 
Cat Cee 
+ ...H4 yan 
| 2 {} n(n — 1) 
aa CO ae ae f) mi 
22-1) n\n Bee gee a Wea cos 4 
n(n — 1)(n — 2) | 
Ge NG Gre ae J 
In a similar manner we have 
an+1 jam +1 { 3 n 
cos (ae ee in| es a cos @ -++ z cos 30 
n(n — 1) n(n — a)(n — 2) 
" (aaa a ae 


(n + 2)(n + 3)(m+ 4) 


fone th 
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and likewise 


sin?” § = eet ee, nen) 
22 nh 2 meri Gath CATES) 
cos 40 — nn = 1m = 2) 6 
4 ESN ee ea 2 a 
and 
| 
Site ace O a {sin Gee sins) 
2°" |n|n + 1 hoe 2 
n(n — 1) ; nn — 1)(n — 2) 
+ -—— —— sin 50 — = ve : 
(H+ 2)n+3)° > Get an 3+ 4) 
Sie /0 retain 
In an exactly similar manner it may be proved that 
cos? 4 cos 10 = ts + 7 cos 20 + mn = 1) cos 40 
HE : 
i Maa Ne Secon e eee } 


whether 2 be even or odd; and in a similar manner 
like 2xpressions can be obtained for cos” 6 sin 10, 
sin” @ cos ”@ and sin” @ sin v0, these expressions taking 
different forms in the latter two cases according as 
nm is cdd or even. 


Problem IIT.—Let us represent the displacement of 
an engine piston in terms of the stroke 2a, the length 
of the connecting rod /, and the crank angle 0. The 
expression giving the piston displacement, x, is 

x= acos@+ //?—a*sin?6—1 . . (1) 

Let the ratio a/] be denoted by @ 

Then we have, by the Binomial Theorem, 


x = acosé + at — 12 sin? @ — T 64 sint 
0 2 8 


—— 


LS bee Coa Sah a8 
— 0" sin’ @ = = sin’ Gs 244 — 1 
16° 128° 
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If we replace the powers of the sines by their expan- 
sions as found in Problem II, namely, 
2 $i” Uk== hy COs. 20, 
8 sint @ = 3 — 4 cos 28 + cos 48, 
32 sin°@ = 10 — 15 cos 20 + 6cos 49 — cos 60, etc., 
and arrange the result in the form of a Fourier’s series, 


ay easy Se ) 6 
2° Gree ae 
ais G au Us Sper Seas ) cos 20 
Ton yt 2 
- 3e 0 (& we) 6 
“G56 . cos 4 + eran cos 6 


This expression is chiefly required in practice to 
obtain the acceleration of the piston. If 6 = nt, so 
that 7 is the angular velocity of the crank, we have 


3 5 
dx = es CONGE—= (c -L eC -|- Eh or cos 20 


ae _ (9% ) | 
+(€ ate ee . -)cos 48 (Cee C05 ae 


The ave maximum of this expression is when 
6 =o and will be found to be w?a(1 + @) which can, 
of course, be obtained by differentiating (1), or the 
Binomial expansion, twice. 


Tue Direct SUMMATION OF FourIER’s SERIES. 


Closely allied with the preceding representation of 
an analytic periodic function as a Fourier series is the 
direct summing of a given Fourier series, which is 
simply the inverse process. How to tell whether a 
given Fourier series represents an analytic function 
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or an artificial one will be dealt with in Chapter VI. 
A single example will suffice to illustrate the method 
adopted in either case. 


Example I.—Find the sum of the two series 


c2 3 
eae tele. obs seen 


: co, Cob 
and eee pa ae tee ee 


where | ¢ | or 


If we denote the sums of these series by C and S 
respectively, we have, on multiplying the second one 


by 2 and adding it to the first, 
3 
C+ 18 = ce 1 © p20 4 as aa, 
2 3 


which is easily identified as — log(1 — ce’). 
Taking exponentials of both sides after changing 
their signs, we have 
e~(cos § — zsin 8) = 1 — ce", 


whence 
e-C§cos§ = 1 —ccosé and e-©sinS =csin# 
c sin @ 
and hence tan S = 
I — ¢cos0 
and 


e-2C(cos?S + sin?S) = e-2% = 1 — 2¢cos@ + c? 


and so 


C = — Flog (1 — 2¢ cos @ + ¢?) 
and Seastan os SN a 
1 — ¢ cos 0 
ligase ty these ‘simplity to 
C = log(4 cosec 4/2) (from 8 = 0 toz) 
and Sie. (from Ose O.t0 7), 
Sy 


the latter agreeing with the result in (7) p. 42. 
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EXAMPLES. 


1. The length of the focal radius vector of an ellipse of eccentricity 
¢ and semi-latus-rectum / is given by 


if 


it ee a, 
I — ¢cos 8 


Prove, by comparing this expression with that in Problem I, that 
the Fourier’s series for r is 
r= Isec B{1 + 2 tan B/2 cos 6+ 2 tan? B/2 cos 20 
+ 2 tan? 8/2 cos 30+ .. .} 
where sin 8 = ¢. Note that if a and b are the semi major and minor 
axes respectively, we have b= a cos # and, since / = 6?/a, the factor 
lsec B is simply equal to 6. 


2. Prove that the sums of the series 


c Ge 3 
1+ — cos 6 + -—— cos 20+ — cos30+ ... 
E 2 [3 
and 
ee Cas coe 
— sin §@ + —— sin20-+ — sin 30+ ... 
e 2 [3 
are e° £89 cos (c sin 6) 
and 6 6088 sin (¢ sin 6) 
respectively. 


3. Show that the sums of the series 


aml einige Do GEES. as 
Bey 


and sin §-+ u sin Ay b U ae sree 


|2 


2” cos" 8 cos (n + 1)0 
and 2” cos" 6 sin (n + 1)0 
respectively. 


are 


4. Show that the sums of the series 
1+ ¢cos0-+ c* cos 20+ c?cos30-+ ... 


and esin@-+ c?sin 26+ c3sin 30+ ... 
where | ¢| is less than unity are 
I — ccos0 c sin 8 
See aNd es weet eee 
I — 2¢cos@+ c? I — 2¢ cos0+ c? 


respectively. 
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5. Obtain the results : 


# 


cos 20 , cos 30 


COO es cl: ; + ... = log | 4 cosec 0/2 | 
cos 20 , cos 30 1 
oo ea ae . = — log | $ sec 6/2 | 
cos + SESE 4 SSE _. . = Flog | cot 6/2 | 
sin @ — ot ne ... = 4Flog| tan 0/2 + 2/4) | 


= 4 log|sec@+ tan 6| 
These results should be compared with the series 


eae: 


noe! sin 20 a sin 30 
2 3 
and 
; sin 20 _ sin 30 
sin § — 5 + 3 
which represent the function in Fig. 12, with the origins at the points 
marked o and x respectively ; and also with the series 


SEW Re cia eae rite 
3 5 
and 
an6 cos 30 pe a 


which represent the function in Fig. 9(a), where the origin is taken 
at the point A, and midway between A and D, respectively. 


CHAP TERY V. 


ON THE REPRESENTATION OF EMPIRICAL FUNCTIONS 
BY FOURIER’S SERIES 


Tue empirical function which we desire to represent 
by a Fourier’s series is necessarily given us either : 
(a) by some experimentally obtained graph; or (0) by 
some experimentally obtained set of points. In either 
case, the ideal results that were aimed at have only 
been reached to a certain degree of approximation. 
This limited accuracy of experimental results should 
always be remembered when subjecting them to mathe- 
matical operations. 

When a graph is given, two courses are open to us: 
we may either apply some form of mechanical con- 
trivance called a Harmonic Analyser to the graph to 
evaluate the Fourier Constants, or we may select a 
finite number of points on the graph and proceed as 
if these points were all that we knew about it. 

Several mechanical contrivances for the harmonic 
analysis of graphs have been invented which possess 
very various limitations with respect to the number of 
harmonics they are capable of dealing with; some- 
times, indeed, too severe to make the apparatus of 
much practical value. Space prohibits us from dealing 
with these machines here, since nothing but a com- 
plete explanatory description with diagrams and illus- 
trations would be of any value. A useful account of 
several such machines will be found in the Article on 
Calculating Machines in the 11th Edition of the 
Encyclopedia Britannica. References to the original 


74 


; FOURIER’S THEOREM 75 
description of several of these machines will be found 
in this Article. 

These machines are all both complicated and expen- 
sive, and, apart from the mechanical errors of the 
machine (and sometimes theoretical errors too), the 
accuracy of the results depends upon the accuracy 
with which one can move the pointer of the machine 
along the given graph. We do not consider these 
machines essential for practical use ; and a new method 
of harmonic analysis which is described below, we 
think, renders them all the more dispensable except 
perhaps where a very large number of similar graphs 
have to be analysed for the same harmonic. 

We will suppose, therefore, in the remainder of this 
chapter, that we are given a set of points on a periodic 
function, which we will suppose to be situated on 
equidistant ordinates, and that we require to obtain 
a Fourier’s series which will pass through the given 
points. ‘There are two distinct ways in which this 
problem may be solved. In the first method we find 
the Fourier’s series with the least number of the 
lowest harmonics which will give a function passing 
through the given points. This is the usual method 
of harmonic analysis and will be explained in detail. 

The new method referred to, which has not, as 
far as we know, been published before, is quite distinct 
and consists in joining the given points by straight 
lines (or by arcs of parabolas of the second or third 
degree) and then obtaining the Fourier’s series repre- 
senting this built-up curve from the discontinuities 
of the function defined by this series of arcs. This 
method is somewhat more laborious than the first 
one, but it has the advantage that we know exactly 
what the resulting series represents. ‘There is, how- 
ever, no utility in applying this method when all 
that is known is a set of ordinates of the given 
function; for we could not then tell which of two 
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different curves passing through the given points 
represented the function most accurately between 
the points. But the method is valuable when a graph 
of the function is given if we take ordinates at such 
a distance apart that it is obvious to the eye that the 
form of the curve between two ordinates is very 
approximately parabolic. In this way we may obtain 
a Fourier’s series which represents the given graph 
with much greater accuracy than does the simplest 
Fourier series passing through the given points. Since 
points can be measured on a graph much more accur- 
ately than the pointer of an instrument can be moved 
over the graph by hand; and since, also, with a reason- 
able number of ordinates, the various parabolic arcs 
will agree more accurately (through their entire length) 
with the given graph than a pointer moved over it 
by hand however carefully one does it; we consider 
this method of harmonic analysis as being more 
accurate than any Harmonic Analyser. 


The Simplest Fourier’s Series passing through 2p points 
per period. 


We will suppose that the period is 2% as usual and 
that we require to find the curve with the lowest 
possible harmonics which passes through the 2p points 


( ‘), (ee A,) ee (Ci hab) 


We are at liberty to assume any series with 2p 
unknown coefficients. If we assume 


f(t) = “2 + a cos t+ ag cos 2¢ + ... + ap cos pt 


+ bj sing + b, sin 2#+ ... + bpsin 97, 


we have 2p + 1 unknowns; more than can be deter- 


mined, but we note that when t = Ne pt = Na and 
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so the last sine term is seen to be zero at each ordinate : 

it therefore makes no difference to the curve passing 

through the given points, so we omit it. We will 

assume, then, the expression 

A ay COS mt ap— 1, COS (Pp — 1)t + “2 cos pt 
+ b,sin¢+ ...+ bp-rsin(g—1)t (1) 


oe Ap 
the reason for writing “° and “2 instead of a, and ag 
2 2 


is, as will be seen later, so that all the a’s can be obtained 
from the same expression. 

If in this expression we put t=, ae See or 
sueccssively, the leit-hand «ide becomes.A,,. A, 2 ..; 
Azp respectively by hypothesis. If we then add the 2p 
equations so obtained all terms cancel save the first 
by (F) and (G), p. 9, and so we have 


acy eee XAg. 
Zz 


PeNyemap Al anyiite ay =, == |, |. 270s SUCCESSIVELY 
‘in (1); but before adding the equations multiply them 
by cos ah cos aus pecoseL respectively, where 

i 

n+p and then add, we can see by equations (H) 
to (J), p. 10, that all the cosine terms will vanish 
save an, and that all the sine terms will vanish also ; 
and hence we have 


an& Cos? dis XAg cos 4, 
P p 


By the use of the formula cos? 0 = te we 
easily see that the value of the sum on the left-hand 


side is p, since there are 2p terms, except in the case 
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when ” = p or 2 = 0, when the sum is 29, but since 
we have allowed for this by writing 2 and = instead 
of a, and ap, we see that 


ee vere of he We aCe) 


in all cases. 
In an exactly similar manner, we find 


een Vicrieees Mem soe 2 es 
sah nd 
We have therefore solved the given problem. 


We will now prove that, as the number of ordinates 
tends to infinity, the constants a, and by, tend to the 
Fourier constants a, and Dy. 


Aq is the ordinate at t = 7: we may denote it by 
; : n 
A; instead. In this case the factor cos ?”” must be 
ie 
denoted by cos nz. "The factor 3 is very small, when 


p is large, and is = of the distance between successive 
ordinates; but this distance, in the limit, is dz, and 
hence ~ becomes a in the hmit when the summation 
becomes an integration. The summation extends 
from t=” to t= 2 which, in the limit, is from 


o to 2a. Hence 


27 
I na I 
= 2H Ag COS q”" becomes ‘| A; cos nt dt, 
P P Udo 
which is the value of the Fourier constant a,. Simi- 
larly, by» becomes by. 
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The expressions in (2) and (3) which consist of 2p 
terms can be reduced to one-quarter of the number of 
terms in the case of the four special classes of periodic 
functions. In what follows we shall denote the 
distance z/p between two consecutive ordinates by a 
and by Ag we shall always mean the ordinate at 


t= qa. 

For an odd function, let us divide the half-wave 
inte partseand set.the ordinates at0s,0, 20,05 7 = be 
0,57, Bo... a bp=7,0. Lhe complete set of ordinates 
over a whole period is thus 

O;.Da, B, eo» ie ‘s Bp-1, O, —Bp-1, eo) emo ra Bs, Oo. 


If these be inserted respectively for the A’s in the 
formula b, = . x A, sin qna, and if we note that the 


sine factors for B, and — B, are numerically equal 
but of opposite sign; and that the sine factors for B, 
and Bp: are also numerically equal, and of the same 
sign when z is odd and opposite when n is even, we see 
that 


by mee + (— 1)"-! Bp-r) sin 1% 
ks (3, + (— 11 Bp-2) aid 2na 


{ (8, + ( 1)"-# Bp-3) sin 3ma + .. | 
(4) 


the series continuing till all the B’s have heen taken 
once only. If there are an odd number of B’s, 1.¢., 
if p is even, the last term in this expression will be 


Bopj/2 sin”. The assumed series for the function is, 
2 


of course, 
f(t) = b, sin¢ + by sin 2¢-+ ... + bp-18in(p— 1) 2. 
ayer ee. = (5) 
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For an even function, let the ordinates be C,, C,, 


Cre: OR Fa eh A kn ary GO UN Wine 
ordinates over a whole period may be denoted by 
po ASH - Chega. Connie or Cy, Chee inserting 


these for the A’s in a, = aoe cos nq%, and grouping 


together, as before, the four terms which have the 
cosine factor of the same numerical magnitude, we get 


~ = 3 {Es ae = 1)"Cp bis ( ai 1)"Cp-1) BOoiTe. 


4. ( ie 1)"Cp-2) sie Bs | . .6 


the series continuing till all the C’s have been taken 
once only. The series in this case is 


f(t) = =2 + a, cos t + a, cos 2¢ + TS 


+ ap- AG ee 5 C08 Pt 5 aa) 


For an odd harmonic seal let us take the origin 
at one of the points where the curve cuts the axis, 
which it necessarily does at intervals of x. Let the 
ordinates over a half period be denoted by 0, D,, D,, 

. Dp-1, 0, so that the complete set of ordinates 
over a period are given by 


O, IDE IBY; : Dire Fy O> =a Dy, = Des a Dp- I- 


Senter these ee the A’s,in ae and (3) we 
easily get, on remembering that n is odd, 


an = (Ds = Dp-r) cos ne + (D. = Dp-2) COS 22% 


3 (D. — Dp-s) cos 3ma-+ .. } Se eS) 
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and 


be =={(D, x Dp-r) sin no + (D, + Dp-2) sin 2a 


i (D, cn Dp-s) sin nat. . | Pe) 


the series being continued till all the D’s have been 

taken once only. It must be noted that these expres- 

sions are not to be worked out for even values of n. 
The series representing the function is 


f(2) = a, cost + agcos3t+a,cossi+ ... 
a Dpsit-rop sb, sin gf =p. oP 3 (10) 


If p is even the cosines stop with the term 
ap—1 cos (p — 1)¢ andif pis odd with t cos pt, and the 


sine terms stop with bg-1 sin(p—1)t or bp-2sin(p—2)t 
respectively. It can easily be seen that the total 
number of terms in both cases is p. It may be noted 
that since f(o) = o the a’s satisfy the condition that 


Ayeajadg|-0on oie sO 


For a Class IV function, we need only measure 
ordinates over a quarter of a period. Let us divide 
the quarter period into p/z parts, p being assumed 


. “4 
even, by ordinates at — = a apart, as before. 
For an odd Class IV function, let the ordinates be 


7 
eet te, 2.) kupys At. O,g 20, . ~ Then from 


(9) we have 
lo aE sin ma + 2E,sin2ma+... 


+ Epa sin) ae et) 
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the factor 2 inside the brackets being associated with — 
all the terms except the last. ‘The series 1s 


f(t).= b, sint + bs sin 3#-++ ... + be-1sin(g — 1)t 
(12) 
For an even Class IV function, let p again be even 


and let the ordinates over a quarter period be F,, F,, 
wid pie2as 0. 9 Lhen-trom(6). since now, 


Cpe Ou Perey 5 Ci etc 
we have 
n= AF + 2F, cos na + 2F, cos2ma+ ... } 


(13) 
the series being 
f(t) = a, cost + ag cos 3t#-+ ... + ap-1cos(p— I)é 
(14) 
The student must remember that formule (11) and 
(13) are only to be worked out for odd values of n. 


There is still another great simplification which 
may be made in the practical calculation of these 
coefficients. Let us change m into p — m in any of 
the expressions for a, or by. Since gna i.e. gnx/p 
then becomes 94(p — n)x/p or ga — gnu we see that 
a cosine term will be unaltered if ¢ is even, and merely 
have its sign changed if g is odd, and vice versa with 
a sine term. So that the series for apn or bp-—n is 
always the same as that for an or by with the signs 
of alternate terms changed. Hence we add up the 
odd and even terms separately ; and the sum of these 
two sums gives an or by and the difference gives 
ap-n Or bp-n. In order to avail of this simplification 
in the case of odd harmonic functions, it 1s obvious 
that p, the number of ordinates per half-period, must 
be even. 
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Schedules for Harmonic Analysis. 


We are now in a position to construct Schedules 
for Harmonic Analysis. Such a Schedule for a Class I 
function when using ro ordinates per half-period is 
shown on page 84. ‘The given ordinates are denoted 
by the y’s; both yo and yro being necessarily zero. 
The sums and differences of these when arranged as 
shown are denoted by the U’s and V’s respectively. 
The first column contains the numbers from o to 5, 


; 5 oe ti) 
and the second contains the value of the Se sin 5 


when p = 10, for the values of m from o to 5. Nine 
further columns for the nine harmonics are taken, 
arranged and headed as shown. We then multiply 
the U’s in turn by the successive sine factors and 
place the products, which are denoted by the u’s, in 
the columns for the 1st and gth harmonic; placing 
them alternately in the third and fourth columns as 
shown. If the sums of these columns by denoted by 
Pac we have P, 4.0, band P, — ©, s=:b,. 
In filling in the w’s for the 3rd and 7th harmonics 
we place u, in column § in line 3; wu, should go in 
column 6, line 6, but the sine factor of this line is the 
same as for line 4, so we place it in this line; us should 
go in column 5, line g, but line 1 gives the correct 
sine factor, so it is placed there; ,, which should be 
in line 12, has the same sine factor as line 2, with its 
sign changed, so it is placed as shown. ‘The directions 
for placing these u’s may be stated as follows: starting 
with , in line 3 go up and down the lines from 0 to 5 
in steps of 3 at a time placing the w’s alternately in 
columns 5 and 6, and changing the sign of the term 
each time the sine factor passes through zero. For 
the 5th harmonic, 1, “,, and u, occupy line 5, while 
a and “4 occupy line 0; but since the sine factor of 


% 
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SCHEDULE III.—FOR ANALYSIS OF ODD HARMONIC OR CLASS III 


FUNCTIONS. (14 ORDINATES PER HALF-PERIOD.) 
Ordinates © Y,; Ye Ys Va Vs Vo Vr 
rey BA Biby BA 34) 376 
Sumas @ Wy, Wh Wh Wr We WG We 
Difference o Vy V2. Vs Vag Vs Ve V7 
F Harmonic. 
@| Sine | Cosine 
% factor. | factor. 
ist (& 13th). 3rd (& 11th). 5th (& gth). 7th. 
fo) fo) 1429 +0, 
—v 
I] :0318 | -1392 |u, v,|—4Us —v,|—Us —0; : 
2| 0620 | -1287 Ug | Ue U,| Uy +v,|+U6, 
3] ‘0891 | -1117 | 2, U3] Uy V1 |—4s —v, 
Ane DEL esOSIL U,| U4 —U,|—U, U,| Us 
Gil anekeyy, 0620 |\u; Us| Us Us| wy Vy 
6] -1392 | -0318 V6 | Ue Os|| Op —U4|—M%y 
7| +1429 fo) Uy 0 |—u, ©|| We O| U,;—Us 
+Us—Uy 
Sum of w’s P, Q 12) Q 1B 122 
Sum of v's Bul ces: | aR eS eae Be GUL 
Pi+Qi=b, Ps + Qs = bs P, + Qs; = b; P, = by 
P,—Q:= bj; P; —Q3 = by P; — Q; = by 
R,+S,;=a,/ R,; +S; =a, R, + Ss = a5 R, =a, 
R,—S,=443| R; —S3; = ayy R,; — Ss = ay 
The w’s are the U’s multiplied by the sine factor, and the v’s are the V’s 


multiplied by the cosine factor, of the line on which they are situated. 
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SCHEDULE III.—FOR ANALYSIS OF ODD HARMONIC FUNCTIONS. 
(14 ORDINATES PER HALF-PERIOD.) 


Ordinates oO 104 I92 264 320 360 384 392 
O 104 192 264 320 360 384 


SUM eee 0 208 384 528 640 720 768 392 
Difference o o fo) fo) fo) fo) Oo 392 
Harmonic. 
>| Sine | Cosine 
‘| factor. | factor. 
1st (& 13th).| 3rd (& 11th). 5th (& gth). 7th. 
) fo) ‘1429 
(| -0318 | «13902 lose) —22'9 —16:8 
»| 0620 | -1287 23'8 39°7 47:6 
BOSOl eran 147-0 18°5 —64'1 
B |) -1IL7 | 0891 715 — 85:8 42°9 
Mlci2o7 | -0020)) 92-8 68:0 26°8 
| -1392 | 0318 106'9 53°4 — 89:2 
ml 2b Aa2O fo) 56:0 — 56:0 56:0 29°7—70'4 
102:9—5'6 

Sum of “‘ sums.”’ | 202-4 | 262°2 7:6 FB} I'9 13 
Sum of “ diffs.’’ 

b, = 404°6 Da 14-0 [Ds = BA b; =1'2 

Dis = by = <3 bs = °6 

a = fo) y= © An = © a= 0) 

aaa = fo) aS Ae = @ 


go A PRACTICAL TREATISE ON 


TABLE OF SINE OR COSINE FACTORS FOR HARMONIC 
ANALYSIS. 


p or Number of Ordinates per Half-period. 


4 6 8 10 12 14 16 18 20 22 24 


fe) fo) fe) fe) fo) fo) 

*353 | +167 | -096 |:0618 |:0431 |-0318 |-0244 |-0193 |-0156 |:0129 |-0109 
-500 | +289 | +177 |:1176 |:0833 |-0620 |-0478 |-0380 |-0309 |:0256 |-0216 
333 | -231 |:1618|-1179 |-0891 |-0095 |:0555 |:0454 |:0378 |-0319 


*250 |*1902 |*1443 |*III7 |:0884 |-0714 |-0588 |-0492 |:0417 
*2000 |*1610 |1287 |-1039 |-0851 |:0707 |:0595 |:0507 
+1667 |*1392 |*1155 |:0962 |-0809 |:0687 |-0589 
+1429 |1226|-1044 |-089QI1 |-0765 |-0661 


"1250 |-1094 |:0951 |:0827 |-0722 
*IIII |-0988 |-0872 |-0770 
*1000 |:0g00 |:0805 
*0909 |:0826 

0833 


this line is zero, they are omitted. ‘The even harmonics 
are calculated from the V’s, but the arrangement of 
the v’s, which are the V’s multiplied by the sine 
factor of the line on which they are situated, being 
exactly similar to that for the u’s, need not be further 
explained. 

In an exactly similar manner, a Schedule for the 
analysis of any odd function using any number of 
ordinates per half-period can easily be constructed. 
A table of the sine factors required for a Schedule for 
any even number of ordinates from 4 to 24 per half- 
period is given above. 

Page 85 shows Schedule I worked out for a par- 
ticular example. The equation of the resulting 
Fourier’s series there obtained is 


f(t) = 19°43 sin ¢ + 3:08 sin 2# — 1-43 sin 3¢ — 1-38 
sin 4¢ + +73 sin 6+ -37 sin 7t — +33 sin 8¢ — +49 sin Qf. 

The given points will be found to lie on two straight 
lines and the Fourier analysis of this function is 
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; in t : 1 
Ge ee {sin = = + sin + ee | 


or 
f(4) = 19°27 sin 242-98 sin 2#—1-32 sin 3f— 1-20 sin 4¢ 
+°53 sin 6¢+-24 sin 72 —-18 sin 8¢—-24 sin 9f+ ... 
with which the Harmonic Analysis is in very good 
agreement especially for the lower harmonics. 

Page 86 gives the Schedule required for a Class II 
function, using Io ordinates per half-period. It is 
constructed similar to Schedule I and is an embodiment 
of formula (6), just as Schedule I is an embodiment 
of formula (4). Page 87 shows the working out of a 
concrete example. It will be seen that the given 
points lie on two straight lines. ‘The Fourier analysis 
of this function is easily obtained from its discontin- 
uities. It has four discontinuities, of slope only, at 


o,a and + ma From this we easily find that 


100 2n1 
ln 5 — I — cos mz + 2 cos —— 
BEA 5 


from which we obtain 

F(t) = 17 + 626 cos ¢ — 9:16 cos 2t — 1°82 cos 32 
— +87 cos 4z + -81 cos §¢ — -39 cos 6f — +33 cos 7# 
— *57.COS-8t —-08 COs Of +. ... 

as the Fourier analysis, in contrast with 

Ta) = 17 +631 cos f — 9°47 cos 27 — 1°96 cos 3% 
— I-99 cos 4t + I-00 cos §¢ —-+53 cos 6¢ —-51 cos 7t 
—1-o1 cos 8¢ + +16 cos gt 

as the simplest function passing through the given 

points. 

Schedule III, on page 88, for the analysis of Class III 
functions, is also constructed in a very similar manner. 
This Schedule, as far as the sine harmonics go, is exactly 
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like Schedule I for the odd harmonics. A separate 


Schedule can, if desired, be constructed for the cosine 
harmonics, but these two Schedules may be put 
together, as it will be found that the terms in one 
Schedule occupy the spaces left vacant in the other. 
For this purpose, it is necessary to have a column of 
the numerical cosine factors adjacent to that of sine 
factors. Further, it is necessary that the v’s should 
be distinguished from the w’s in the same column so 
that both sets may be added up separately. It is best 
if all the v’s, and the cosine factors, are written in red 
ink; or, alternatively, the u’s kept to the tig t of each 
column and the v’s to the left." €are must be taken 
that the v’s are not multiplied by the sine factors 
by mistake, which is not likely if different coloured 
inks are used. 

Class IV functions present no difficulty whatever. 
If arranged as an odd function, they may be equall 
well worked on Schedule I or HI and only half fill 
the Schedule as the “ differences” vanish. If arranged 
as an even function they may be equally well worked 
on Schedules IT or III, and in this case the ‘‘ sums ” 
of the ordinates vanish. Page 8g illustrates an odd 
Class IV function worked on Schedule III. The | 


: : : oi Gan 
given ordinates lie on the curve f(#) = ——* u(n—2), 
7 


from which it will easily be found that the Fourier 
analysis of the function is 


6 ye : 
fC) = oboe _ a sin ¢ + oe = a i Sra + 


or 
f(4) = 4046 sin ¢ + 15-0 sin 34 + 3:24 sin 5¢ 
+ 1-18 sin 7¢ +°56 sin gt +-30 sin 11¢ +-18 sin 13¢ 
sae S415 
with which the result 
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f(D) = 404:6 sin ¢ + 14:9 sin 3¢ + 3-2 sin 52 

+ 1-2 sin 7t +--6 sin of + +3 sin 11¢ + -2 sin 137 
obtained by Harmonic Analysis isin excellent agreement. 
The agreement with the Fourier analysis is seen to be 
very much better than in the case of the examples 
worked on Schedules I and II. This is partly due to 
the fact that we have used 14 ordinates in place of 
10; but is chiefly because, with this function, the 
discontinuities are in its second differential coefficient 
and not its first. The same function analysed by the 
use of 10 ordinates per half-period, but with the 
maximum ordinate reduced to 250, gives 

f(t) = 258-0 sin ¢ + 9°5 sin 3¢ + 2:0 sin Sz 

+°7 sin 7# + +2 sin Of 

against the Fourier analysis of 

J(4) = 258-0 sin ¢ + 9:60 sin 3¢ + 2:07 sin St 

+75 sin 7¢ + °35sinot+ .. 

There remains now only the general periodic func- 
tion to be dealt with. This is easily analysed by 
resolving into its odd and even components. We 
measure ordinates over a range from —a to +2. 
If yz and y-g are any ordinates at ¢ and —¢ then 


Yq + V-g is the corresponding ordinate of the even 
2 


component and ya So that of the odd component. 


These two functions are then analysed on Schedules 


II and I respectively. 
We frequently, in Applied Mathematics, when we 


desire to calculate a curve to pass among a number of 
experimentally obtained points, have a much larger 
number of points than disposable constants in the 
equation of the curve selected. In this case we cannot, 
in general, make the curve pass through all the points : 
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instead, we determine the constants so that the sum 
of the squares of the differences between the calculated 
and observed values of y, is a minimum for all possible 
variations in the constants. If the curve selected be 
represented by a Fourier’s series containing all har- 
monics up to the th, when ~ is insufficient to make it 
pass through all the points, it may be proved that the 
values obtained for the coefficient of any harmonic in 
this way is the same as if we had assumed an equation 
with sufficient harmonics as in the first part of this 
chapter. The proof of this will present no difficulty 
to the student familiar with the Method of Least 
Squares. Hence in dealing with periodic functions 
the difficulty of determining the constants of an 
equation (when we have less constants than points) 
so as to represent the observations as closely as possible 
does not arise: we simply adopt the ordinary method 
of harmonic analysis for the lower harmonics and just 
neglect those harmonics which are higher than we 
want to consider. 


A Graphical Method of Harmonic Analysis. 


Harmonic analysis may be performed graphically 
by the aid of diagrams which save the necessity of 
constructing Schedules. The accuracy in this case, 
though rather less, is often sufficient for practical 
purposes. Fig. 19 shows the method applied to 12 
ordinates per period. A large sheet of millimetre 
paper is taken and 6 lines inclined at 30° to one another 
are drawn as shown. 

These lines must be drawn as accurately as possible, 
their inclinations being set off by the tangents or 
cotangents of the angles, using as large a base as the 
paper permits, and not by a protractor. Three 
concentric circles are then also drawn as shown. ‘The 
intersections of the radii with the smallest circle are 
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then numbered 0, 1, . . . 11 consecutively. The 
intersections with the intermediate circle are alter- 
nately numbered 0, I, up to 11 and left blank. Finally, 
every third radius is numbered on the largest circle 
from o to 11 with the result that each numbered 
radius in this case bears three numbers. 


FIG. 19. 


The harmonic analysis is then performed as follows : 
We measure the ordinate A,, given cither by figures 
or graph, along Ow and read its length in millimetres ; 
then we measure off A, along radius marked 1 on the 
small circle and read off its projections both on Ox 
and Oy; we likewise measure A, along the second 
radius to the same circle and read its projections, and 
so on up to Aj. The sum of all the projections 
along Ox divided by 6 gives a,, while the sum of all 
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the projections along Oy divided by 6 gives b,. If 
instead of adding the sums of the odd terms to those 
of the even ones, as we have just done to get a, and b,, 
we subtract them, we get a; and b;. 

If we measure A,, A,, A, ... along the radii 
marked 0, 1, 2, . . . on the intermediate circle, and 
add up the projections and divide by 6 as before, 
we get a, and b,, while if we subtract the sum of the 
even terms from those of the odd terms we get a, 
and b,. Likewise a, and b, are obtained by using the 
radii to the largest circle, while a,/2 is one-twelfth 
the sum of all the ordinates. 


The Error of Harmonic Analysts. 


It is important to have some idea of the deviations 
of the a’s and b’s found by harmonic analysis from the 
true Fourier Constants the a’s and the b’s. We shall 
now show that the a’s and b’s can be expressed in terms 
of the a’s and b’s by means of the following 


Theorem. If 2p ordinates per period be taken at 
o, 7, ™ . . . and if each ordinate be multiplied by 
Rael te of cos mt at that ordinate, the sum of the 
products (symbolically denoted ener cos ann |p) 

g=t 


is equal to 


plan + aap—n + Gopin t+ ayp—n + aaptnt soe } 

cae (15) 
and if each ordinate be multiplied by the value of sin nt 
at that ordinate, the corresponding sum is 


Pion — bop—n + bop tn — dap—n t+ byptn— ~~. 3 
see (16) 
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# 
For we have 


Ag cos gon|p = |%2 + @ cos qz/p +... 


+0; sit guip—- ... } cos qu |p, 


and if the trigonometrical products be replaced by the 
sum or difference of two terms in the usual way, we see 
that the coefficient of 4, in the given sum is always 
zero, and the coefficient of ay is zero unless r — n or 
r+ mis zero or a multiple of 2, when it is equal to 
%.2p; from which it easily follows that the required 
sum is equal to p(an + dap—n + dopin+...). 

The second result follows in the same manner, 
but since the product of two sines is equal to the differ- 
ence of two cosines, the coefficients of the #’s are 
alternately + and — giving (bn — dDop-n + bop4n 

Since the left-hand side of (15) gives pan, we see at 
once that, in particular, 


a 
Sr dab 4p Gb6poam e 
2 Z 
ay = A, + Gap-1 1 Aopti tT A4p-1 + 4ypti T -- 


250 ag + Azp-2+ Gapt2+ Gr 


(17) 


Ap—1 = Ap-1+ Gp4+it G3p-1+ 43p4+1-+ Asp-i1 +... 
ap = Ap + Ap + 43p + 43p + asp 4 
this last can be written 


<P = ap + asp rasp + agp + 
Similarly, 

by = by — dept doptx— byp-rt dypti— -. | 

by = by — dap-2+ Doptan— ... - (18) 
bear bp—1 — Dp41 + b3p-1 — 43p41 + <enemne | 


H 
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If the a’s and the b’s diminish regularly and pro- 
gressively, we see from these equations, that the 
accuracy of the a’s and b’s regularly diminish from 
a, and b, till ag-1 and bp-13; which latter two cannot 
be very accurate, since we have, roughly, 


ap-1 = Ap-1 +- ap+t and bp-1 = bp-1 — bp+1- 
The remarkable accuracy with which 2 represents dp 


should be noted. : 
It should also be noted that in the equations (17 
each a occurs once, and only once, and similarly with 
the 6’s in (18). Hence each harmonic in the function 
analysed affects one, and one only, of the coefficients 
of harmonic analysis. For instance, if we are using 
10 ordinates per half-period, and if there is a large 
thirteenth harmonic while all the other harmonics 
beyond the ninth are small; all the a’s and b’s will 
be nearly correct except a, and b,; for we have in 
this case a, = dy + a3 + dog + Ggg + Qypt ... 
and by, = b, — bys + bey — b33 + Dy — «eS 


or, approximately, a, = a, + 4,3 and b, = b, — As. 


Determination of the Amplitude of any Selected Harmonic. 


We sometimes desire to find the amplitude of some 
particular harmonic, say the mth, which we are interested 
in, without troubling at all about any other harmonic. 
If the harmonic 1s completely unknown, three ordinates 
per period of the harmonic is the least which will 
suffice to determine it, although if we only want either 
the sine or the cosine component two ordinates per 
period will suffice. With a larger number of ordinates 
per period, the harmonic can naturally be determined 
more exactly. We shall give formule for the cases 
of 3, 4, or 6 ordinates in the period of the harmonic 
respectively, which will at the same time indicate the 
error of determining a harmonic in this manner. 
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If in the formula 
2p-1 nn 
2 Ag cos ee = plan + dap-n + aopin + ...} (19) 


we put 2p = 3”, so that the ordinates are at intervals 
of 22/32, we get 
Fee AA cy 

2 2 2 


* 2 
3n 
== as ayn ao aon +- Agn a Asn + aan +4 ar og 3 


and if in the formula 

ty F UN 
a Ag sin t = P{On ae bop-n + bapin—.. "| (20) 
q=o 

we make the same substitution, we have 


MSLA, = Ap + Am Ag+ Ay — Ag+ af | 


= 3 bn — bam + ban — Bon + Doe — Serie i 


Efsnow the 27th, anth, Svth'...-. harmonics are 
negligible compared with the uth, these formule give 
us approximately ay and dn, using three ordinates per 
period of harmonic. 

If in the formule (19) and (20) we put 2p = 4n, 
that is, if we use four ordinates in the period of the 
harmonic, we get 


A,—AgtA,y—Agt -- - =2n{an+agnt+agnt+am+ ...} 
and 
A,—A3;+A,;—A z+ . .. =20{bn—bsntdgn—Dant+ .. .} 


which give a, and by provided the 3th and higher 
harmonics are negligible. 

If we put 2p = 6m, that is, use six ordinates per 
period of the mth harmonic, we get similarly 
Oa pee eerie, va 

pp Z 2 2 


= 3n{an + den + An + Ayn + Asn+ . - es 


Wail, of cone We 
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and W3/A. + Ay — Ay— As sear te eel > + 
== 3n{by Den Da Dyin = Pisn ahora +S 


which formule give us a, and b, provided the 5zth 
and higher harmonics are negligible. Similar formule 
for eight or any other number of ordinates per period 
can easily be written down. The simplification of 
these formule for the four special classes of periodic 
functions may be left to the student. 

As an example, suppose we wish to find the amplitude 
of the third harmonic in the odd Class IV function 


Gerned. by (Ge ree ontade = 
Ut 


use four ordinates per period of the third harmonic, 
that is twelve per period of the function, the ordinates 
over the first quarter of a period are 0, 3, 6, 9, so we 
easily find 

= O'==.010, == 0g pi bye . 20} | OF Oy == —— Prapprom 


Using six ordinates per period of the harmonic we 
find that the ordinates over half a period of the 
function are 0, 2, 4, 6, 8, 8, 6, 4, 2, 0, from which we 


find - 
Y3(— 8) = 9[0s — dig + bar — ue < 


or b, = — ‘ 3 IO 


Actually, 4; for this function is — eS = — ‘SII. 
7 


A Simple Method of Harmonic Analysis. 


There is also a modification of this method of har- 
monic analysis, which we do not consider very satis- 
factory, but which may be employed when we are 
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given a graph of the function to be analysed. In 
this modification we draw different sets of equi- 
distant ordinates, and use one set of ordinates for 
determining one coefficient and another set for 
determining another, choosing the ordinates in each 
case so that we have a definite (small) number of 
ordinates to the period of each harmonic. The work 
involved in placing and carefully measuring such a 
large number of ordinates renders, we consider, this 
method more laborious than the one described, while 
biev-tesults are: generally less. accurate. Ihe only 
saving is in the simplification of the sine and cosine 
factors—or their abolition in the case when only two 
ordinates per harmonic are chosen and placed at the 
positions of the maximum and minimum values of 
such harmonic. 

This method may conveniently be developed from 
the formule connecting the ordinates with the co- 
efficients. If, for instance, we have a general function 
and neglect all harmonics beyond the third and take 
12 ordinates per period it can easily be seen that we 
have the following equations for the coefficients : 


ee a A a 
a, = H(A, — Ap + Ay— Ag + Ag — Aro) 

bg = (A, — Ag+ A, — A, + Ay — af 

Ge ee A ee A) 

by HAxs — Ags + As — Aro's) 


by — b, = #(As — Ay) 

Here Ax-; means the ordinate half-way between 
A, and A, and similarly for the other terms in the same 
line. ‘To use these formule, we thus have to place 
and measure sixteen different ordinates, while the same 
number of equally spaced ordinates would, with the 
ordinary method of harmonic analysis, enable us to 
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determine the first seven harmonics. All that has 
been saved is the trouble of multiplying the ordinates 
by the sine factors—a matter of very little trouble 
with a slide rule. 

For the general periodic function, a sufficiently 
accurate result for rough work can sometimes be 
obtained by assuming that all the harmonics beyond 
the third are negligible. Eight ordinates are then 
sufficient and it will be found that the ordinary 
formule for harmonic analysis give: 


8.2 =A, tA, + A+ Ag+ Ar t+ Ast As + A, 
Aaa e wget 77 aN heey tee Nee) 


Aaa fey ig Ne 

Ae Oy le eee eee) 

8. = (A, + Ar + Ay+ As) — (Al + Ast As + Ay) 
Ab =A, — A, 2707(A, A, = Ae =o A) 

be Ae te Ay A A 

fa = (Ay AD a) 70h eee) 


A New Method of Harmonic Analysis. 


We will now explain a totally different method of 
harmonic analysis which, though it is rather more 
laborious than the ordinary method, possesses consider- 
able advantage from the point of view of accuracy 
and is quite easy to apply when it has become 
familiar. As mentioned above, the method consists 
in joining the tops of the given ordinates with straight 
lines or by parabolas of different degrees. This can 
be done in quite a variety of manners, but the three 
varieties of the method given probably include the 
most useful that can be obtained. 

Let us, as usual, take 2p equally spaced ordinates 


Yoo Vina. - ~ per period: atio, wip. . Guoand dor) brevity, 
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let /p be again denoted by a. Let us join the tops 
of the ordinates by straight lines; then the slope 


between y, and ¥, is 1 20 and that between y, 


and y, is Je . ya ; so that there is a discontinuity of slope 


at y, of Y2— = ame that is, of Pea aanere Ary, 
MLA 


is that second difference of the y’s which is vertically 
below y,. Substituting this value for I’, in the 


formule (15) and (16) on p. §9* we get 


an = — PsA, cos 
mn 
(21) 
. 
and bn = Ree Sine 


I Gis 
But —2A*y,cos— is the an of ordinary harmonic 


analysis for the function whose ordinates are A*y.,, 


Zz 
A Vigne s,s 
Hence we have 


2 
; Sites age 2) 
and b, = — oe bs| 
2 n? 


The quantities an and by can be calculated in the 
usual manner on one of the given Schedules so that, in 
labour, this method of harmonic analysis exceeds the 
ordinary method only by that involved in forming the 
second differences of the ordinates. The Schedules, 


* The student must not confuse the present use of « with its use 
in these formule. In (15) and (16) « was used to specify the position 
of one of the discontinuities; now it is being used to indicate the 
width between successive ordinates; but no confusion need occur 
through this double use of this very convenient letter. 
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of course, only give an and by» when + 9, but if we 
write p + in the expression for an we get 


: gp + nn 
Abin == =n dey, COS 2 
Pies ? Vq ? 


which is the same as if we had written p — n, hence 
we have ap+n = ap—nand similarly bpin = — bp—n, 
and we see also that @ep+n =dn and Dien OF 
hence we are able to write down the complete analysis 
of the function defined by the straight lines joining 
the tops of the ordinates. In actual practice, of course, 
the function to be analysed does not consist of these 
straight lines and so the harmonics beyond the pth 
are too inaccurate to trouble about. 

As a very simple example, let us take the function 
analysed in Schedule I, page 85. If we form the 
second differences of the ordinates we find them all 


27 27 
zero save at —, where we have — Io, and at — —, 
where we have -++ Io. 

Hence we get 
fe) . 20n : 200 
b-b, = — ted — 10 sin —— + I0 sin (- =). 
Un | 5 5 
200 . 2nx 
or bn = —, sin ~_, 
al ad 5 


agreeing exactly with the harmonic analysis given on 
page 91, which it obviously should do. 

In a second manner of applying this method we 
draw a parabolic arc through the tops of the three 
ordinates ¥y,, 1, Ye; another one through the tops of 
Yes Ya, Ys and so on; we thus get a function with dis- 
continuities at the even ordinates y,, Ye, yy. . -, and 
at these ordinates there will, in general, be a discon- 
tinuity both of slope and of curvature. We must 
find an expression for the magnitude of these dis- 
continuities. 
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If we take the origin at y, it will easily be found that 
the equation of the parabola through the tops of y.,, , 
and vy, is 

APE M4 + 2 Je, 9 2, 


2a” 


where a is, as usual, the distance between the ordinates. 
Similarly, the equation of the parabola joining yp, Ys, 
and y, referred to an origin at 4, is 
> ia PY, 9° 

y= yp ta Mp Ma Da TI ya 
20" 
Differentiating and nas, t= — a, we see that this 
parabola has a slope of 


ia Salant aise 2 


20 


at Vo. 
Similarly, it will be found that the first parabola 
has a slope of 
Ja — 41 + Yo 


20 


at the same ordinate; and so the discontinuity of 
slope at yz is 


ee Ave Oe 4 eau) 


20 


4 
or — ap where A*y, is the fourth difference of the 
a 


y’s which occurs vertically below y,. 
We must now find the discontinuity in ay For 
the second parabola we see that a sans = ay ye 


: : : A 
or in our difference notation, —<*. Similarly, gy 
: a 
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for the first parabola is abt Lig and hence the discon- 
a 


Oat rer i : . A*y, — Ay, ; 
tinuity in = at the ordinate Yo 1s ce ok that is, 


the difference between the second differences which 
lie each side of that one below y,. Replacing a by 
z/p and substituting in the formule (15) and (16) of 
page 59, we get 


an aa eA Me cos oa 

a FRA, — A%y,) sin aa io = ee 
and 
be Fed J Snit oF 

= Fay, — Ay,) cos a oe ey 


The summation extends to all the ordinates at 
which there are discontinuities. ‘The expressions actually 
written under the 2’s represent the discontinuities 
at V2 only. 

If a’n and b’» are the ordinary harmonic analysis 
coefficients of the function with ordinates A*y, at Ya, 
etc., and a”, and b”, are those of the function with 


ordinates Ay, — A*y, * at y2, etc., we may write these 
coefficients 


pee ?* (¢ 3 " 
an = an =. tp nm. : A (25) 


27242 
and 


by = a bf = Paen . A . (26) 


nN 
2n2n? men 


* With the obvious notation of using the suffix 1°5 to indicate a 
term standing mid-way between terms of suffixes 1 and 2 we have 


A*y, a A*y, = A®yo.5 + Ay,.5 while Aty, = A¥yo.5 — AS y 1.5 - 
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As an example of this, we will take the case of an 
even Class IV function whose ordinates are 2-000, 
7 7 and™ respectively. Owi 

6 3 nd respectively. Owing 
to the symmetry, we need only obtain the differences 
over a quarter of a period, but to do so we must write 
down two ordinates before and after the quarter 


period chosen. We thus get the Table 


1-732, 1:000 and o at o 


r Te ™T T Te wT 2 Sr 
Saar 6 2 6 aaa 3 6 
y 7000 1°732 2:000 1°732 17000 0 -—I1:000 —1°732 
Ay 732 +268 -~—+268 -—-732 —1:000 —1:000 —:732 
A*y —-464 —+536 —-464 --268 0 +268 
A3y E*072 “072. “196 ~-268 -268 
BAY. "44. “124 -072 
A?y3—A’yy 0 268 +464 +536 


Thus, in a quarter of a period, we have two discon- 


. esd . . fA . 
tinuities in Ay, viz., at oand~. Ina complete period, 


there are two discontinuities similar to the first, 
and four similar to the second. In A*y, — A?y, we 
only have one discontinuity in the quarter period, 


Wize; AL 5 Substituting in (23) and remembering 
that 7 is odd and p = 6, we get 


= [ 
on ann? | 


2 * +144 cOsi0 + 4 xX +072 COs =m 
a2 Zal4 x +464 sin =| . 
Working out this formula, we find 
}() = 1°9975 cos t — -0097 cos 52 + -O081 cos 72 
— +00033 cos II¢ + -00163 cos 13¢ 
+ +00007 cos 17# + 00064 cos I1gt+ ... 
which agrees remarkably closely with the simple curve 


f(t) = 2 cos t which the given points evidently lie on. 
As an improvement in this method we may also 
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draw parabolic arcs joining the ordinates 4, yz, v3 and 
Ys, Vay V5 Ctc., and then add the analysis of this case 
to that of the former case and divide by two. This 
will be allowed for if we sum the expressions in (23) 
and (24) for every ordinate, instead of only for the even 
ordinates, and divide the result by two. We leave 
it as an exercise for the student to apply this method 
to the above example and find that the agreement with 
2 cos ¢t is made considerably better. 

A third manner of applying this method of harmonic 
analysis is to draw a curve joining the tops of ordinates 
y, and y, such that: at y, the curve is parallel to the 
straight line joining the tops of y, and y,, and at y, 
it is parallel to the straight line joining the tops of , 
and yz. Since this is virtually drawing a curve through 
four points, the curve must be a parabola of the third 
degree. ‘The tops of all the other consecutive ordinates 
are joined in a like manner. This function then has 
no discontinuities either in magnitude or slope; but 
it has discontinuities in its second and third differential 
coefficients. 

We leave it to the student to prove that the equation 
of the curve joining y, and yz is 

t 2 
y = Yo + (Ys — OR Cg Va 21) 


3 
(yg > 39a =r 3¥a— Vises 


the origin being taken on y,. The equation of the arc 
joining , and y, can be written down from this by 
writing y, for 3, y, for y, and changing ¢ into — ¢; 
and so is 


t 2 
De Pat Oe) ts 92 WAV Sack 2s) 


23 
Tina) ee ea, 


referred to the same origin. 
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The third differential coefficient of the first curve is 
constant and = 3(y, — 3s + 3%. — y,)/a3, while for 
the second it is 3(ys — 32 + 39, — yo)/a®. Hence 
the discontinuity in this coefficient at Ye 18 3(Va — 43 + 
Oye — 491 + ¥.)/a2; that is, 34%y,/a% in our previous 
difference notation. 

The second differential coefficient of the first curve 
at t= 0 Is (— ¥,+ 4ys — 52 + 24,)/a?, while for 
the second curve at the same point it is (2s — 5¥_ + 
4), — 4'.)/#. Hence the discontinuity in this co- 
efficien: at yz 1s (— yy + 273 — 2, + y,)/a?; that is, 
it is (-- A?ys + A*y,)/a?. If we reduce these suffixes 
by unizy so as to get the discontinuities at the ordinate 
y, and then substitute in the formule (15) and (16) 
of page 59 we get 


an = — £2(A%y, — 4*y,) sin mi 

a 3 As : “7; em aac, 
and 
oe yg — A®y,) cos i 

ge BPE Aty, sin a (28) 


Or if a’, and b’» are the ordinary coefficients of 
harmonic analysis for the function with ordinates 
A4y, at 44, etc., and a”, and b", similarly for the 
function with ordinates A?y, — A®y, at 4, etc., we 
have 


3 4 
" 3 / 
et. eats 
and 


3 4 
Dy = 4 avn + aE Oe ae ante (30) 


110 A PRACTICAL TREATISE ON % 


It is noteworthy that we are concerned here with 
exactly the same differences as we were when dealing 
with ordinary parabolic arcs joining three ordinates. 

Applying this case to the same example as we took 
previously (page 107), we get 


ay, = 


e I, x +268 sin + 4 x “464 sin 


eo 
3n3 | 
+ 2 X +536 sin =} 


NIC 


6 

+" 4X *072 COs a) 
Working this out, which is not the trouble it appears 
to be, for after having worked out a, there is very little 
to do save divide the same figures by v® and u4, we get 
f(t) = 1-99814 cos # + -00002 cos 5¢ + -00138 cos 7t 

+ sOOOOOI cos IIT¢ + *00003I1 cos I4t 
+ -000000 cos I7#+ .. « es eee) 


a remarkably good approximation to 2 cos ¢. It may 
be objected that the ordinary method of harmonic 
analysis would have given an even better approximation 
as it would have given 2 cos ¢ exactly. Quite so, but 
that is accidental and due to the fact that the selected 
points lie exactly on the curve 2 cos ¢ which is one of 
the terms of the assumed harmonic expansion. Had 
we, on the other hand, used the very similar function 
8i(m — t)/x? and taken 12 ordinates per period and 
joined successive even ordinates by a parabolic arc 
passing through the intermediate odd ones, our 
present method would, in this case, have given exact 
agreement with the Fourier analysis of the function, 
which is 


f(t) = 2:0641 sin t + -07627 sin 3¢ + -o1651 sin 5t 
++ 00602 sin 72 + -00283 sin gt + . : 


+ Oe X +144 cos 0 + 4 X +124 cos 


e895 
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whereas the ordinary method of harmonic analysis 
with the same number of ordinates would have given 
us 


J (4) = 2:0577 sin t + -07407 sin 3¢ + -01630 sin 52, 


the higher harmonics being quite undeterminable. 
The accuracy of this is “‘ miles behind ”’ the accuracy 
with which (31) represents 2 cos ¢. 

In the last manner of applying the present method 
the labour is about twice that of ordinary harmonic 
analysis apart from the labour of finding the fourth 
differences, since we have two different discontinuities 
at each ordinate and each discontinuity has to be 
treated just as if it were an ordinate. The present 
method is not put forward for use where any rough 
analysis is good enough, but only when the highest 
accuracy is required, and especially in those cases 
where we require the amplitudes of one or more of the 
lower harmonics very accurately, together with as 
high an accuracy as is obtainable, on the higher 
harmonics. ‘The last variety of the method described 
should then always be employed. 

The student who ever has any practical harmonic 
analysis to do will soon find he has had ample exercise 
in the subject. We would recommend all students, 
however, to construct one or two Schedules for a 
different number of ordinates per half-wave-length 
than those we have given. 


Interpolation Formule for Periodic Functions. 


When the value of a periodic function is known for 
non-equidistant ordinates the harmonic analysis can 
be performed by first writing down an interpolation 
formula and then transforming it by the method 
described in Chapter IV. 

If we have points, (a,A), (b,B) .. . (k,K) (LL) 
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the equation of the parabola of the m( — 1)th degree 
passing through the points is evidently 
a NO eee) 

i) tiem ee = yy 
(ied) Mia 0) eno ate) 
(TaN yas 

If now we imagine all these points lying between 
o and 2m and imagine the whole set repeated an infinite 
number of times at intervals of 2a to both — o and 
+ oo, we find that in the limit (32) becomes 


a 


relia (2) eens Ste area) 
I) SE) PANTO n 
sin HW 2) ee SI ee) 
aid —a)...sin#(l — k) ma 
when x is odd and 
ap NPOMLA Oe 0) onan SIDE Aad) 
Coa at) sin $(4—b).. .. sin 4(a—)) eras 
Fig esha saa ae LO STEAE 20) 
TACOS BU ont): eee ec ee 
(34) 


when 7 is even. 

The proof of these formule, though rather tedious 
to write out, is quite easy if the following results be 
borne in mind : 


RET I oe m2 

sin —X — —X I eT NO: I SS <a enue un 6) 
2 2 4n? 167? 
I I I I I 

$ cot —x¥ = —-+ - -- 
2 Xx 27 x — 2% x -- 4a 
I I I I I 
$ COSEC —¥ = +. 

Z G Rete Bing KN 206 x + 47 


The expressions in (33) and (34) obviously satisfy 


the conditions of passing through the given points 
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rd 
and, since each of the factors in the numerators 
merely changes sign when x is increased by 2% and 
there is in both cases an even number of such factors, 
the period is evidently 27 in each case. 

We know that if in (32) each given ordinate be 
increased by the same amount, f(x) will be increased 
by this amount for all values of x; that is, the formula 
(32) still represents the same curve : we conclude there- 
fore that formule (33) and (34) likewise satisfy this 


condition. 


UX 
If in (33) and (34) we write ¢2 =./u and arrange 
the result in powers of uw it will be seen that no 
fractional powers occur and that further, if the result 
is expressed in sines and cosines of multiples of x, 


ie) 
th 
iz 


the highest harmonic when 2 is odd is the 


5 5 nN 5 
and when x is even is the —th. But these are just the 
2 


harmonics which we found it necessary to assume present 
for our harmonic analysis when taking equidistant 
ordinates: we infer, therefore, that for equidistant 
ordinates (33) and (34) give the same result as is 
attained by ordinary harmonic analysis. 

When the ordinates are equidistant the expressions 
(33) and (34) may be greatly simplified. Let the 
points be (a,Y,), (2a,Y,) . . . (27,Yn), where a = 20/n. 
Then it may Be shown, by a well-known factorisation 
result in trigonometry that (33) becomes 

nx 


K(*) = on za — Y, cosec #(” — a) + 


Y, cosec #(% — 2a) — ... — Yn cosec #4 — 2n)| (35) 
which holds only when is odd; while similarly (34) 


becomes 
ji 
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f(x) = = sin a Y, cot Ha — a) + 


Vg COC 4 (= 20) eae ee COU ae 22) | - (354) 


which holds only when n is even. ‘That these expres- 
sions satisfy the conditions of passing through the 
given points and of having a period of 2m can be very 
easily verified. These formule, then, are mathematic- 
ally the same as the results arrived at by ordinary 
harmonic analysis: they suffer from the results not 
being in so useful a form from a physical point of 
view, but they possess the advantage that they can 
be immediately written down without any calculation. 

Formule (35) may be usefully employed for inter- 
polation with non-periodic functions when these are 
given for a fairly large number of equidistant ordinates. 
They are much more easily written down and calcu- 
lated from than (32) and give practically the same 
except near or beyond the end-points. 

We will now simplify (34) in the cases of functions 
belonging to one of our four Special Classes. 

First let f(%) be an odd periodic function and let the 
n given points (4,A) . . . (J,L) be considered as lying 
between o and x. If we make (34) in addition to 
passing through these points also pass through the 
nm points (— a, — A), etc., it will be found to become 


one sin # (cos~—cosb) ... (cos x—cos 1) 


sin a (cos a—cosb) .. . (cos a—cos 1) sane 
sin % (cos x — cos a) << . (cos x) — cosh) 
sin] ‘(cos 1 — cos a) . . . (cos 1 — cos k) 


(36) 
Similarly, if f(x) is an even function and if we make 
(34) in addition to passing through the m given points 


also pass through the m points (— a,A), etc., it will 
become 
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J 
De (c0st 0s 0). 2 (cos # — cos 1) 
I) (cos a — cosh) . . . (cos a — cos J) 
Rc08 A 005 4) (COS 2 COS TR) 
(cos 1 — cos a) . . . (cos 1 — cos &) 
(37) 
Again, if f(~) is an odd harmonic function, and if 
we make it pass in addition through the m points 
(a — x, — A), etc., it will become 
Rep chi Plc sin (4—" 2) 
ee ee Sina = 
sin(# — @)... sin(# — k) 
Per ene 


if m is even and 


Arcata” 


as Rea esi) 
ee ein = ante 


i hanes (39) 
if 2 is odd. 

Finally, if the f(x) is a Class IV function and if the 
nm given points be considered as lying between 0 and z/2 
and if we make (34) pass through similar points in the 
other three quadrants we get 
fy) = sin x (cos 2% — cos 2d) . . . (cos 2% — cos 2/) 

sin a (cos 2a — cos 2b) . . . (cos 24 — cos 21) 
eed Art ca eeee AO) 
for an odd Class IV function and 
fle) = cos x (cos 2” — cos 2b) .. . (cos 2% — cos 2/) 
cos a (cos 24 — cos 20) . . . (cos 2a — cos 2l) 
Ste PCR TACO) 
for an even Class IV function. 

These results are most readily obtained from (36) and 
(37) respectively. All the formule (36) to (41) 
evidently satisfy all the conditions imposed upon 
them. 
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The formule (36) to (39) can be further simplified 


to half the number of factors in each term if the 
ordinates are chosen symmetrically so that 1 = a — a; 
k =a —b, etc. In this case, we will suppose the 
ordinates at a and a — a to be A, and A, respectively 
and similarly for all the other ordinates. 

We then find that for an odd function (36) becomes 


om { + A, sin x rf A, — Ag sin 2%) . 


2 sin a 2 sin 2a) 
(cosi27 — cos 2h) Sins (ces 24-1 cou2)) 3 
(cos 2a — cos 2b)... (CEE SEA) mate!) 


Likewise for an even function (37) becomes 


doa {& — A, cosx , A, + A, cos 2*| . 


2 cos 4 2 cos 24 
(COS 227 COS) 20) een (COSA COG Ze) 
(cos 24a — cos 2b) . . . (cos 2a — cos 2/) peers) 


While for an odd harmonic function both (38) and 
(39) give 
nee {= + A, sin x 4: A, — Ag cos = . 


2 sin a 2 COs @ 
(cos 2% — cos 2b) . . . (cos 24 — cos 21) 
(cos 2a — cos 2b) . . . (cos 2a — cos 2)) +... | (44). 


In these three formule we have supposed that all 
the ordinates signified by the letters from a to 1 occur 
between o and 2/2. The term written down is the 
one for the ordinate at a and a similar term must be 
added for each of the other ordinates from 6 to J. 

In any of the formulz (33) to (44) we can take two 
of the ordinates undefinitely near together and so 
obtain the equation of a curve which has a given slope 
at a given point. By taking three ordinates indefinitely 
near we can obtain a curve with both a given slope and 
a given curvature at a given point. ‘Two examples 
will suffice to illustrate this point. 
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Example 1. Find the odd Class IV function of 
period 2 passing through the points (6,$6) and 


G H), when 6 is infinitesimal. 


Formula (40) gives us at once 


vas sin % COS 2% — COST 
sin 6 cos 26 — cosz 
sin x cos 2¥ — cos 20 
sin 7/2 cosa — cos 20 


or 
fiay— Sein costa 4- Hi sinex 
which may be written 


(2). ots sins + E Fe Ss Bn. 


Example 2. Find an even periodic curve of period 
2x through the points (0, 0) ( ) and (z, H), when 
6 is infinitesimal. 

Here (37) gives us at once 


S(*) 


_ (cos * — cos oXcos « — cosa) Ré? 
(cos 6 — cos 0)(cos d — cosa) 2 


(cos « — cos 0)(cos  — cos 9) 
(cos % — cos 0)(cos z — cos 0) 


or Eo We sin? a+ At — cos x)?. 


Problems of this kind, however, are often best solved 
by writing down a Fourier series, with the least number 
of the correct kind of harmonics which are seen to be 
necessary, and then determining the coefficients to 
satisfy the given conditions. 
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EXAMPLES. 
1. In an odd Class IV function, if 0, E,, E,, E, are the ordinates at 
vt ™ respectively, prove that 
Oe 5 tespectively, p 


by = 3(E, + E, V3 + Es), 
bs = 3(2E, — Es), 
bs = ME, — EpV3-+ E)). 
2. Show that b,, bg and b, in Ex. 1 are respectively, 
by = by — by + yg — O93 + 295 — 
bg = 43 — bg + P45 — ba + O97 — 
Be ee 


3. Analyse the odd Class IV function defined by f(#) = 2" from 
™ 


LO Wo - by using 3, 4, 5, 7 ordinates per quarter-period ; and show 


that in the first case the coefficients obtained satisfy the relations in 
Ex. 2 and satisfy similar relations in the other cases. 


CHAPTERS Vi 
THE THEORY OF FOURIER’S SERIES 


Tue present chapter may be entirely neglected by 
those whose interest in the subject is purely practical. 
Nevertheless, we hope we have made the subject so 
far sufficiently interesting to induce many, who have 
no intention whatever of making mathematicians of 
themselves, to peruse this chapter out of an interested 
curiosity in the subject. 

For the possible benefit of a few readers we will 
first recall one or two theorems on the convergence of 
infinite series. First, any infinite series 5, — 53 + 53 — 
5, +... in which the terms are alternately positive 
and negative, and in which each term is numerically 
less than the one before it and such that the terms 
decrease ultimately to zero, is convergent. For, 
writing the series in the form (5, — 52) + (53 — 54)+... 
we see that the sum is positive; and writing it in the 


form 5, — (5, — 53) — (4 — 55) — + » « we see that it 
is less than s,. Hence the series is convergent. 
Secondly, the series 1+4+4+4+... 1s 


divergent. For if the reader constructs a series of 
adjacent rectangles of unit width on the same base and 
of heights 1, 4, 4, etc., and also draws the curve 


he * through the vertices so as to lie inside the 
Be 


rectangles he will see at once that the sum of all the 


, GX ae 
rectangles is greater than \ ire: which equals log « 
i 
which is infinite. 
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Thirdly, the series I + oy + 7 + ; Ee 8 


n 


convergent when »>1. For if we draw rectangles 
as in the last case, and a curve through the vertices 
lying outside them, we see that the sum of all the 


Pde 
terms after the first is less than \ ee which equals 
I 


I 
n—1 

We will next prove that if 5,, 5, 53 . .° be any 
series of regularly decreasing positive terms decreasing 
to zero ultimately, no matter how slowly, then the 
series 


5, cos B + 5, cos 28 + 53 cos3B+ . . (1) 
and $, SID Por sy SIMIZ Ra kk, eee 


are convergent whenever f can be expressed in the 
form px/q, where p and g are integers of which p must 
be odd. Whatever the value of f this can be done to 
any desired degree of accuracy. If, ¢.g., B = 42/7, the 
theorem will hold for 6 = 4,000,0012/7,000,002, so 
the two series can only be divergent for certain isolated 
values of B. 

Consider the first 2g terms of either of these series 
and divide these terms into two groups, the positive 
and the negative terms. Let the sum of the first be 
denoted by P, and of the second by —N,. If 
Sy cos rB, say, is any term in P, then s,44 cos (r + 4g), 
that is, — 5,44 cos rf is a term in N, which by hypo- 
thesis is numerically less than the corresponding term 
sy cos rB in P,. 

Similarly, all the terms in N, are less than the 
corresponding terms in P, and hence N, <P,. If 
now we consider the terms from 52741 to 54g and denote 
the sum of the positive and negative terms by P, and 
— N, respectively, we see, by a similar argument, that 
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Po N, and that IN, =< P, and so on. Hence the 
series of terms, 

P,—N,+P,—-N,+P,;—..., 
are such that each is numerically less than the preceding 
and of alternate sign. It follows in this manner that 
both the series (1) and (2) are convergent. 


Writing 6 + ¢ for B in (1), we see that 


5, cos B cos t + 55 cos 28 cos 2¢ + “1 (3) 

— 5, sin Bsin t — sy sin 2B sin 2t — he 
is convergent as a whole, for all values of ¢. But the 
first line is an even function of ¢ and the second one an 
odd function. Now the infinities of an even function 
cannot everywhere neutralise the infinities of an odd 
function, for if they neutralised for positive values of t 
they would re-inforce for negative values. It follows, 
then, that the series for both the odd and the even 
functions must be separately convergent; and also, 
that a simular result must hold for the two series, 


5, cos PB sin t + 5, cos 2B sin 2t+ . +} (4) 
and 5g Sin B cos ¢t + 5g sin 28 cos 2t + 4 4 
obtained by putting B = 6 + ¢in (2). 

One obvious exceptional case must now be noted, 
that is, when 6 = oora multiple of 2x in(1). All the 
terms are then positive and the series 5, + 5,-+53+. 
may be either convergent or divergent, according to 
the rapidity with which the s’s decrease. We conclude 
from this that either of the series in (3) may be infinite 
when $8 + ¢ is zero or a multiple of 27, but must be 
finite for all other values of 7.* 

Each term of the series (2) is identically zero when 
p=. eput. the series may then behave@in a very 


* If the reader draws, from an origin, radii of lengths 5,, 5... ., 
in the directions 8, 28 . . ., respectively, it will be nearly (if not 
quite) obvious to him that the sums of (1) and (2) are necessarily finite 
for all values of 8 save zero or a multiple of 27. 
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peculiar manner which must now be investigated. 
Let us suppose that £ is indefinitely small so that we 
can replace sin mB by nf, the series then becomes 


SB tre 25a 8 t= 38 e0/te ee 
or Als, + 25. + 3se+ ...]. 

Now the series in brackets may be convergent or 
divergent: if convergent there is no doubt about the 
sum of (2) being zero when f = 0, but if divergent 
the result takes the form 0 X « when B = Oso that 
in this case we cannot say anything about the value 
of the series 

5, sin B + sg sin 264+ . 


when 6 =o except that it is formally zero on the 
strength of the fact that each term is zero. 


s I . . 
Let us suppose that 5, = —, so that the series is 
n 


(8) =F nae f ne 2. GAN) 


We may write this 
sinf , sin 2B , sin 36 sin NB | 
| Dare oo Gy he 


The student must now recollect that the limit of 
the sum 


BUI(B) fet . +. - CNB)] 


when 8 = 0 is 
NB 
10) 


by the definition of a definite integral. 


Hence the limit of the sum of the above series up to 
the Nth term is 


x. 


\ sin x 


O x 
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If now we make 6 = 0 and N = «x, we do not get 
a unique value for the result, but a value which is a 
function of the product Nf. We conclude that the 
given series, when 8 = 0, may have any value that 
the integral 


can take. This integral clearly increases with X from 
X =o to X =a, when it assumes the value 1-8519; 
after then it decreases till X = 2%, owing to the 
additional elements being negative, when it takes the 
value 1-4181; it then goes on oscillating with a 
diminishing amplitude till X = o, when it has, as we 


IU 
saw on page 11, the value of St 1°5708. 


This is the behaviour when £ has an infinitely 
small positive value. Since the series in question is an 
odd function of f, it behaves in an inverted manner 
for negative values. 

Hence the series (5), when B =o, can assume 
any value between — 18519 and + 18519. We 


shall call the value 2 the principal value of the series 


when # is +0, and — x the principal value when 


B = — o, and zero the formal value of the series when 

B =o. The principal value corresponds to the value 

when Nf = «, and is thus the value of the series 

when f is an exceedingly small fimite quantity, or 
. 7t “A 

as we may express it, ¢(-++ 0) = 5 and ¢(— 0) = — a 

We can similarly investigate the behaviour, when 


fp = 0, Of the series 


4 sin B+ 5s,sin2B+... 
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when sy = 5 and p is any positive quantity. The 
n 
series becomes in this case 
sin B , sin 26 sin NB 1 
en cee an ee 


Fic. 20. 


which, when £ tends to zero, becomes 


Ne .: 
a ae (6) 
Oo 


xP 


The student familiar with definite integrals will 
recognise that this integral is finite and convergent 
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whenever p<2. So if 1<p<2 the result of (6) 
is zero on account of the factor f?—*. It can be 
proved that the result of (6) is zero even if p> 2, 
since the zero of the factor f?~—" is of higher order 
than the infinity of the integral owing to the factor 


xP in the denominator. If, on the other hand, 9 lies 
between o and I the result of (6) is infinite when B = 0, 
unless the range of integration Nf is infinitely small, 
owing to the B?—* factor which is then infinite. In 
this case the series can take any value from o to « 
when 8 = + 0, and any value from o to — « when 
8 = .— 0, since it is an odd function. 
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The appended diagram (Fig. 20) shows qualitatively 
the behaviour of the series 


Siiep usin 2p 
{P- T 2p 
for the values of p of 4, 1, 3/2, 2 and 3, in the neigh- 
bourhood of B = o. 


Fig. 21, similarly, illustrates qualitatively the 
behaviour of the series 


cos 2p 
2P 


cos B 
1? 


+ : 


for the same values of p in the neighbourhood of 
pie 0: 

In both these diagrams the value of the series is 
infinite when p <1 and also when p = 1 in Fig. 21. 
The length of the heavy vertical line in Fig. 20 shows 
the extreme range of undeterminateness of the series 
(5) when ~ = 0, Both: arcs, for the value p ==73/2- 
in both these figures, have an infinite slope but a 
finite radius of curvature at B = 0; in both diagrams 
the slope “at 8 = 0 isi finite ip = 2, ana initesron 
‘p= 2in Fig. 20: thercuryaturesat p= Ors nites 
both cases when p > 3 and for p = 3 in Fig. 20 if 
the arcs on the two sides of the origin are considered 
as separate curves.* 


The Relation between the Fourier Coefficients and the 
Properties of a Function. 


It is not necessary that a periodic function /f(/) 
should be everywhere finite for it to be expanded in 


* ‘The arcs on the two sides of the origin are necessarily represented 
by different analytic functions, and each arc has a finite radius of 
curvature at the origin. ‘This must not be confused with the point 
of inflection of an analytic curve for which the radius of curvature 
increases to infinity as the point of inflection is approached. 
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a Fourier Series; but it is necessary for all the Fourier 
constants to be finite, and for the higher ones to 
decrease indefinitely, if the series is to be convergent. 
The first condition will obviously be satisfied if 


Taalat ey eae ee) 


2m 
is finite; for the integral | f(#) cos nt dt consists of 
ie) 


elements none of which are larger than corresponding 
elements in the first integral; while some are taken 
positively and some negatively in the second case, and 
all positively in the first case. 

We will next prove that if f(#) is everywhere finite 
the Fourier constants an and bd, cannot be greater, 
when 71s large, than c/n, where ¢ is some finite constant. 
In other words, the Fourier constants, 4,, ad), a3 . . + 
must decrease, ultimately, at least as fast as the series 
he ba se ar 

The function f(z), in the interval from 0 to 2m, can 
be divided into a finite number of ranges in each of 
which range it is of one sign and either constantly 
increasing or constantly decreasing (or stationary). 


Let us consider the value of the integral \ (2) cos nt dt 
over any of these ranges when z is fairly large. Over 


a short range of ~ the integrand will be positive and 
n 


will contribute s,, say, to the integral. Over the next 
range of a/n the integrand will be negative and give 
— 5g, say, and so on up to 5, where s, corresponds to 


the last complete half-wave in the range considered. 
Now the terms 


Sms Spe Sg Sas ES; 


are alternately positive and negative and either steadily 
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increasing or steadily decreasing. In either case we 
can see by a similar argument to that used on page 119 
that the sum cannot exceed either 5, or s,, whichever 


is the greater. Also, it is clear, neither of these terms 


2 a OMI ; 
can exceed “. MM.“ = ~~, where M is the greatest 
7 n 


n 
value of f(#) in the interval considered. To this result 
there may have to be added the contribution from any 
fractional half-waves at the beginning and end of the 


: : 2M 
interval, but neither of these can exceed ——, so that 
n 


the value of \ 70 cos nt dz over the range in question 


cannot possibly exceed ——. Summing for other ranges, 
n 


we see that the whole integral from 0 to 2m cannot 
exceed c/n where c is some finite constant. The same 
result obviously holds for the sine integral. 

From this important result we can easily prove 
that if f(z) has no discontinuities in magnitude, and if 
f (2) is everywhere finite,* then the Fourier constants 
an and by, cannot exceed c/n? where c is some constant. 


For 
; rea Wee 
an + iby = - \ f@e ae; 
ie) 


IU 


integrating this by parts, we get 
aah itm = [EY 
aL in 


If f(t) has no discontinuities in magnitude the first 
term vanishes, since it is periodic. ‘The integral gives 
the Fourier constants for the function /(#) and these, 


2m 2 
Sanlg \ (dent dt. 


og tse 


* The first condition has been stated to prevent misconception : 
the second condition includes the first. 


FOURIER’S THEOREM 129 

- 
we have seen, are <c’/n where c’ is some constant, 
since f’(t) is everywhere finite. Hence ay-+ ibn, and 
therefore a, and bn, are < c/n? where ¢ is a constant. 

In a similar manner, we can prove that if there are 
no discontinuities in f(t) and f’(4), and if f(z) be 
everywhere finite, then the Fourier constants a» and 
bn must be < c/n? where ¢ is some constant, and so 
on for discontinuities in the higher differential 
coefficients. 

We will now investigate the order of magnitude of 
the higher Fourier constants when f(t) becomes 
infinite at one or more points. Ifit becomes o at more 
than one point in the range 2” we can resolve it into a 
number of functions each of which becomes infinite 
at only one point. Further, in considering any of 
these functions, we may take the origin at its point of 
infinity. We will accordingly suppose that any one 
of these functions can be represented by 


fi) = 5 + 40) or by f() = klogt + 40, 


where ¢(#) is under no restrictions save that it is every- 
where finite; so that its Fourier constants, when 7 is 
large, cannot exceed c/n. The Fourier constants for 
the first term are given by the real and imaginary 
parts of 


27 
k git 
do tP 


Let us write nt = x so that the limits are now o and 
2an; the expression then becomes 
2TN 
k er dx 


SO ied ee 


If p lies between o and 1 the integral is finite when 
m is made oc in the upper limit, so we see that the 
EAS 
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Fourier constants are of the order of magnitude not 
greater than c/nt~? when n is large. Compared with 
this, the constants due to the finite function (#) are 
ultimately negligible. 

Hence we see that the order of magnitude of the 
higher Fourier constants depends on the nature of 
the infinities of f(z) if it has any, but if not, it depends 
upon the nature of its discontinuities in magnitude ; 
if it has neither, they depend upon the infinities or 
discontinuities of f’(t) and so on. 

It should be noted that if p<1 the integral 


2m 
\ | {(d| dt is finite; for then 
Oo 


i_ ky o k(20)1~P 
fo) iP Lic 9 


We thus see that whenever the condition that 
Pigs 
\ | f(#)|dz is finite is satisfied, the Fourier constants 
(@) 


decrease indefinitely with n, since both these conditions 
are satisfied at the same time, viz., when 7 lies between 
oand tI. 


The Convergency of a Fourier Series. 


a 


The convergency of a Fourier series, representing 
any function /(#), say the series 
a, COS t+ ad, COS 2% +... ., 
only needs to be considered when the series 


|| + |42| + |a@s| + 
is divergent, for if this series is convergent the first 
series is obviously convergent. 
It follows from this, and from what we have said 
above, that it is only the effect of the infinities and 
discontinuities of magnitude of f(#) on the Fourier 
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constants that need to be considered here; for 
discontinuities of slope or curvature only contribute 
a term of the order ¢/n? or c/n to the Fourier constants 
which could not affect the convergency in any way. 
Let the function f(z) be represented as the sum of 
several functions each of which has only one point 
of infinity or of discontinuity of magnitude. When 
any of these functions are being considered let this 
point be chosen as the origin. 

Now any periodic function of period 22, which 
has no infinities, and only one discontinuity, and 
that at the origin, can only differ from a multiple 


of the odd function defined by /(# =" from 


t= 0 to z= 22, by a function whose Fourier series 
is absolutely convergent and therefore need not be 
considered. But the former function is represented 
by the series 

si site ze 


+ 


I Z 


(see page 42), and this series, we have seen above 
(page 120), is convergent for all values except for t= 0 
(or a multiple of 27) when it is indeterminate and can 
assume all values from — 1-8519 to + 1°8519. Since 
the function at t = 0 is indeterminate within the range 


sin 3¢ 


“- Zs 


7 5 4 
from —~ to + -, the range of indeterminateness of 
2 2 


the series is 1°179 times the magnitude of the dis- 
continuity of the function. The formal sum of the 
series is zero, which is the point midway between the 
extreme values of the function at its discontinuity, 
and this property will obviously be preserved if this 
series is added to another series which represents any 
other function with no discontinuity or infinity at t= 0. 

Similarly, if the function has a single infinity at 
t = o, and in the neighbourhood of ¢ = 0 approximates 
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to k/t?, where p lies between o and 1, we have seen 
that the series representing it approximates, in its 
higher terms, to 


COSGi mr COS hay eCOsiat 


or to 


p[citieeen SIDE Ze sine 3 

\y1-2 r GLA? a giSe nm a I 
where ¢ and c’ are constants, according to whether 
the function is odd or even: if it is neither, it can be 
represented by the sum of both these series. It is, 
of course, only the behaviour of the higher terms of a 
series which affects its convergency. Now both these 
series, we have seen, are convergent for all values of 
Esave t = 0 (page 120). 

When # = 0, the first series, we saw, was infinite; 
while the second was indeterminate over an infinite 
range, which is exactly how the graph representing 
the function behaves in the two cases. 

We have thus seen that for any function, with any 
finite number of discontinuities, or infinities which 
are such that the integral 


J ifolae 


is finite, the Fourier series 75 convergent and can thus 
be used to calculate the function. 

We may note here that it will easily be found that 
the series 


cos ¢ COs 22 Cos 3t 
fp eG ee ae 
I 2 2 
which we have not mentioned represents the artificial 
periodic function defined by f(t) = log |scosec¢|. 
When fis small this becomes — log | z| approximately ; 
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which is infinite when ¢ = 0, but is a lower order of 
infinity than is represented by 1/#?. 

Since a periodic function (whether even or odd) 
which has a single infinity of the order k/i? at t = 0, 
where p lies between o and 1, has its higher Fourier 
constants approximating to c/n"~? we see, by integrat- 
ing both the function and the series (s — 1) times, 
that a series in which the higher Fourier constants 
approximate to c/n’? represents a function whose 
(s — 1)th differential coefficient has an infinity of the 
order k/t? att =o. Also from the foregoing results, 
we see that if the higher Fourier constants approximate 
to c/n’ where s is an integer, then the (s — 1)th 
differential coefficient will have an ordinary discon- 
tinuity in magnitude if it is an odd function and will 
have an infinity of the nature of log|z| if it is an 
even function. 

We may remark that the functions in the last 
paragraph are not necessarily artificial functions 
though they usually are. For instance the function 
defined by y= sin(¢+ 4) is an analytic function. 
It represents an odd periodic function of period 2z ; 
a distorted sine wave which has an infinite slope at 
the origin and attains its maximum value of unity 
when t= 2/2 —1. Near the origin we have t = y3/6 


tl 
approximately, or y = 61/371/3 s from which = = ft-?/3, 


so that the higher Fourier constants of the differential 
coefficient approximate to c/n'/3; and so the higher 
Fourier constants of the function itself approximate to 
c/ntls, 

The Fourier constants may, of course, decrease 
much more rapidly than is represented by the law 
c/nS; an and by may, for instance, when 1 is large 
approximate to ck” where k <1, as in Problem I 
of Chapter IV. In this case, if we differentiate the 
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series s times, a, and bd» for the differentiated series 
approximate to ck”n‘; and the ratio of one coefficient 
: : nN iS C : Hb 
to the previous one is (22) , which, in the limit, 
n 

when 7 is infinite, is &, so that the series is still con- 
vergent. Hence such a series represents a function 
in which none of the differential coefficients are ever 
infinite or have any discontinuities. Such a function 
is necessarily an analytic function.* 


Proof of Fourier’s Theorem by Direct Summation of the 


Series. 


We have proved in Chapter II that if any periodic 


function of period 2x can be expanded in the series 
Gee gcse Gd Ob, SU eee ee 
2 


then the coefficients must have the values there found 
for them. We also gave later a proof that the expan- 
sion was possible (page 23). Just how far the reader 
will regard these two proofs taken together as constitu- 
ting a rigorous proof of Fourier’s Theorem will depend 
upon what we may call his mathematical temperament. 

Another proof of Fourier’s ‘Theorem has virtually 


been obtained in Chapter V. We there obtained a 


* ‘The student familiar with the elements of the Theory of Functions 
will recognise that the Fourier series in this case will be convergent 
for complex values of ¢ as well as real ones. If we write t= u-} iv 
pees CE EE G ame) which is the ratio of one harmonic to the pre- 

ay, COS nt 
ceding, we find that when 7 is large it rapidly approaches a,4yel?!/a,,. 
Hence the series will be convergent so long as el¥| < k, 1.e. so long 
as |v | < log 1/k. 

If the series of Fourier coefficients is still more rapidly convergent, 
it is possible for the Fourier series to be convergent for all complex 


values of ¢. Such will be found to be the case when a,, or 5, approxi- 
nN 


kn ‘ : 
mate to ia even when & is greater than unity. 
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sine and cosine series which passed through the tops 
of a number of equidistant ordinates and we showed 
(page 78) that when the number of ordinates increased 
indefinitely, the coefficient of each sine and cosine 
term moved to a definite limit, and so the finite 
number of terms became in the limit the infinite 
Fourier series. ‘The fact that the coefficient of each 
term approximates to a definite finite limit as the 
number of ordinates is made infinite, taken in con- 
junction with the convergency of the series, which 
we have proved above, is a satisfactory proof of the 
validity of Fourier’s expansion. 

If we took any infinite series of different functions, 


Pib,-oo(1) . . . . we could of course determine the 
constants in the expression 
61Pi(t) + Cobalt) + . « » + enbn(?) (8) 


so as to make it pass through the tops of any # given 
ordinates; but it would not in general happen that 
the c’s tended to any definite limit as » was increased 
to infinity. Nor would it generally happen that the 
graph represented by (8) would, when the ordinates 
were sufficiently close together, approximate to the 
straight lines joining the tops of successive ordinates. 
In general, the curve would make a half oscillation 
in running from the top of one ordinate to the next ; 
a fact which is consistent with, and explains, the fact 
that the expression (8) has no definite limit when x 
tends to infinity. The fact that the sine-and-cosine 
series on the contrary approaches a limit, means that 
the graph represented by it approaches a limit, and 
this limit cannot be anything else than the limit of 
the straight lines joining the tops of consecutive 
ordinates. 

We will now sum directly the whole of the first 
harmonics of Fourier’s expansion including the 
constant term. We have 
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fi) = 2+ acostt ... 


+b, sint+ .. 
27 
where Ayes { J(24) cos nt dt 
ce 10) 
ie ee 
and bn = ai F(2) sin nt dt. 
ut 10) 


Let us use the letter x for the variable in these 
integrals and, since the starting-point of the range is 
immaterial so long as the range is 27, make the limits 
t—zaandt-+<2; so that 


I t+n 
an = — \ f(*) cos nxdx. 
Uw Ji_x 

The apparent dependence of a, on ¢ in this result is 
illusory, for, as we have just said, the value of z does not 
affect the integral. 

Replacing bn by a similar expression and adding 
together the two terms constituting the mth harmonic, 
we have 


1 ttn 
an cos nt + by sin nt = — \ f(x) cos n(x — t)dx. 
: 
t—T 
On writing « = u + ¢ this becomes 


=: i T(t + u) cos nu du. 


Hence the sum, S», of the first ~ harmonics plus 
the constant term is 


EV y+ w{2 + cos w+ cos re ee + cos nud. 


aed V) 
Replacing each cosine term by its exponential value, 


we easily find that the value of the series is see 
2 sin Su 
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Hence we have 


= = \ ie + 2) sa (ies) du. 


sin 4u 


Now change the variable again by writing (v+4)u=v, 
and for brevity write n’ for n+ 4. We then get 
I P(e 
Ss 


santo tc 


pa zy na iced: sane) 


Uv 
2n sin 
2n' 


ifn large, f(: + “7 , for moderate values of uv 

say for v <V, is sensibly equal to f(t) unless f’(z) be 
Ae 

very great ; also 2n’ sin Sat then sensibly equal to »v, 


hence, taken over for the range from — V to V, the 
integral is sensibly equal to 


(Os ae. Sa) 
I 
3 -V 

Now this converges when V = «, the contribution 
from large values of uv being really less than seems to be 
the case, since each positive half-wave is very nearly 
neutralised by the following negative half, when v is 
igige.e lt Y is moderately large, then, the value of; 
this integral is sensibly that of 


“sin v 
S 
\ = Av), 


which we have seen (page 11) is x. There is clearly 
no difficulty in choosing a large finite value of V which 
will make the integral in (10) as near to a as we choose 
and then in taking n»’ much larger than V so that the 
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approximations (of 27’ sin fe v and f(: a :.) ry) 


in (9) are as accurate as we wish them. 

We conclude, then, that the limit of S,, when 1 1s 
infinite, is f(2) provided f’(#) is finite. ‘This furnishes 
a rigorous proof of Fourier’s expansion. 

If the function f(z) has a discontinuity in magnitude 
at the particular value of ¢ considered, the result (9) 
breaks up into two, giving 


nt “ 
I v sin v 
Sn *{ (: ) dv 
aa jt ty! v 
(@) 


2n sin —; 
2n 


0) e 
I v sin v 
eS es dv. 
al C= 5) a 
—nNTN 


20 
2n 


As previously, the first integral becomes in the limit 


F(t + 9) 


Lh 


{ S88 = 4/6 + 0) 


fe) 


and the second integral becomes likewise 


PUES 


IU 


Ome 
(ie ee, 
v 
—“ 
so that the sum to infinity of the Fourier series is 
Hitt 0) + ft — 9} 


at a point of discontinuity, and f(#) at other points.* 


* The student with a little mathematical courage will probably 
obtain this result by boldly putting 2’ = oc in (9) when it becomes 


S.= 2] 70 BY do=* fey, 


and this is not so reprehensible as some authors would have us believe. 
It is necessary, however, before writing oc, or any other particular 
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At a point of discontinuity this investigation clearly 
obtains what we have called the formal sum of the 
series; and we have seen before that at such points 
this sum of the series is midway between f(t + 0) 
and /(z — o).* 


value for a variable, to make sure that the expression assumes a single 
definite value when this is done. A case in point is the series 


sin 1/f-++- $ sin 2/44 sin 3/t+ .... 
As t increases towards +- oc, this series, as we have seen, approaches 


. . ri . . . 
the limit of + —; while when ¢ is negative and moves towards — oc 
2 


the series approaches the limit of ae, while actually at t= +c the 


series is indeterminate and can assume any value from — 1°8519 to 
+ 1°8519. Such behaviour, however, is not difficult to detect, and 
in the case of infinite series, is always bound up with the fact that at 
such a point the series does not begin to converge till after an infinite 
number of terms. 

* It is clear also that the limit of (9) when n’ is oc will only be 


tf 22 do or £0) 
id es Rae 
provided that f(¢ + 2) is sensibly constant for all values of v below a 
n 


reasonably large quantity, say V, when n’ is very large. This is 
generally expressed by saying that f(¢) must not have an infinite number 
of maxima and minima or discontinuities in a finite space in the 


neighbourhood of the value of ¢ considered; for if so, f(e+ a) 
n 


could fluctuate in magnitude during a finite change in v when n’ 
was infinite: and so we could tell nothing about the value of the 


integral 
| f(¢ 4 =) sin Y dv. 
n Vv 


This is the reason why so many books dealing with Fourier’s series 
continually repeat the condition that the function must not have an 
infinite number of maxima or minima. We have generally omitted 
specifying this condition, since no practical function ever does behave 
insuch a manner. Such behaviour is exclusively confined to functions 
invented by mathematicians for the sake of causing trouble. 
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EXAMPLES. 
By putting t= o in the formula (21) of page 30, and by 
eee the a’s by their integral values, prove that 


fo +42) +4(8) +... +4) = Ff "pome 


a = [10 cos ntdt + *[ 10 cos 2nidt-+ .. . 


[The reason for writing }f(o) + 4/f(27m) instead of ether f(o) or f(27) 
is in case the periodic function f(t) has a discontinuity at t= o which 
will of course be the case if f(2m) += _f(0).] 

Be eke by repeated integration by parts, that 


uC cos ntdt = — ae) ees Sr 


Dees that all the differential coefficients are finite and the series 
is convergent. 
3. Prove, by substituting the result of Ex. 2 in Ex. 1, that 


“hod = “(aso + 5(F) +». + User) 


2 fo) I I I 
fen) — oe a 
"f"en) — FON, 

nt Ve 4 ke 34 [ss 

The numerical series in brackets are known to have the values 
nm nt 7S 78 10 
6° 90° 945° 9450° 93555) cae 

4. By changing the range of integration in Ex. 3 from 2m to nh 
prove that 


—— 


+ 2 


. respectively. 


nb 
SOat = HII + SOF ICD + «+ fod) 
Bes) n ae \ Sess We, mr 
he ro) irs —'"(@)} 
Vv Vv 38 VII /, Vil 
5M FO} + Elsen) — FH} 
cpt 
Lea 


provided the series is convergent. ‘The series can often be used when 
it is not convergent, provided the terms decrease to a small value before 
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they finally increase, for it may be proved that in many cases the 
error made by stopping at any term is less than that term. 

5. By taking the limits of integration as « and oc prove that the 
result of Ex. 4 becomes 


[ fod = safe) + fet + fet 2) 4.64 


Sh, Patch ane sie! beats 
37 ~~ eyef Coe cranes Co ae 


|7 


oe te 


} 
ie 
supposing the integral convergent. 
: ile cee : : 
6. By taking f(‘) = Zin Ex. 4 and the limits of integration as m 


and , which are integers, prove that 


I I I 

Se meaner mes 
ee ie me Es ue i) 
12\n2 m2) ' r20\n* m*/ — 252\n® ~— m8 


ft I 
a sae a = Put 
Calculate log,2 to seven decimal places by taking = 10 and 
m==20. [Ans, *6931472.] 


7. Prove, from the result of Ex. 6, that the limit when m is infinite 


m I I 
log — = — : 


n 2n N 


of “ — : eo eee + —— log m is a constant. Show, by 
taking = 4, that this constant is 577215 approximately. 
8. By taking f(t) = log t and the limits m and n as integers, show 
that 
m log m/e —n log nle =} logn+ log (n+1)+... 
BAI eet \ eles hl Sh ae 1a 
ee rar a nJ " 360 \m3 8 J 


12 \m 


Le fat I 
lg ar 
Regarding m as a variable and n as a constant, show that this result 
can be written 
log | m= Const. + (m+ 3) log m —m-- —- — sien 
I I I 
1260m®  1680m? ' 1188m® 


=I 


By taking m= 10, show that the constant is ‘91893856... It 
can be shown theoretically that this constant is } log, 27. 
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Deduce from the last result, that when m is large, | m approximates 


to a ratio of equality with 
(2) 
Vanm\ =) 


g. From Ex. 5, page 73, obtain the following results 


bd ; sin2x” . sin4x , sin 6x 
—| log sin « dx = x log 24 { 2 + “ f- Pe +- tay 
(eo) 
wheno< x <7; 
bs ; : ay 
| log cos « dx = x log 2 — (Be SEY an oh + 
° 1 2 3 
Tv Tv 
when — —< 4 < — 
2 par; 
and 
— [log tan x dv = 8" ! gh Oe neeinee om 
° I 3 5 


Tv 
when OI 


10. Prove, by differentiation or by replacing the sines and cosines 
by their power series, that the curve represented by the series 
cos B 


- 4 COVA tn Con Bed hae 


2P Be 


is finite when 8 = oif p > 1 and hasa finite radius of curvature at 
8 =oif p> 3. Also prove that the series 
sin , sin 2B 


a p+. 


has a finite slope at the origin if p > 2 and an ordinary point of 
inflection at the origin if p > 4. 
11. If (x) = o when x = 0, prove that the series, 


p(x) , (2x) , 9(3-) 
pee eee sea ee 


can have any value when w = 0 which the integral 


[a 


o « 


can assume, where X is any positive quantity. 


CHAPTER. VI 
FOURIER’S INTEGRAL THEOREM * 


Fourigr’s theorem, for a periodic function f() of 
period T, is that f(#) can be expanded in the form 


(see page 29) 


1G). = S + a, cos ia + d, Cos ssa Ae 
+ b sin FE + by sin FE +. ele ears Li) 
where 
Te 
2 2nit 
an = -ra\ f(t) cos dt . . . (2 
I \ 9 L @) 
T/2 1 
and by = at f(z) sin oie Ai tee oe 
iba ) at (3) 
If we change the variable under the sign of integra- 
tion from ¢ to 4, we may multiply (2) by cos pei and 
(3) by sin sie and put these terms under the 


sign of integration; since they are mere constants 
during the integration with respect to 4. The two 


2nmt » Piling 
terms dp Cos ar + Dy sin 


T constituting the nth 


* This chapter is intended for the student of physics; it has no 
engineering applications; but we hope it will not on that account be 
found uninteresting to all engineering students, 
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harmonic then become 


211 co 


—T/2 


, 2nd — bay, 


: . 270 . 
We will write <= = 6 and now suppose T increases 


48 
without limit so that « becomes indefinitely small ; 
the above expression then becomes 


avo cos no(i — t)da. 


—t/o 


Now replace all the harmonics in (1) by expressions 
of this a ; we get 


ji Oe a Ef(A)da spe “hay cos o(A— t)dh 


— 16 ae 
+ rh K) cos 20(4 —f)dd + . 
Now we know that if F(«) is any function of , the 
limit of of F(o) + F(o) + F(2c) + ... + F(x)} when 


x 
C= Onis \ F(a)dx; and, stmce each single termes 
oO 


indefinitely small in the limit, it makes no difference 
whether the first term is oF(o) or $ oF(o). Applying 
this to the above case, we see that if we write 6 = da * 
and replace the limits +a2/o by +, which they 
then become, we have 


Mo wh af 0 cos «(4 — f)ddA . . (4) 


This result is known as Fourier’s Integral Theorem. 


* This is quite a different use for « than that in the earlier part of 
this book : we do it for the sake of uniformity with other writers. 


r FOURIER’S THEOREM 145 
The theorem, of course, only holds provided the 
integrals are convergent. The integral 


\ {0 COsta(An— stdin i x) (5) 
= 
will necessarily be convergent provided 
{ata 
—o 


is convergent. We will accordingly suppose that this 
condition is satisfied. ‘This requires that f(A) should 
decrease faster than 1/4 when 2 is large, for both 
positive and negative values. ‘The student will also see, 
in a similar manner to that adopted in the last chapter, 
that if f(A) is finite everywhere, but has discontinuities 
in magnitude, each integral with respect to 4 in the 
equation 


fc) = #\ 


CC co 

[cos zal T(A) cos ad di 
fe) —% ’ 
f(2) sin a a ae 


ee 


+ sin zai 


in which form (4) may clearly be written, 
will tend to proportionality with 1/« when « is large ; 
and if f(A) has discontinuities in slope but not in 
magnitude, each integral will likewise tend to propor- 
tionality with 1/7, and so on. The integration with 


ice) 
respect to « is always convergent when \ | fA) | da is 
—%4 


convergent owing to the cos «t or sin at factors, which 
virtually make the integral a series whose terms are 
alternately positive and negative such that, when « is 
large, each term is less than the one before it. 

L 
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Like Fourier’s series, Fourier’s Integral ‘Theorem 
applies both to analytic or artificial functions, subject 
to the integral 


) 1x | da 


being convergent. 

We can have functions of two distinct types; either 
functions gradually increasing from zero at t= — 
and decaying to zero again as ¢ tends t+ © ; or those 
which are zero from t= — © to t=4,, say, when they 
suddenly become finite and remain finite from t= 4 
to t = ft, say, when they suddenly become zero again 
and remain zero up to t = o. A function of either of 
these types is conveniently called a “ disturbance ”’ or 
a “ Pulse.” It is important to bear in mind that it 
must necessarily be defined from t= — ~tot=+ 0 
or Fourier’s Integral Theorem cannot apply to it. If 
we take the limits of 2 as ¢, and #, in Fourier’s Integral 
Theorem, instead of — o and + o, we are simply 
applying the theorem to the pulse which is equal 
to f(t) from ¢=7, to %="% and zero-outsiae mt uese 
limits. 

Fourier’s integral theorem takes a simpler form 
when /(#) is either an odd or an even function, as might 
have been expected from the behaviour of Fourier’s 
series. If, for instance, f(#) is an odd function of 1, 
the first integral with respect to 4 in (6) vanishes ; 
since f(A) is an odd and cosad an even function of 4 
and the product is odd: also f(A) sin ad being an even 
function, the integral from — « to +cc is equal to 
twice the integral from 0 to cc. 

- Hence, if f(z) is an odd function of t, we have 


Oe 2 {sin at] {4 sin al al day 
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Similarly, if f(2) is an even function of ¢ we have 


fi) = 


If we write 


2 { cos at] { cos ad ai Atop e an Ss) 
{ £0) sin at di = 4a), 
and 


\ 4 cos ad di = y(a), 


these results may be written in the following form : 


that if YF) sin at di = $0 i ae 
hes [so sin ate = 2/0 | ee 
and if \ © ee) | 

ine [ve cosatde =) | Ye 


In the first of these f(A) is odd, and in the second it 
is even; and, since we see that the relation between f 


and ¢ and f and y is reciprocal, (save for the factor =) 


we see that, if we wish to interpret > and y for negative 
values of a, ¢ is to be taken as odd and yw as even. 

Asean example, let f(t) = ¢-¢ from’ ¢ =e to ¢ + ec, 
Then we have by (9) 


ee sin at di = Ko), 


fe) 


but this integral is ——,, ; 
§ cz La? 4 
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whence we infer that 


Ce 
a . U 
+——5 sinatda = —e-%; 
gE 2 , 
\. co +a 2 
or, changing symbols, that 


SEDI 


~ ~ sin mx n 
mney 
Cab 2 


Similarly by (10), since 


= C 
e-cosad dA = =; 
\. Fae 


we have 
oe 
C 7 
Tey COS OL dd el, 
ce? +a 2 


or, in other symbols, 


oC 
\ COST 7 cm. 


Sieh me Shyer 
Fourier’s integral theorem in the forms (g) and (10) 
thus becomes a useful theorem in the integral calculus. 
The way in which Fourier’s integral theorem per- 
forms the “ harmonic analysis” of a pulse can be 
seen by the following considerations. ‘The ratio of the 


period of the mth harmonic to the (nm + 1)th is Psa 
n 
which tends to unity when 1 is infinite. Also, when the 


period is made infinite, all harmonics of finite period are 
necessarily harmonics of an infinite order; they are 
thus infinitely close together, or in other words, there 
are an infinite number of harmonics whose periods lie 


between two given finite limits. Now the equation 
(6) can be written 


> 
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iE) = a ae T(4 cos ad a cos at da 


ce 


IU 


ate a3 [ wo sin ad a sin at da 


Dropping the integration with respect to a, we see that 


I Cc 
- [\ F(A) cos ad a] cos at da 
—o 


I 


ae Tl) sin ad ai SIN aida: Ae Gr) 


is the part of f(#) which consists of harmonics whose 


270 


periods lie between — and , or whose fre- 
a a 


If we denote 


; é a 
quencies lie between — and 
2 


the integrals within square brackets by P and Q 
respectively, (11) can be written 


(=) cos at + eS sin at. 
TU TU 


Here a is not quite constant but varies over the 
small range da from a to a + da, so we see that the 
amplitude of the cos at term at a may be regarded as P/x 
per unit range of a or 2P per unit range of frequency. 
Similarly, the amplitude of the sin at term at a may 
be regarded as Q/z per unit range of a, or 2Q per unit 
range of frequency; and the amplitude of the com- 
plete harmonic as 2./P? + Q® per unit range of 
frequency. 

It will now, we think, be best if we adopt the 
language of that branch of physics in which this sub- 
ject has its chief applications. ‘The pulse is generally 
supposed to be some disturbance travelling through 
a medium ; it may be, say, an air pulse caused by an 
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» 
explosion, or an electromagnetic pulse due to the 
stoppage of an electron by collision. In the first place 
we should take f(#) as proportional to the condensation 
in the medium at any given point; and in the second 
case to the electric (or magnetic) field at the point 
considered. In either case the energy per unit volume 
of the medium at time fis proportional to f(2)?, and sowe 


ec 
may take | frat as a measure of the energy of the 
—% 


whole pulse.* 

Now on page 34, in finding an expression for the 
mean square of the value of a periodic function in terms 
of the coefficients, we were really expressing the fact 
that the energy in the function was equal to the sum 
of the energies in each of its harmonics. An exactly 
similar result holds in the present case, which we will 
now obtain. ‘Taking the equation 


CC eg 
je =| aa{ F(A) (cos ad cos at + sin «A sin at)da 
oO TPES, 
or, with our previous notation, 
Ce 
Har “| da(P cos at + Q sin az) ; 
[@) 


let us multiply both sides by f(#)dt and integrate from 
—« to-++«. Remembering that P and Q are func- 
tions of a only, we have 


j seo bs z| a [ef cos at dé +. of Ke sin at at] 


a | fore iS | @ + Ode. 


* ‘There is no need to introduce the factor 4, which the student may 
be tempted to supply, as we have dropped other constant factors 
depending on the properties of the medium. 
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¥ 

This shows that the total energy may be divided 
up among the different values of «, the energy lying 

2 2 

between « and « + d« being mn "ae that, 1s, the 
P? + Q? 
== Ole tlic 
TU 


energy per unit range of « at « is 
energy per unit range of frequency at the frequency 
ies i . . . 

3 8 2(P2 + Q?). Since the energy of an oscillation 


is proportional to the square of the amplitude, this 
result is consistent with the one arrived at on p- 149. 


When the different frequencies in a pulse are separ- 
ated in any manner the result is called the “ spectrum ” 
of the pulse. The light given out by a piece of red- 
hot metal or carbon, for instance, consists of the 
“pulses ’ produced by the collision of the electrons 
with one another or with the atoms of the substance. 
There is not necessarily any periodic motion whatever 
associated with the origin of the light from an incan- 
descent solid body. A glass prism or a diffraction 
grating analyses this succession of pulses into its 
constituent harmonics in exactly the same manner as 
Fourier’s integral theorem analyses it. The fact that 
the spectrum, in this case, shows the presence of 
harmonics of all frequencies within the range of the 
instrument, shows not only that the origin of the 
light is in pulses, each containing all frequencies in 
itself, but that these succeed one another at irregular 
intervals. For a succession of exactly similar pulses 
at exactly equal intervals of time would obviously 
constitute a periodic disturbance, from which all 
harmonics would be absent except those whose periods 
were sub-multiples of the given interval; hence the 
spectrum would be a “bright line” spectrum and 
not a continuous spectrum. The distribution of 
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energy in a succession of similar pulses can only be 
the same as the distribution in each pulse if they 
follow one another at perfectly random intervals of 
time. 

If we make a slit in the screen on which the spectrum 
is focussed we can let through any desired small range 
of frequency at any point, and by absorbing these rays 
we can measure the energy in this small region of the 
spectrum. In this manner we can get experimentally 
the distribution of energy in the spectrum.” 

We have seen above that if we knew the form of the 
pulses we could calculate this distribution of energy. 
We shall now show that, knowing the distribution of 
energy, we can find the shape of an odd pulse, and 
also of an even one, which would give such assigned 
distribution of energy in the spectrum. f 

Let {F(«)}#dx be the energy in the spectrum be- 
tween «anda+ da. Then we have aF(a)? = P? + Q?. 
Accordingly, if the pulse is even so that Q = 0, we 


have: \/#i (a) =| P= AO cos «4 di, where f(#) is 


the form of the pulse. 
Equations (10) at once give us the result that 


| 27) cos at dt = 7 f(i) 


I 


oe joe val FO) bbs ee 


* It is, of course, necessary to correct the result for the light lost 
by reflection or absorption in the prism or grating. 

} The student familiar with radiation from electrons will recognise 
that if an electron is projected against another and is brought to rest 
by the interaction and repelled back along the line of advance, the 
pulse will be an even one; while if it had passed on without being 
brought to rest the pulse would have been an odd one. This is easily 
seen, since the disturbance at a great distance is proportional to the 
acceleration of the electron. 
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wr 
Similarly, for an odd pulse its form would be given 
by 
me 


fd) = Jah FO See: 


Example. Word Rayleigh proposed the expression C%x2e~ 2d for 
a+ da 
, the 
2 


constants C and ¢ depending upon the temperature of the body 
emitting the light. Here, F(a) = Cae-, so that an even pulse to 
give this distribution would be given by, 


, Shed 
the energy in the spectrum between the frequencies — and 
2m 


G a 
J) Jal oo cos at de 


and an odd pulse by, 
Bd as : 
Ai) = val @ sin at da. 


Multiplying the second by / — 1 and adding it to the first we get 
ei Gath C Cit + 02 
5 a (it-c)ady — = en ees 
Fo(2) at tf, (2) Jal oe V rit = c)? V/ 1 (#2 = ote 
And so 


Ce — 7?) 
na VJ 1 (c2 + 22 


F(t) 


and 
2Ccr 


A® ae a/ 7 (02 a 12)2" 


Similarly, we could find the shape of the pulses to produce any other 
law of distribution of energy in the spectrum, such as Planck’s, whose 
formula agrees much better with experimental determinations than 
the one used above; but this calculation has not any very great 
physical value, since white light does not consist of a succession of 
pulses of exactly the same shape, so what we have arrived at is only a 
a sort of mean shape of the pulses. 
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EXAMPLES. 


Tv 
Taek | F(x) cos nxdx = ¢(n) 
° 


prove that 


[ cos nxdx = 4e(0) + (1) cosa + 


2 @(2) cos 2n-+ @(3) cos3m-+ ... 
providedo <n<n. 


[For within this range each side is equal to ~F(n)] 


2. Prove, similarly, that if 


[Fe sin nxdx = y(n) 


then 


[ve sin nxdx = y(1) sin 2 + (2) sin 2n-+ (3) sin 37+ . 


provided thato <a <7. 
3. Verify the result in Example 1 when F(x) = nr — x. 
4. Prove that the energy, lying between the frequencies «/27 


and («+ da)/2x, in the pulse given by f(t) = 


to e262, 


5. Prove, similarly, that for the pulse f(#) = e~/*", the distribution 
of energy is proportional to e~"/2*, 
6. Prove, similarly, that for the pulse f(t) = sech pt, the distribu- 


: : : Te 
tion of energy is proportional to sech? — 


¢ ‘ : 
———, is proportional 
py arg ee Oe 3 Y 


7. Show that the distribution of energy in the pulse, which is 
equal to unity from ¢=— 1 to t=1 and to zero outside these 
toute es : sin «\2 
limits, is proportional to 5 

Ci 

8. Show that the distribution of energy in the pulse, which is 
zero when ¢ is negative and equal to ¢-*sin pt when ¢ is positive 
(where & is very small compared with 7), is very approximately given 
by 


I 
2+ (au — p)® 
Hence note that a damped oscillation, when examined in a perfect 
spectroscope, will not show a perfectly sharp line, but will show a 
maximum of energy when « = #, and the energy will be appreciable 
so long as « — p is not many times greater than &. 


FOURIER’S THEOREM 155 


e 
g. Show that the forms of the odd and even pulses, which will 
give a constant distribution of energy over the frequencies between 


(« — 6)/2m and (« + 6)/2m, are 
j= oS sin at 
and (Oe ie of) cos at 


respectively. 

10. The energy lying between the frequencies «/2m and («+ da)/27 
is proportional! to («?/~2¢-%«)dx ; prove that the forms of the odd and 
even pulses to give this distribution of energy are 

__ sin (p tan”? z/b) 
iG) sand (B2 es #2)p/2 


and 
__ cos (p tan™1 ¢/b) 
SO= ep PPR 
respectively. 
it. If f(z) represents any convergent pulse, prove that the succes- 
‘ tek P y BEneP Pp 
sion of pulses 


Efe — am) +fO+fe+ amt fet a+... 


is represented by the Fourier series 
“2 G@q_COSt--. . 
+b,sint+... 


where 


ip ies 
On = | FS cos nt dt 
-—o 


ih GS 
mid he teahe ft) sin nt dt. 
UF Mace 


CHAPTER VIII 


THE SEARCH FOR PERIODICITIES IN OBSERVED DATA NOT 
PERIODIC 


Tue physicist, astronomer, or meteorologist— 
especially the latter two—frequently obtain curves 
which are suspected of containing one or more periodic 
components. If the curve is found to be definitely 
periodic, and the observations have extended over a 
sufficient number of complete periods to establish 
the fact beyond doubt, recourse can be had to the 
methods of harmonic analysis previously given when 
as accurate an analysis of the curve can be obtained 
as its own accuracy permits of. More generally, how- 
ever, the curve contains several harmonics the periods 
of which are not in simple ratios to one another, so 
that the period within which the motion recurs is very 
large—if, in fact, it recurs at all. For instance, a 
curve composed of two harmonics of periods 12 and 
17 days has a period of 204 days, but if the ratio of the 


periods, instead of being 1-416666 . . . , was in error 
by about 1 part in 760, so that the ratio of the periods 
was actually 1-414213... or »/2, the “period,” 


instead of being 204 days, would be infinite; since 
there is no common multiple of two numbers bearing 
this ratio to one another less than infinity. Thus to 
calculate a period in cases like this, where the periods 
of some of the harmonics are only known to a limited 
degree of accuracy, is to tread on far too uncertain a 
ground, 

We have thus to determine the presence of a har- 
monic constituent in the curve, or set of observations, 

156 
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in the absence of any indications of a definite period. 
This can be done by the application of Fourier’s 
integral theorem explained in the previous chapter. 
Let the given curve or observations be represented 
by f(i) from ¢ = — T to ¢ = T, and let us apply the 
theorem to the pulse defined to be equal to f(z) from 
t= — T to z= T and to be zero everywhere outside 
this range. Then we see from page 151 that the 
energy in this pulse lying between the frequencies 
a a+da. i 
--and ———_ 16 proportional to 
27 27 


[psa pseal|e- 0 


or, say, R@da. ‘These integrals can be evaluated by 
Simpson’s rule or other suitable method of quad- 
ratures; such evaluation must be done for a fairly 
large number of values of a. 

If R? is plotted against a the result is called the 
“* Periodogram ” of the original curve or set of obser- 
vations. We thus see that the periodogram is nothing 
other than the distribution of energy in the spectrum 
of the pulse defined as being equal to the observations 
over the range for which they extend and zero outside 
that range. Practically it is advantageous to split 
any given function into its odd and even components ; 
R for either of these functions is then either 


2 1) COS CHAN « . <5) geet) 


or 0 4) SIL GA dA 1 eee eee en) 


It is further advantageous to plot |R|, instead of R, 
against a and to call these diagrams the periodograms 
of the functions in question. 
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We shall now show that the periodogram possesses 
pronounced maxima at each value of a which corre- 
sponds to a harmonic constituent in f(#) provided that 
the range of observations extend over several such 
periods. As an example, let us suppose f(#) consists 
of the harmonic ¢ sin pz only and that the observations 

Na 
extend from — Me tos-+ =— so that they covers 
periods. In this case the first integral in (1) vanishes 
and we have 


Nand 
R= ac\ sin pA sin ad dd 
e) 


es gate — a)Nx/p _ sin (p + ae . es 
p—a pra 

The first term attains its numerical maximum value 
of 2cNx/p when a = p, while the second term neces- 
sarily lies within the limits + 2¢/(p + a), and lies within 
the limits + c/p in the neighbourhood ofa = p; that is, 
it is only +1/2Na of the maximum of the first term. 
The first term vanishes when (p — a)Nz/p is a multiple 
of x other than zero: that is, it vanishes when 
a= (1 +1/N), ora = p(1 +2/N), etc. It will have 
subordinate numerical maxima at roughly half-way 
between these points; the value of the maxima when 


(a — p)[p = + (n + ¥)IN being 


oo 
(n + 4)p/N 
approximately, which is Cane: of the maximum 


when a = 9. 

Fig. 22 shows a graph of |R| when N= 5; the 
extent to which the values of |R| are uncertain, 
owing to the neglect of the second term in (2), is 
shown by the distances between the two dotted curves 
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and the line midway between them. It is clear that 
such a relatively small uncertainty can do nothing to 
mask the subordinate maxima, let alone the very 
pronounced central maximum. 


a 
ett ee ee 
wert mane mew wo mw www ee mew men ns 


10 . ° . 8 a iP 


Fic. 22. 


If we had considered the harmonic c cos p¢ instead of 
c sin pt, we should have had 


Nx/P 
Kis 2c\ cos pd cos ad dd ; 
ie) 


which would have given the same as (2) with merely 
the sign of the second term changed, so that the graph 
of {R| would still lie within exactly the same limits as 
previously. It should also be noted that all these 
results hold whether N is an integer or not. 

Let us now suppose that f(#) contains two harmonics 
of different periods, say, csin pt-+c’sinp’t. If the 
two periods p and p’ are widely separated there 
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will be no difficulty in recognising the two central * 
maxima corresponding to a= p anda=p’. If p 
and p’ are nearly equal, the two maxima become 
confused; but if the amplitudes of the two harmonics 
are at all comparable, separation of the periods is 
always quite easy when the two central peaks stand 
just clear of one another: that is when (p’ — p)/p = 
+2/N; or when the number of periods of p’ com- 
prised within the range of observations differs from 
the number of periods of p (which is N) by two. 
Hence to separate harmonics with nearly equal periods 
an extended series of observations is necessary. 

Separation is generally easy when the maximum 
of one central peak coincides with the zero of the other ; 
so that one of the harmonics has only one more period 
in the range of the observations than the other. In 
this case, if c and ¢’ are of the same sign and com- 
parable, the two central peaks will combine into a 
single peak whose width is three times that of the 
subordinate peaks. The subordinate peaks for the 
two periods will, in this case, reinforce one another 
and so will be prominent. If c’ = — c the peak for 
p would appear negative if we plotted R but since 
we plot | R| it will readily be seen that the combination 
of the two peaks will give two equal peaks each of 
a width 1-5 times the width of the subordinate peaks. 
In this case the two systems of subordinate peaks will 
largely neutralise one another and so may be indis- 
tinguishable ; nevertheless their width is known theor- 
etically (viz. a variation of 1/N in the ratio «/p), and 
the presence of two equal peaks each 1-5 times this 
theoretical width is conclusive proof of the presence 
of two harmonics of nearly equal periods and of 
amplitudes of opposite signs. 

In our search for periodicities by this method we 
are not interested in the exact shape of the graphs: 
we merely use them for the positions of their maxima : 
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it follows that we could use instead any other curve 
which possessed equally distinct maxima in the same 
positions. A much more easily obtained curve, and 
one which is quite as satisfactory, can be found by 
replacing the cosa/ and sina factors in the trouble- 
some integrals 


T ae 
\ F(A) cosad did and \ F(A) sin ad da 
fe) ie) 


by another periodic curve of the same period which 
will make the integrations easier. Such a curve is that 


Fic. 23. 


in Fig. 23, which is alternately + 1 for half a period 
and then — i for half a period. When the origin is 
chosen so that it is an odd function, its equation is 


$;(ad) = 4 {sin ah me is oe ai woe ihe | 


and when arranged as an even function, it is given by 
cos 3aA , cos Sad _ | 
3 5 J 


Substituting ¢,(a/) for sinad in the second of the 
given integrals gives us 


Pet... xf oe. @ 


° rele 2nila 
M 


dead) = {eos ah — 
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while the other one becomes 


n/2a e3n/2n e5n/2« a 
[f- f+ | ROT: + | ] Koa. . - (4) 
fo) m/20e oJ 3n/2a 
Each of these component integrals is very easily 
obtained, since we have the graph of f(4) (or can easily 
get it by plotting the observations) and from this we 
can read off the ordinates required in evaluating the 
integrals by Simpson’s rule. How the graphs of |R| 
are modified by using these functions ¢, and ¢, instead 
of the sine and cosine may easily be seen from their 
harmonic analysis. When using the function ¢,(aA) 
instead of sin af, the first term in (2) becomes 


4 (sin(p —a)Na/p , rsin(p—3a)Na/p , 
x| p—a 3 p—3a J 
instead of 


sin (p — a)Nx/p 
eer ee 
the second term of (2) being just as negligible in this 
case as in the former case, 

Thus the shape of the maximum in the neighbour- 
hood of a = p is unaltered, but we have introduced 
into the graph for |R| another maximum at a = 9/3 
which is 3rd of the height of the maximum at a = p 
and which need not correspond to any frequency in 
f(é) at all. Similarly, there is another maximum at 
a = p/5 which is {th of the maximum at a = p and 
which also does not exist in f(#). It follows, therefore, 
that we can use this graph just as the correct one, 
provided we neglect these maxima at sub-multiple 
frequencies. Should, however, the maximum at 
a = p/3 be much more than ird of the maximum 
at a = f, it is clear that a harmonic of the frequency 
corresponding to a = p/3 must exist in f(Z), since 
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e 
there is more of this frequency present than we have 
introduced by using the functions ¢, and 4,.* 

We may vary the functions ¢, and ¢,, so as to reduce 
still further the labour of evaluating the integrals 
involved, without affecting their effectiveness in 
detecting the various harmonic constituents. If we 
add unity to each of these functions, for instance, the 
result will only be to add to |R| the constant quantity 


ae 
\ f(4)di; and will therefore make no difference to 
—T 


the peaks of |R|; yet it will have reduced the labour 
of finding |R| to one half. Hence we may use for 
é, and ¢, functions which are zero for half a period 
and unity for the other half. It can readily be seen 
that we can go further in this direction and use 
functions which are zero for three quarters of a period 
and equai to unity for the other quarter. Going to 
the limit, we see that we can use functions which are 
zero throughout the whole period, save in the imme- 
diate neighbourhood of one point. In this case the 


trouble of integrating ¢$,(«4) . f(A) and 4,(ad) . f(A) is 


* The student of optics will readily see the optical interpretation 
of this. The first periodogram corresponds to the analysis of a pulse 
by a grating which only produces a spectrum of one order—the first. 
The second one corresponds to a grating which not only produces a 
first order spectrum but a third order one of $th the energy of the 
first one and a fifth order of ,!,th the energy and so on, the image 
in the third order underlying the image in the first order correspond- 
ing to 4rd the frequency or three times the wave-length. No 
diffraction grating, only giving one order of spectrum, has ever been 
able to be constructed, although it is theoretically possible. We 
have to be content with one behaving more or less in the second 
manner. There is, similarly, no reason why we should insist on 
employing the analogue of the first case in constructing the periodo- 
gram. ‘The student may also note that if the intensity of the image 
in the first case is unity, the intensities of the different order images 
in the second case are 8/n, 8/gn*, 8/25x%... which add up to 
unity. 
Wee 
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x 
done away with, as the result simply reduces to the 
sum of those values of (4) for which the function 
¢, or ¢, is not zero. This is made the basis of the 
practical method described in the next Section. 
From the above we deduce the following method 
for searching for periodicities in any given series of 
observations. We take the origin of time at the 
middle of the range and represent the observations as 
the sum of an odd and an even function by the formula 


fy Oa fe) ee 


For each of these functions we plot |R| and find its 
central maximum peaks. Each of these corresponds 
to a definite harmonic in the series of observations. 
In this manner all cosine harmonics are found from 
the even function and all sine harmonics from the odd 
function. 

Having found from the two graphs for |R| the 
frequencies that exist in f(#), we must next find their 
amplitudes and phases. Let a, be the value of a 
corresponding to any frequency present in f(/). Let 
M, be the central numerical maximum at a, of the 
graph for |R| constructed with the use of the cos ad 
factor for the even component of f(#), and let M; 
similarly denote the maximum for the odd component 
of f(¢) constructed in a similar manner with sin ad. 

Then if a, cos a,t + 0, sin ayt are the two correspond- 
ing harmonics in f(#), we have from (2), (page 158), 

Nis rae 
Oy 

But 2Nz/a, is the total range of time, T, over which 
the observations extend. Hence we have the simple 
relation 

Ca WVineis 
Similarly, a, = M,/T. 
If the graphs for |R| had been constructed with the 
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aid of the functions ¢,(a4) and ¢,(ad) and if M; and 


M, had been the values of their central maxima at 
a, respectively we should have had 

Gat NAN ee ue Ae aie BCS) 
and Oreo Nia Ta. > res 15, (6) 

We must similarly write down the harmonic terms 
corresponding to any other plainly indicated central 
maxima flanked by their subordinate maxima and 
minima which occur in either of the |R| graphs, 
neglecting, as we have explained above, central 
maxima corresponding to sub-multiple frequencies if 
the functions ¢,(a/) and ¢,(a4) have been used unless 
such maxima are unduly prominent. 

Having found in this manner a number of periodici- 
ties in f(#) and the amplitudes of the sine and cosine 
components of each of them, we may subtract the 
sum of all these terms from f(t), when we shall be able 
to see whether the residue consists of periodic terms 
whose periods were outside the region tested, or 
whether it consists of a non-periodic function, or 
whether it appears to be nothing more than the errors 
of the given observations. 


A Simple Way of Applying the Method in Practice. 


The foregoing method may often, in practice, be 
greatly simplified—so much so that it hardly seems to 
depend on the above Se att making use of the 
fact that if we add the values of a function for equi- 
distant intervals, all harmonics, save those for which 
the interval is a period, tend to cancel out.f 

Let us suppose that the given observations can be 


* Similar easily obtained results will hold when any other periodic 
functions are employed in place of 9, and 9p. 

+ If the accuracy of the given observations is poor they should be 
plotted, and smoothed values read off from the graph in the ordinary 
way. 
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represented by several harmonic terms together with a 
slowly varying residual (if any), so that we may write: 


f(t) = & cos (at — 71) + C2 cos (at — 72) +... + y(Z). 


If we form the sum 
K+ f(e+™ +f(e+ #)+ pre) 


when a = a, taking as many terms, say N, in the sum 
as the extent of the observations permit, the first 
harmonic will become Ne, cos (a,t — 7,). The other 
harmonics will be in different phases in different terms 
of this sum unless their frequency is a multiple of a, ; 
and so their resultant will in general be small com- 
pared with that of the harmonic a,. The effect of 
any harmonic, say, ag, will not be small if its extreme 
phases in the sum (7) do not differ by more than z 
from one another, for then every term has a positive 
component along the line of the mean phase. This 


will be so whenever a,( 2") differs from a multiple of 
Qy 


2n by less than a/N; that is, when a,/a, differs from 
‘an integer by less than 1/2N. On the other hand, if 
the total range of the observations—that is, 2%N /a,— 
is an integral multiple of the period of any other 
harmonic, such harmonic will absolutely cancel out 
in the sum unless its period is also a sub-multiple of 
a,. The effect of the slowly varying function y(#) 
on the sum considered will only be to add to the 
harmonic Nc, cos (a,t — 7,) another slowly varying 
function which can do little to obscure it. 

If this sum be plotted against ¢, the result is thus 
approximately the sine curve N«¢, cos (a,t — 7,) plus 
something which is approximately constant. ‘The 
amplitude of this harmonic can be found by measuring 
the total variation of height, V, of the graph which 
is equal to 2Ne,; also, if the position of the maximum 
value is at 7= #4, we have y, = ayty. 
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If this curve is not very approximately sine shaped, 
the presence of harmonics whose frequencies are 
multiples of a,/2~ are indicated. If this is so, these 
harmonics can be determined by subjecting this curve 
to the usual process of harmonic analysis. 

The only difficulty about this method of analysis is 
to know when the value we take for 2a/a in forming 
the sum (7) agrees with the period of one of the 
harmonics present in f(2). 

If we use a value of a slightly different from a, the 
term ¢, cos (a,¢ — y,) will not be in the same phase 
in all the terms of (5); and if 27Na,/a were equal to 
(N + 1)2z, the sum of all these terms would vanish, 
since their phases would be uniformly distributed over 
arange of 2m. So it follows that the range of variation, 
V, has a central maximum when a= a,. In fact, 
the variation is clearly according to the now familiar 
euiyey == = ; or if we always consider V as positive 


(and it is not easy to tell when it should be considered 


iy ‘sin x 
as negative) according to the curve y = oe 
We accordingly work out the sum (7) for several 
different values of ¢ for each different value of a taken 
and plot the total range of variation V for changes of 
t against a (or 2/a) and pick out each central maxi- 
mum. When drawing a graph through the points of 
this curve, unless we clearly see that one particular 
point lies right az the peak of one of the central 
maxima, we must re-work the sum (7), for the value 
of a which the graph indicates will be the highest 
point of the peak; for since the peaks are steep (the 
steepness increases with the time over which the 
observations extend), a point only a little way off the 
true centre may give a value of V much below the 
maximum value which is sought. 
Having found these peaks, we can easily write down 
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the harmonic term corresponding to each peak. 
Subtracting the sum of these terms from f(t), we can 
then examine the remainder and see whether it 1s 
constant, or some slowly varying (generally a decaying) 
function; or whether it consists of harmonic terms 
whose periods lay outside the region examined, or 
perhaps were missed, by taking too great a gap between 
the different values of a in the graph for V. 


An IxuustrativE EXAMPLE OF THE SEARCH FOR 
PERIODICITIES. 


We will suppose that some phenomena have been 
observed at 120 equal intervals of time (which we will 
call one day) with the following results : 


Day. f(t). Day. f(t). Day. f(t). Day. f(t). 
fo) 19 31 3 62 24 93 14 
I 25 32 9 63 25 94 24 
2 28 33 17 64 25 95 34 
3 30 34 27 65 25 96 32 
4 30 35 37 66 17 97 29 
5 29 30 36 67 sae) 98 22 
6 21 37 32 68 6 99 14 
7 14 38 25 69 5 100 4 
8 II 39 17 70 5 IOI 4 
9 9 40 7 71 15 102 5 
TO 9 41 7 72 24 103 9 
II 19 42 8 73 30 104 14 
iw 28 43 12 74 34 105 23 
13 34 44 18 75 35 106 24 
TA 37 45 26 76 25 107 2 
15 39 46 26 77 15 108 22 
16 29 47 26 78 6 bgore) 18 
17 18 48 24 79 fe) 110 1092 
18 10 49 21 80 —4 III 13 
19 3 50 15 81 4 TE 13 
20 fe) 51 16 82 14 II3 33 
21 8 52 16 83 24 114 14 
22 18 53 16 84 33 I15 18 
23 28 54 177, 85 39 116 14 
a4 37 55 21 86 33 117 13 
25 43 56 77 87 24 118 13 

26 36 57 16 88 14 11g 13 
27 28 58 16 89 4 120 I2 
28 18 59 16 go —5 
29 7 60 14 gl =I 
30 cal 61 20 92 5 


FOURIER’S THEOREM 169 


all 


Here we see successive maxima at intervals of 10 or 
11 days, so there is certainly a component harmonic 
whose period is about that. We will test first of all 
for periods of 8, 9, 10, 11, 12, 13 days respectively. 

The work for the 8-day period is shown in the table 
below. We have only made use of the observations 
extending over Io periods, since the observations do 
not permit of us using a greater number in the case 
of the longer periods. 


TABLE A. 
ime) 25 28 30 30 29 21 14 
II 9 9 19 18 34 37 39 
2 18 10 3 fo) 8 18 28 
37 43 36 28 18 7 al 3 
2 17 27 37, 36 32 25 17 
9 7 8 12 18 26 26 26 
24 Ze 15 16 16 16 17 20, 
17 16 16 16 14 20 24 25 
25 25 17; Io 6 5 5 15 
24 30 34 35 25 15 6 fo) 
202 211 200 206 I9gI 185 178 188 


eieRance, V, is 201 —9178 = 33. 


Proceeding in this way for the other periods, we get 
the following results : 
Period 8 9 Io II 12 13 
Ne 33 32 269 35 202 Be) le 


No satisfactory graph can be drawn from these 
irregular figures, so we must test similarly for inter- 
Mieaiate petiods, say, of 8-5, 9°5, 10°5, 11-5 and 12-5 
days. ‘To test for 84 days we write down in the first 
bern the value of j/(t) foro, 8, 17,'25, 34, 42, 51, 
59, 68, 76 days respectively, and in the successive 
columns for 1 day greater than each of these numbers 
successively. The intervals between these numbers 


-* In this last case we can only test over a range of 9 periods and we 
get V= 34. We therefore increase the result in the ratio 10:9, 
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are 8 and g days alternately. Had a complete graph — 
of f(z) been available for the whole of the 120 days, 
we could have picked out ordinates at intervals of 
8-5 days for this test, but ordinates at intervals of 
8 and g days alternately are quite as satisfactory, since 
the difference of phase from the mean phase is only 
2/17, the cosine of which is +93. 

Similarly, if we were testing for a period of 84 days, 


Ce See ee 
as 
Peete es. 
i) ae 
Pie Rea 
7 ee 


we should write down in the first column the values 
Gh f(t) forkt = 0,8. 16, 25,133.40. SOsn oe eolee ce 
sively, where the intervals are 8, 8, 9, 8, 8,9 . 

In this manner, we obtain the following results : 

Period 8-5 9°5 10°5 II°5 12°5 

bees 58 112 167 184 116 

Fig. 24 shows the graph obtained by plotting V 
against the period in days. 

There is clearly one period of 10 days very approxim- 
ately and another about 12 days. Since the position 
of this latter maximum is a bit uncertain, we calculate 
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V for a period of 113 days, which gives V = 199. We 
then infer a period of about 11-9 days for this second 
peak. 

The heavy dotted curves show the secondary maxima 
corresponding to the peak for the 1o-day period. 
This makes it clear that the distortion of the peak for 
the 12-day period is largely due to the first secondary 
maximum of the 10-day period peak. We may there- 
fore take 12 days for the second period. We cannot 
possibly expect to find the period to a greater accuracy 
than this; for the difference is only 1 day in the length 
of 10 periods, which is only a difference of one-twelfth 
of a period in the whole range of the observations. 

We next plot the sums of Table A for the 1o0- and 
12-day periods against 1. 

For the 12-day period the sums are: 


ts; Sum. ihe Sum. t. Sum. 
fe) 269 4 18r 8 gi 
I 291 5 138 9 139 
Zz 2606 6 92 be) 176 
3 241 7 89 II 235 


These when plotted give sensibly a sine curve of 
which the equation is 


: 2m t 1 
yy = 100 sin (= -b =) 
2 3 
superposed on a gradually decaying function; but 
since each “sum” is the sum of ten items, the equation 
of the harmonic in /(#) is 
: 2m t a4 
y= ro sin ( sl +t ). 
12 3 
Since the range of observations used in obtaining 
this curve—120 days—is exactly 12 times the 10-day 
period, this curve is not at all affected by this latter 
period. It is different, however, with the test for 
the 10-day period. A range of 100 days only was 
used for this test, which is not a multiple of the 12-day 
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period. Since we have sufficient observations at our | 
disposal, we re-calculate this test, using 12 ten-day 
periods instead of only 10. We then get the following 


figures : 
t. Sum. ts Sum. t: Sum. 
fe) 75 4 310 8 187 
I 133 5 369 9 127 
2 191 6 310 
3 252 7 248 


On plotting these figures, we get the result shown 
in Fig. 25, which we recognise as a Class IV periodic 
function * whose equation is of the form 


aa 
REBREEe S| 


Fic. 25. 


ew {cos a + cos 3% , Cos 5x 
a 9 25 

On determining the amplitude, which is as near 

150 as possible, and dividing by 12, we find that the 


expression 


I0of{. 2nt . I 
12-5 — 10 cos Mrs EGS MSU Ts : cos iz ae sare a 
on Ae wate umes lo 0 5 10 J 


is a constituent of f(2). 


+o} 


: Neglecting the gradually decaying function mixed up with it 
which is indicated by the ordinate for the 1oth day being less than 
the initial ordinate. 
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We next subtract the sum of this term and the 
12-day period term from f(z) to study the residue. 
The values of both terms are very easily read off from 
the graphs we have previously constructed of them. 
We do this for every fifth day, thinking that sufficient. 


The following results are obtained : 


ise Residue. t: Residue. 

fo) 10°3 65 5:0 

5 ere 79 22 
10 g:0 75 59 
i5 9:0 80 4°7 
20 8-7 85 4:0 
25 8-0 90 37 
30 LY. 95 4° 
35 7:0 100 4:0 
40 7°O 105 3:0 
45 6:0 IIo 3°3 
50 6:3 II5 30 
55 6:0 120 3°3 
60 5:3 


Since f(t) was only given in the first place to integers, 
these residues must be considered as uniformly 
decreasing. Plotting them reveals a curve of an 
exponential form which seems to decay to zero ultim- 
ately. The time of falling to 1/e of its initial value, 
which is 10, is seen, from the graph, to be about 100 
days. 


Hence the expression 


1oe-O!4 + ro sin (= ae =) + 12°5 — —— 
a Le, Phe A 
Su 19 J 


represents the given observations to well within the 
accuracy to which they are given: viz., half a unit. 


CHAPATIER: six 
A BRIEF HISTORICAL SURVEY 


AutHoucH the purpose of this book is to impart 
such a knowledge of the subject that the reader shall 
be able to make practical use of his knowledge, a few 
historical facts about the development of the subject 
will probably be welcomed for their interest. “The 
history of our subject is very simple : * so much is due 
to Fourier and so little to other, people. 

The fact that an analytic function could be repre- 
sented in what we now call a Fourier series was known 
in 1777 to Euler (1707-1783), who obtained the 
coefficients in the same manner as we obtained them 
on page 22. With an analytic function, of course, 
there is no question about the possibility of the 
expansion. 

Jean Baptiste Joseph Fourier (1768-1830) was, 
however, the first to assert that any finite artificial 
function could be represented, between any two finite 
limits, by such a series, and that the coefficients were 
given by the same formule as for analytic functions. 
Fourier’s assertions were at first emphatically denied 
by the leading French mathematicians of the day, 
including Lagrange, and his proofs were regarded as 
inconclusive. 

Fourier arrived at his results during the consider- 
ation of the solution of problems on the flow of heat ; 


* We are not now considering the complicated theoretical develop- 
ments of the theory of Fourier’s series which have been made in com- 
paratively recent years and with which this Practical Treatise is 
entirely unconcerned. 
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and his results, though published previously, are fully 
incorporated in his famous Théorie Analytique de La 
Chaleur (1822), one of the great classics of Mathe- 
matical Physics, which even now no mathematician 
ever reads without feeling a genuine admiration for 
Fouriex’s genius, or ever refers to except in terms of 
praise. Fourier’s work and Fourier’s theorem will 
always remain fundamental in the solution of that 
great variety of important physical problems which 
can be formulated in terms of partial differential 
equations. No doubt Fourier’s proofs of his theorem 
leave a modern mathematician rather doubtfully 
convinced ; but Fourier’s work is a striking illustration 
of the fact that the intuitions of a genius are of 
infinitely more value to science than the proofs of a 
pedant. 

Fourier also gave the special application of his 
theorem to odd and to even periodic functions, and, 
further, the result now known as Fourier’s Integral 
Theorem, which is what the ordinary theorem becomes 
when the period is made infinite. In fact, in nearly 
every way, Fourier showed a much clearer compre- 
hension of the subject than many subsequent writers. 
Fourier alone, at first, grasped the fact that he had 
found a method of great importance. Applied to the 
analysis of analytic periodic functions only, the theorem 
was of no importance whatever: it was only an 
awkward way of effecting a transformation that could 
readily be effected by a little algebra, as in Chapter IV. 

Poisson (1781-1840) and Dirichlet (1805-1859) were 
among the first to give a proof of Fourier’s theorem 
by the direct summation of the terms of Fourier’s 
series, Dirichlet’s work being the more satisfactory, 
and rendering any further doubt about the validity 
of Fourier’s expansion, for functions which satisfied 
his simple conditions, impossible. 

Dirichlet proved that for any given value of ¢ the 
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sum of Fourier’s series was f(#), provided ¢ was not a 
point of discontinuity, and equal to 3{f(t— 0) + 
f(t + 0)} at such a point. Dirichlet’s proof is repro- 
duced in a simplified and abbreviated form on pp.136-8. 
The fact that, at a point of discontinuity of a func- 
tion, the sum of the series was really indeterminate 
within certain definite limits, and that the range of 
such indeterminateness was 1-179 times the abrupt 
change in magnitude of the function, seems only to 
have been noticed comparatively recently ; though for 
ourselves we can hardly make out how it could have 
escaped detection in Fourier’s day that the series 


SINE ~ Sin 2F 


ae coe ea 
I Ze 


X sin x 


when ¢=o0, can have any value that \ dx can 


fe) 
take ; since everybody knew even then that the limit 
when ¢ = 0 of 


Hf) + fle) + f34) +... +f} 
x 
was | fede, and one has only to multiply both 


numerator and denominator of each term of the series 
by ¢ to make the result instantly obvious. 

The subject of the harmonic analysis of a function 
defined by a set of points was not developed in any 
practical detail by these early mathematicians because 
the scientific observations requiring such treatment 
were not then extant. Lagrange (1736-1813), sub- 
sequent to his doubts over Fourier’s results, obtained 
the equation of a finite sine series passing through the 
tops of a number of equidistant ordinates; while 
Poisson extended this result to a finite series including 
both sines and cosines and so gave the fundamental 
formule for all harmonic analysis with equi-distant 
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ordinates. Lagrange’s and Poisson’s results are exactly 
the same as those given in our Chapter V. No doubt 
Poisson, had he required to undertake practical 
harmonic analysis, would have invented Schedules 
for calculating the coefficients almost indistinguishable 
from those we have given. Several formule for inter- 
polation for periodic functions, when the ordinates 
were not equi-distant, were given by Gauss (1777- 
1855). As far as we know, the complete list of these 
formule, with indications of the types of function 
to which they are respectively applicable, is now given 
for the first time at the end of Chapter V. 

One would think it impossible for mathematicians 
to pay much attention to the harmonic analysis of 
functions defined by a set of equi-distant ordinates 
without soon arriving at the expressions giving the 
coefficients of the different harmonics in terms of the 
actual Fourier constants of the function. It is some 
surprise to us, therefore, that we have not seen these 
relations in print earlier than A. Clayton’s article in 
the Fournal of the Institution of Electrical Engineers, 
Vol. LIX, 1920-21, p. 491 ; though it seems inevitable 
that several workers must have obtained these relations 
for themselves prior to this. 

Coming to the subject matter of Chapter VII, 
which is essentially due to Fourier, the late Lord 
Rayleigh (1842-1919) was the first to prove that if 


{@ Lear 


(o.6) 
and \ fd) sin at dt =B, 
104. 


then \ fl)2dt = =| (A? + B2)da, 
en) e) 


It 


which is the extension to pulses, of the well-known 
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result that the energy of any periodic wave is the sum 
of the energy of its different harmonics. 

Finally, we must state that the very simple method 
given at the end of Chapter VIII for the practical 
determination of harmonic constituents in a non- 
periodic function is the same as that given in Whittaker 
and Robinson’s Calculus of Observations (1924). 
Although much attention has been given to this 
important practical problem, from the time of La- 
grange to the present day, the theory of the problem 
that we have given shows that nothing simpler can 
be obtained than this; while at the same time it shows 
that more elaborate methods are not worth the extra 
labour of applying, and are no more necessary than 
is the ideal diffraction grating, which only gives 
a single order spectrum, to the experimenter in 
spectroscopy. 
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